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ABSTRACT 


A design of a highly programmable, land based, real time synthetic 
aperture radar (SAR) processor requiring a processed pixel rate of 
2.75 MHz or more in a ^our-look system was conducted. Variations 
in range and azimuth compression, number of looks, range swath, range 
migration and SAR mode were specified. A number of alternative range and 
azimuth processing algorithms were examined and analyzed in conjunction 
with projected integrated circuit, digital architecture, and software 
technologies. 

The selected design for the Advanced Digital SAR Processor (ADSP) 
employs an FFT convolver algorithm for both range and azimuth processing 
in a parallel architecture configuration. This overall design approach 
met all of the system implementation and performance criteria for 
programmability, modularity, adaptability to VLSI, low risk, reliability 
and cost. 

The report provides algorithm performance comparisons, a detail design 
of the selected system, implementation tradeoffs and the results of a 
supporting survey of integrated circuit and digital architecture technologies. 
Cost tradeoffs and projections with alternate implementation plans are 
presented. 
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1.0 INTRODUCTION 


A general practice In synthetic aperture radar (SAR) systens has been 
to use a specially tailored processor to convert the raw radar signals 
to Imagery. The advent of practical digital signal processing tech- 
niques. brought about by LSI technology, has made high performance, 
multi-mission programmable SAR digital signal processors feasible. 

This Is the final report of a study to quantify various Issues related 
to the hardware development of such a SAR processor. The generic name 
for the processor Is the Advanced Digital SAR Processor and It shall be 
referred to from this point on In the report as the ADSP. 

A very bn ad set of performance goals have been set by the Jet Propul- 
sion Laboratory (JPL) for the ADSP. These are discussed In Section 2.1 
together with the system performance requirements. 

Two principle design facets are examined, processing algorltfms and 
technology. The algorithms are further divided into time and frequency 
domain techniques; while the broad Issue of technology Includes digital 
integrated circuits, digital architecture and software. Section 2.2 
gives a comprehensive description of the candidate processing algorithms 
together with hardware implanent.atlon approaches for them. 

A key factor in the selection of a processing algorithm is the quality 
of imagery it produces. Extensive computer simulations have been de- 
veloped and run to evaluate the performance of candidate algorithms. 

Thii data, given in Section 1.3, was then applied as appropriate to the 
hardware sizing of the various implenentations. 

Technology is addressed in Section 2.4 with the current state of the 
art and projections provided for digital integrated circuits, archi- 
tecture and software. 

A design recommendation for the ADSP is synthesized and described in 
Section 2.5. It meets all of the performance and implementation goals 
of the ADSP. 

Key issues in the selection of an algorithm for the ADSP were its 
effectiveness in compensating for range migration, progammabil ity for 
multiple modes, performance level and the inherent computation 
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r^v-ili^efnents and memory storage. The selection of an algorithm was 
closely tied into the technology issues, particularly digital inte- 
grated circuits and architecture. The design selected minimizes 
memory storage and offers simple memory management techniques. It 
provides a progratranable architecture in modular form which offers 
both incremental development and growth potential with little waste 
of effort. Finally, the design has inherently high reliability 
and maintainability features. 

Cost factors, alternative development plans and tradeoffs are pre- 
sented in Section 2.6 and 2.7. 

Software developed under the program is provided in the appendix. 


ORIGINAL PAGE 

2.0 TECHNICAL DISCUSSION OF POOR QUALimT 

2.1 ADSP SYSTEM REQUIREMENTS 

A satellite SAR system Is depicted In Figure 1. As the radar moves 
along Its flight path, multiple radar pulses are coherently processed 
to achieve a resolution corresponding to the length of a synthetic 
antenna. This length Is determined by the footprint of the real an- 
tenna on the surface. The limiting resolution of the SAR Is thus 
determined by the size of the real antenna and is equal to one half 
its diameter. To achieve radar returns which depict the contours of 
the surface, an oblique Incidence angle is employed. In addition, a 
forward look angle off broadside Is commcn. Resolution In the cross 
track, or range, dimension Is achieved by using conventional radar 
pulse compression techniques. The total range swath covered Is also 
constrained to the dimensions of the radar footprint. 
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2.1.1 Functional Requirements 

The ADSP requiranents have been established In anticipation of pro- 
cessing a wide variety of SAR data. SAR missions which represent 
a comprehensive set for design parameter purposes are given In Table 1. 
Both satellite and aircraft SAR are represented. 

Table 1 

SAR Mission Set 


MISSION SET 

RES^IjTION 

SWATH 

WIDTH 

Ikml 


INCIDENT 

ANGU 

FREQUENCY 

POLARI- 

ZATION* 

AiTnyoE 

SEASAT 

25 

100 

4 

22.54 

1275 

S 

800 

SIR REFIIGHT 

40 


7 

20-70 

1275 

S 

210-280 

AIRCRAFT 

L 

10-20 

■BB 

6 

0-40 

1275 

S 

3-12 

X 

10-20 


8 

o-eo 

9400 

S 

3-12 

VOIR 

HI-RES 52® 

75 

TOD 

4 

52 

1275 

S 

250 

LO-RES 52® 


30 

30 

52 

1275 

S 

250 

ICEX 

150 

m 

6 

3&42 


S 

700 

ERSAR f 

15 

i 

40 

4 

eo® 

1275 

D 

300 

ERSAR X** 

15 

40 

4 

fiOO 

1 

miQiiiiim 

D 

300 


• S • SIMGIE, 0 • DUAL 
•• PREUMINARY 


The most stressing signal processing requirement of t^e missions in 
Table 1 is the ERSAR (Earth Resources SAR) which has been selected by 
JPL for the baseline design mission for the ADSP. Functional re- 
quirements for the baseline mission are given in Table 2. 

Both range and azimuth processing requirements are included. Speci- 
fications are also provided for range migration correction, reference 
function update; integrated sidelobe level and output data format. 
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Table 2 

AOSP Baseline Functional Requirenents 


RANGE CORREIATOR 


INPUT SAMPLE RATE iCOMPLEX WORDS) 

IIMHZ 

RANGE COMPRESSION RATIO 

410 

RANGE PULSE LENGTH 

47S (COMPLEX SAMPUS) 

RANGE INTEGRATED SIDaOBE RATIO DUE TO PROCESSING 

'ISdB 

AZIMUTH CORREIATOR 


SWATH WIDTH (OUTPUT PIXaSI 

4.000 

RADAR PRF 

Z.SOOHi 

NUMBER OF AZIMUTH LOOKS (CONCURRENTLY PROCESSED) 

4 

AZIMUTH BANDWIDTH PER LOOK 

500Hz 

AZIMUTH COMPRESSION RATIO PER LOOK 

350 

RANGE WALK COMPENSATION ACROSS 1 LOOK 

Z8 SLANT RANGE BINS 

RANGE WALK COMPENSATION ACROSS 4 LOOKS 

86 SLANT RANGE BINS 

RANGE CURVATURE COMPENSATION ACROSS 4 LOOKS 

7 SLANT RANGE BINS 

RERRENCE FUNCTION UPDATING RATE 

CROSS-TRACK EVERY 8 
SAMPLES: jaONG TRACK 
EVERY 12.000 RADAR PULSES 

AZIMUTH INTEGRATED SIDaOBE RATIO DUE 
TO PROCESSING 

-ZOdB 

OUTPUT PiXa FORMAT 

(BOTH OPTIONS A AND B ARE REQUIRED) 

A) SELECTABLE 12 BITS IN 
AMPLITUDE OVER 80dB 
DYNAMIC RANGE 


IZBITS^O 


Changes In operating mode are anticipated to occur within a SAR mission and 
these changes could Involve, altitude, look angle, resolution, number 
of looks, swath width and the range and/or azimuth weighting functions. 

The processor Is required to operate In a continuous or burst mode. 

The former has the sensor emitting non-changing parameters. The latter 
Interleaves alternate bursts with different- parameters, yet permits 
contiguous Imagery formation. Alternate bursts may Involve different 
radar frequencies, different polarizations or different ctntenna pointing 
angl es. 

Design objectives of the ADSP Include programmable features 
to permit range compression from 10 to 600, azimuth compression from 
20 to 350 per look and swath width from 500 to 4000 pixels. A desired 
feature provides the capability to tradeoff looks, resolution 
and swath width within the capabilities of the processor. 



Modularity In design was a key functional design objective. This 
feature could provide greater swath width (up to four times the base- 
line), more parallel looks and multiple frequencies and polarizations. 

2.1.2 ADS P System Scope 

The principal elements of the ADSP as shown In Figure 2 are an Input 
data handling facility to permit either real time or recorded data 
to be processed, a range correlator, an azimuth correlator and a multi- 
look Integrator. All of the high speed processing hardware Is encom- 
passed In the pipeline from data Input to output. The host computer Is 
used to set up conditions for mission operation. These will Include 
SAR parameters and overall control commands. The hardware Itself 
will respond to the program commands and under firmware 



Figure 2 
ADSP Elements 


control will complete the required operations. A goal of the system 
design was to maintain software simplicity for the host computer by 
Incorporating as much as possible of the SAR variables In the hardware. 
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2.1.3 Analysis of ADSP Requirenents 
2. 1.3.1 SAR Parameter Variation 

The azimuth resolution for operation near broadside Is given by: 

®A2 “ 2V T- 
s L 

where Rg ■ slant range. A" wavelength, Vg ■ satellite velocity and 
Tl = coherent look time. The azimuth beam width factor Kg, Is asso- 
ciated with the receive weighting required to achieve a desired side- 
lobe level. For -36 dB sidelobe level. Kg "1.2 for a Taylor fi ■ 8 
weighting function. This weighting produces an Integrated sidelobe 
level, ISL, of -25 dB. When other sources of sidelobes are Included 
such as digital granularity and the lack of Interpolation In range 
cell migration compensation, an overall ISL of -20 dB Is obtained. 

The chirp bandwidth associated with the coherent look time Is given by: 
given by: 

2Vs^Tl KaVs 


For the ERSAR case Vg » 7.73 km/sec, <5^^ ® ^ 

B|^ = 618.4 Hz for each look. At X = 0.235 m, Tj^ = 0.686 sec., 

®L^L “ angle compression during the look angle = 

424/1.2 = 353. 

For four looks the total processed bandwidth => 4Bj^ ■ 2473.6 Hz. 

The azimuth beamwidth that just matches this spectrum Is given by 
^B “ ^'^s (2.15°). The antenna length required 

for this beamwidth Is = .886 A/ = 5.55 meters, assuming uni- 
form Illumination. The PRF of the ERSAR = 2500. The ratio of the PRF 
to the four-look bandwidth = 1.01. This ratio Is lower than It should 
be for low doppler aliasing (angle grating lobes). A calculation shows 
that the Integrated grating lobes from the four looks Is only -8.2 dB. 

A higher ratio of PRF to four-look bandwidth would be required for 
reduced grating lobes; however, the range grating lobes would then be 
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required for reduced grating lobes; however, the range grating lobes 
would then be excessive unless the swath were reduced by the same 
ratio that the PRF was increased. The sampled range to allow for a 
40 Ion swath at 60° incidence and a 35.2 jisec transmitted pulse requires 
266.1 |isec. The ratio of the 2500 PRF period to 266.1 jjsec ■ 1.503 
is a reasonable value to inhibit range grating lobes adequately. A 
lower ratio of periods caused by a higher PRF would cause the range 
grating lobes to appear. 

SAR mission sets given in Table 1 can be evaluated as a function of 
the umber of looks, swath width and required resolution for the 
optimum PRF and antenna parameters. For the case of dual frequency 
such as ERSAR-L and ERSAR-X, the real antenna lengths, PRFs, and look 
bandwidths would be identical. However, the antenna widths and angle 
compression should be scaled by the wavelength. Dual polarization 
reception could be simultaneous, if desired. However, dual polari- 
zation transmission would have to be orthogonal either by frequency 
or time. The possibility for change in look angle depends upon an 
adequate ratio of PRF period to A/D sampling time for the desired 
swath. Electronic scanning or beain switching can be employed to change 
the look angle if the ratio is less than 1.5 The requirement for a 
"burst mode" involving interleaving different frequencies, polarization, 
and pointing angle may also require azimuthal electronic scan. This 
is discussed in Section 2. 1.3. 2. Of course, several PRFs may be 
required to center the swath window between the ambiguous transmission 
pulses to avoid eclipsing if the look angle or orbit attitude is 
changed. 

2. 1.3. 2 Mode Variations 

The ADSP is required to have the flexibility to process data from a 
wide variety of different missions described in Table 1 and to accom- 
modate continuous or burst modes. 

Figure 3 shows the case for four looks in the continuous mode integration 
for azimuth cell j, where Look 4 coherently integrates the first 
samples and Look 1 coherently integrates the last samples over 
the azimuth beam. The four looks are noncoherently integrated and 
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registered at the correct range and azimuth cell. As each new azimuth 
sample Is received, a beam Is completed for each look and one final 
Integrated multllook azimuth Image sample Is completed. 

ORIGINAL PAGE 
OF POOR QUALITY 





Figure 3 

Continuous Mode Integration Requirements 

In the burst mode, a number of variations are possible. When bursts 
are Interleaved and do not Involve frequency chanyes, the Integration 
patterns shown In Figure 4 are obtained. In the case shown in Figure 4a, 
two polarizations are Interleaved by bursts in a four-look mode. The result 
is that the output Imagery is available on a two-look two-polarization 
basis. A polarization comparison of Image samples can then be made. 

The processing can be Identical In form to the continuous case except 
that Inefficiencies occur In that not all of the antenna coverage Is 
Integrated and used In the final Imagery. This is because the field 
must be Illuminated for a full look before the data can be processed 
to obtain a full integrated signal-to-nolse ratio with the desired resolution. 
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A large number of looks can be transmitted on an Interrupted basis as 
shown In Figure 4b. In this case, 32 looks are processed. Look 
Integration procedures may be varied as a function of time. If the 
burst sequences are short. In the limiting case, each look can be 
an unfocussed segment of the entire beam. The beam Is then split In 
angle Into the number of "looks" processed. The Multiple burst looks 
must be placed In proper range registration prior to Integration. 


The conditions change If bursts of different frequencies are transmitted. 

If we assume that the antenna size remains the same, the high frequency 
has a proportionally narrower beam. Thus, if full coverage Is desired, 
multiple bursts at different pointing angles must be transmitted. This Is 
Illustrated In Figure 5. For simplicity, the example of high frequency is 
four times the lower frequency and four looks are employed. The first 
burst at frequency Fi provides the four-look coverage Indicated. Quadrupling 
the frequency narrows the antenna beam proportionally and also Increases the 
ground resolution by the same factor. To obtain the same final Image 
resolution, only one fourth as many pulses need be transmitted per burst. 

But since the beam is narrower, it must be steered to a total of four 
angles with a separate burst transmitted at each angle. Four looks are 
processed at the high frequency also to maintain the final image resolution. 

A strip map can be generated in this way which has one look at each 
frequency for each image location. 


UPJ 12PI 
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Figure 4 Interleaved Polarization Burst Mode 
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Figure 5 

Interleaved Bursts At Two Frequencies 


2. 1.3. 3. Range Cell Migration 

Range cell migration occurs during the synthetic aperture integration 
due to the radial range change between the radar and a ground scatterer. 
This migration consists of a 1 Inear^component with time, rt, and a 
quadrature ratio component, 1/2 r‘t. The range cell migration function 
Is analogous to the phase focusing function and can be derived In 
the same manner. The time, t. Is generally referenced to the point 
of closest approach, r^^, for a non-rotating earth. A rotating earth 
component of velocity toward the radar will actually make the point 
of closest approach occur slightly later In time. An orbital altitude 
rate associated with eccentricity will also affect the time of closest 
approach. The radar beam pointing departure from broadside will shift 
the time of signal return and thus Introduce an r term equal to cos0v« 

where Is the angle between the velocity vector and the ground 
scatterer. 

Table 2 Indicates that for ERSAR there are 86 slant range bins of 
range walk to be compensated for across tSe four looks. Assuming 
an Incidence angle of 60°, a ground bln of 15 meters results In a slant 
range bln of 15 sin 60° =12.99 m. Thus, the range walk Is 86 x 12.99 = 
1117 m or 7.45 ;jsec. The range curvature Is given as 7 slant range 
bins = 90.9 m. The maximum range v/alk across the farthest look = 1117/4 + 
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3/4 X 90.9 ° 347.4 m or 26.7 slant range bins (slightly less than the 

28 specified).. The 4«1ook time Is 4 x 0.686 ■ 2.744 sec. Thus, the 

maximum range rate Is assumed to be 407.1 m/sec. The corresponding 

2 

acceleration Is 96.6 m/sec . The coherent Integration of each look 
must take place along the curved range cells to achieve the desired 
resolution. The range cell migration Infers a two-dimensional matched 
filter requirement for which the processing algorltiwi alternatives 
In Section 2.2 have been developed. 

2. 1.3. 4. Clutter Lock Requirements 

The approximate satellite position, velocity and altitude Information 
Is provided to the SAR processor. The SAR antenna pointing (roll, 
pitch and yaw) will be controlled within + 1^. Uncerta Ini ties In the 
data could easily cause the signal to depart In center frequency by 
more than the per look bandwidth. For example a yaw angle of 1*^ 
translates to a cross-cone angle of the slant range vector of 0.83° 
from broadside. This results In a shift In doppler of 950 Hz for the 
L-band 300 km ERSAR satellite. At X-band the same yaw shifts the doppler 
by 7153 Hz. Earth velocity for a satellite polar orbit has a radial 
component at the equator of V^q sin (incidence angle) ‘ 463 x sin 60° 

= 401 m/sec, causing a doppler shift of 3413 Hz at L-band. For 
accurate ERSAR mapping such as 10 m absolute location (AX), the 

doppler shift due to the earth must be known to, 

2V Ax 2 X 7730 x 10 

Af = -irl * = Hz cr 0.034% 

564.4 X .235 

Clutter locking should be maintained to within 100 Hz to keep the 
per look tuning within 20% of the look bandwidth (4% of the PRF). 

Errors greater than 50% of the look bandwidth will degrade the azimuth 
sidelobes due to azimuth spectral aliasing. Clutter locking errors 
due to Satellite yawing do not cause a mapping error since the doppler 
used to tune the processing bandwidth to the signal spectrum is em- 
ployed in forming the map location. 

2.L3.5 Automatic Focusing Requirement 

The doppler chirp rate f^ of a target viewed by the SAR beam must be 
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matched by the azimuth correlation process to achieve the azimuth 
beam compression and desired target resolution. It Is equal to ■ 

2r /A • The radial acceleration r Is a function of the satellite 
velocity Vg, the slant range Rg. the ratio of the target distance y 
from the orbit axis to the orbit radius, Rg + h^, for near broadside 
operation for a non-rotating earth. The relation for a circular orbit 
Is 

Vg^y Vg^ Rg cos 4 

** Rg (Rg + hg) Rg (Rg + hg) 

where a Is the earth angle between the satellite and the earth target and 

hg Is the orbit height; Rg Is the earth radius. The nominal value of 

r » 100.84 m/sec^ at the midpoint of the ERSAR swath (Vg ■ 7.73 km/sec, 

R^ » 6373.95 km, he » 300 km, R„ » 564.4 km, ■ 4.2°). Over the 40km 

swath, r varies from 104.03 m/sec at the near edge to 97.8 m/sec at 

2 

the far edge. The variation Ar » 6.23 m/sec occurs across the 40 km 
swath, resulting In a change In doppler rate 


^ ^ . 53 Hz/sec 

^ X .235 

The depth of focus, expressed In terms of allowable change In doppler 
rate, can be determined assuming a quadratic phase error allowance of 
+ ir/2 at both ends of the synthetic aperture relative to the center 
of the aperture. The phase error 

^ T/2 T/2 . y dY • 2 

-2ir j- I Af^ dt'^ . J^Afd Tl 

0 0 

The term Af^j * t Thus, the depth of focus, DOF * (4/Tg^), 

where T^ 1s the time to form the synthetic aperture for each look. 

For an azimuth compression ratio of 350 and a bandwidth of 500 Hz, 

Tl = 0.7 seconds. Thus, the depth of focus * 4/. 7^) = 8 Hz/sec. From 
the point of view of a single look, the reference function does not 
have to change more often than 40 km (8 Hz/sec 7 53 Hz/sec) = 6.04 km 
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swath. However, the reference function wilt have to change more often 
In a cross-track direction to permit the registration of the four looks 
within 0.2 azimuth cell . 

The presence of a quadratic phase error corresponds to a mismatch of the 
reference function chirp rate with the target chirp rate. This mismatch 
has the effect of a doppler error between the four looks, 4f, producing 
a misregistration of the linages, Ax. 


Af • 


Ar. 


For X ■ 0.2, 9f^2 * ^ 15 


2 X 7730 X 3 

3 meters, Af 

.235 X 564.4 km 


0.35 Hz 


The time Interval between the four looks ■ 3Tj^» 2.1 sec. The allowable 
change In doppler rate cross-range ■ (0.35/2.1) »0.17 Hz/sec, Thus, 
for registration, a new reference function should be applied more often 
than 40 km x (0.17 Hz/sec ~ 53 Hz/sec) » 126 meters of swath, corre- 
sponding to 8.4 resolution elanents. The specification calls for a 
new reference function cross-track every 8 samples, which Is reasonable. 
Fortunately, the chirp rate variation across the swath Is deterministic 
from the geometrical parameter variations and may be computed from 
the one or two measured chirp rates refined by auto-focusing 
techniques. 


I 
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2.2 PROCESSING ALGORITHMS 


ORIGINAL PAGE IS 
OF POOR QUALITY 


This section describes the candidate SAR processing algorithm 
considered for the ADSP and gives pertinent Implementation data on 
than. Performance levels are treated separately In the following 
section (Section 2.3). 

2.2.1 Basic SAR Processing Functions 

The formation of Images from synthetic radar arrays follows the same 
principles found with optical lenses or real radar antennas. Images 
can be focused In the far-fleld when the energy reaches the antenna 
as a virtual plane wave relative to the antenna diameter or In the 
near field. The far-fleld SAR Is generally ten'ed an unfocused .SAR 
and, being relatively simple to process. Is not treated In this 
study. Processing In a focused SAR is essentially a two dimensional 
matched filter process as Indicated in Figure 6. 



Figure 6 

Focused SAR Image Processing Functions 
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The range, or cross-track, dimension may be processed by conven- 
tional rcdar pulse compression techniques, the method used depending 
on the type of radar waveform anployed. The azimuth, or along track, 
dimension can be considered In the same functional sense as the 
range dimension. Here, the radar footprint of the real antenna with 
Its doppler frequency response forms a quadratic phase (or linear FM) 
function across Its width. As the real antenna moves with Its scan- 
ning platform the radar footprint "waveform" moves along the Image 
plane In azimuth In the same manner that the uncompressed radar 
pulse moves across the Image plane In range. Thus a filter matched 
to the radar footprint "waveform" operating In the azimuth dimension 
will compress the radar footprint to the resolution defined by Its 
time and bandwidth parameters. The limiting azimuth resolution for 
a focused SAR Is one half the azimuth dimension of the real antenna. 



Figure 7 

Satellite Radar Swath 

Two other functional requirements are important In SAR processing, 
range migration and multiple looks. As a SAR satellite platform 
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moves along Its track, the earth velocity causes the radar return to 
move linearly with the motion of the spacecraft (S/C) (See Figure 7). 
In addition, the range of a fixed point on the earth will vary In a 
quadratic manner as the S/C radar footprint moves past. 

Radar returns from objects vary In amplitude as a function of the 
aspect angle of the radar. This tends to give radar mapping return 
a speckled appearance. The effect can be reduced by non-coherently 
summing coh«!rently processed segments of the radar footprint "waveform 
This multi-look processing reduces the final Image resolution In 
proportion to the number of looks processed. The bandwidth of the 
processing required for each look Is also proportionally reduced. 

The combined effect of range migration, linear and quadratic range 
walk and multiple look processing Is Illustrated In Figure 8. 



Figure 8 

Azimuth Resolution Element Motion 
Relative Fixed Antenna 
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2.2.2 Time Domain Processing 

2. 2. 2.1 General Time Domain Processing Consideration 

The general form of a time domain processor is shown in Figure 9. 

It essentially implements the SAR processing equations with a time 
domain analog of the functions. After range correlation, which 
consists of a tapped delay line finite impulse response filter in a 
time domain implementation, the signal is filtered (presummed) into 
bands corresponding to the looks to be processed. This serves to 
greatly reduce the computation requirements of the azimuth correlators. 
A simple example will illustrate the effect. Consider an azimuth 
correlator with 1000 azimuth samples over the total beamuidth with 
an azimuth doppler bandwidth (and sample rate) of 1000 Hz. The 
number of complex multiples per second to perform a time domain 
correlation is then 10^ per range cell. If the signal is processed 
as four looks, the bandwidth per look becomes 1000/4 and the number 
of azimuth samples covering each look is 1000/16 giving a total 
computation rate for four parallel look rilters of ( 1000 X 1000 X 4) » 
10®/16 per range cell. Thus, the computation rate is ^inverTIly 
proportional to the square of the number of looks processed. The 
net complex computation rate for a time domain correlator is 


INPUT 


LOOK 1 



OUTPUT 


Figure 9 

Multi-Look Time Domain Processor Functions 
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where; is the synthetic aperture integration time, is the 
doppler bandwidth, is the number of looks, and ot is the over- 
sampling factor. 

After prefiltering and azimuth correlation, the multiple looks are 
summed. A time delay must be inserted in the look filter prior to 
summation to compensate for the time offset of the look energy from 
particular looks. 

Two general types of architecture have been identified for time domain 
processor and have been evaluated in Reference 1. They have been de- 
signated Type A and Type B and we will continue to use those designations 
for clarity. Type A (Figure 10) is a classical tapped delay line 
matched filter in which the intertap delay lines connected serially 
contain the entire SAR range swath. 


INPUT 



Figure 10 

Serial Time Domain Azimuth Processor 
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Type B (Figure 11) performs azimuth filtering by the construction 
of parallel azimuth beam correlators, each accumulating the processed 
beam over the entire range swath. The general simplicity of form for 
these processors plus their minimization of memory storage makes them 
candidates for the ADSP. We will develop sufficient detail on these 
processors to permit evaluation on an equal hardware technology basis 
with the frequency domain methods. 



Figure 11 

Parallel Time Domain Azimuth Processor (Type B) 

2. 2. 2. 2 Serial Time Domain Processor (Type A) 

Figure 12 Is a more detailed sketch of a Type A or serial processor. 
The processing is segmented Into processing units labeled as corre- 
lator #1, correlator #2 ... up to correlator N. Each processing 
unit has a FIFO (first-in, first-out) memory large enough to hold a 
single pulse radar sweep. Since the range processing Is completed. 
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the FIFO memory covers the range swath plus any uncompensated range 
migration. The FIFO memory can be implemented with shift registers; 
but random access memories which have more capacity per chip, consume 
less power and are less expensive are a preferred choice. In addition, 
the storage length of a FIFO random access memory is easily changed 
by address control. The address control mechanism consists of a 
counter, a reset control word and a comparator. On each count new 
data is read into the memory and the data stored is read out. When 
the count reaches the value of the reset control word, the comparator 
will output a signal which resets the counter to zero. Thus, the 
delay (FIFO) capacity is easily changed. The memory read -mod ify-write 
cycle time must be less than the clock period. If it is not, multiple 
memory units can be multiplexed to meet the speed. 



• OUTPUT TO 
MULTIIOOK 
INTEGRATOR 


Figure 12 

Detail of Type A Processor 
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Range migration compensation Is accomplished within each correlator 
by a combination of variable clock interval (range cell) shifts (Tv) 
plus Interpolation. Simulations Indicate that a two point interpolator 
will give adequate performance, but for the highest quality Imagery 
It may be desirable to use a three or four point Interpolator. The 
focusing function storage In each correlator unit Is related to the 
required Increments of range migration compensation. Each significant 
variation In focusing function, nominally every 8 range cells for the 
baseline SAR mission, has an associated range migration compensation. 
The total reference storage In each correlator unit is at least 4000/8 
- 500 words of at least 34 bits. These bits would be assigned as 22 
bits for focusing function, 7 bits for Incremental range migration, 
and 5 bits for the Interpolator. 

In addition to the variable time delay In each range migration compen- 
sation memory, a fixed Incremental delay Is Inserted to aid In con- 
veniently obtaining an output sunmiatlon. In effect, each correlator 
segment Is operated offset In time by one range cell. Then, as the 
partial sums are passed from one unit to the next through register 
delays they are automatically added with the appropriate delays. 

In the baseline SAR system requirements, the four look prefiltering 
will give a clock rate of 11/4 ■ 2.75 MHz for each azimuth correlator. 
Obviously, the computational hardware In each correlator unit can 
be operated at a higher rate and this Is reflected in the three circuit 
tabulations given in Table 3 for the Type A system. A typical con- 
figuration fcr a 4X multiplexed correlator is shown in Figure 13. 

Note also In Table 3 that the number of Interconnections Is modest 
for the Type A correlator modules. 
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Figure 13 

Multiplexed (4X) Type A Correlator Channels 


Table 3 

Type-A Circuits - Interconnections 


Ewmai 

CIRCUIT TYPES 

U&l 

2X 

S.5MHI 

4X 

llWMl 

PRI STORE (8K X 221 

4RX4RAMS 

12 

24 

48 

PRI MEMORY CONTROL 

MISaUANEOUS 

ID 

10 

10 

RANGE COUNTER 

COUNTERS AND REGISTERS 

4 

4 

4 

RMC AND CONTROL 

MISIXUAMOUS 

60 

70 

90 

FOCUS MULT 

4 MULT. 2 ADDS 

IS 

IT 

19 

ADDER 

ADDERS AND REGISTERS 

16 

24 

26 

CONTROL INTERFACE 

MISCELLANEOUS 

JO 

J2 

14 


TOTALS 

12T 

161 

211 


I/O INPUT 22 

OUTPUT OEIAVID S 

OUTPUT SUM 30 

CONTROL IN-OUT _H. 

TOTAL 106 
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2. 2. 2. 3 Parallel Time Domain Processor (Type B) 

A detailed diagram of the Type B time domain processor is shown in 
Figure 14. Because of the nature of the architecture, a single 
common Interpolator unit can be used for all of the correlators. 
However, this step creates a large I/O pin requirement and the 
accuracy that can be used is therefore limited. The example shown 
has eight, two point interpolators feeding the input bus. Range 
migration compensation is accomplished for each correlator by se- 
lecting the appropriate interpolated point and gating the desired 
range cell. The entire interpolation control sequence is generated 
and/or stored in the Common Control Unit. As the first correlator 
starts its integration process, the required RCM controls as a 
function of range for the first radar PRI are applied. This control 
sequence is then passed to the second correlator and the RCM 
control sequence for the second PRI is applied to the first corre- 
lator. This process is repeated until the last PRI to be integrated 
by the first correlator is received. A similar process is used for 
the focusing function store. The reference words pass serially 
through the string of correlators. Wnen a correlator unit completes 
its integration process, its output is selected and it starts the 
integration process again. 
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Table 4 lists the circuit counts for various multiplexing alter* 
natives for the Type B architecture. 

Table 4 

Typ6«B Circuits » Interconnsctions 





CIRCUI^IMODUli 

4X 

BMiai 

CMMITJXaS 


llMHt 

SQECT 

SilSW 

22 

s 

22 

RCM MEM AND can 

MISCEUANEOUS 

14 

28 

96 

FUNCTIOI STORE 

MBWORY 

6 

12 

21 

ADDRESS GEN 

MISCEUANEOUS 

10 

18 

26 

RANGE raiACCUM 

4KX ARAMS 

M 

32 

64 

COMRiEXMULT 

4 MULT 

B 

IT 

19 

ADDER 

ADDERS 

14 

18 

a 

CONTROL INTBIFAiX 

MISCEUANEOUS 

JO 

JO 

JO 


TOTAL CIRCUITS 

IfO INPUT 

FOCUS IN*OUT 
INTERPOIATE IN*OUT 
OUTPUT 
CONTROL 

109 

1T6 

48 

8 

» 

16 

161 

23T 


TOTAL m 


2. 2. 2. 4 Type A, B Comparison 

A qualitative comparison of the Type A and Type B time domain 
processors Is shown In Table 5. There Is not a significant difference 
In the estimated cost or performance of the two approaches. Type B 
offers more flexibility for variation of SAR parameters. However, 
this Is only true If the number of looks are held constant and the 
compression factor Is less than the established size. For both 
systans, if the number of looks are reduced it takes significantly 
more hardware correlator modules to achieve real time processing 
rates and the system must be reconfigured. Our preference is with 
the Type B system because of Its greater prograinnabll ity and 
somewhat less fault sensitivity. 


26 


^MMWAL PA'JE 
OF WOR QUALITY 


Typ««A, B Comparltoira 


cwmciaisTie 


TYPE A 


mil 


COST. SIZE 

MOOUARIIV 
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• PARAUaaMREIATOliS 
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HIOtOaFT C0IWH.VQI) 
PARAUELOMREIATORS 

an. 


Hlffil IN SERIAL CCN1R0IS 
SERIAL CONTROLS 


PR06RAMMA0ILITV 
• AZIMUTH COMPRESSION 


• RANQE SWATH 

• NUMBBI OF LOOKS 


• BURST MOOES 

• RQ.IABIUTV 

• IMPUEMOITATION 

PROSIEMS 


POOR-lFGREATBi 

G000~MEMmVC(MTR0L 
POOR-WEOEjmNSIVE _ 

WITH FEWER LOOKS 
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6000 -MEMORY CONTROL 
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ALLY 


WITH FEWER LOOKS 


IN ACCORDANCE WITH AZIMUTH COM>RESSION VARIATIONS 
fPSll.OyE TO SERIAL DEIAYS FAIR OIK TO SERIAL CONTROL 
AND SUM PIKUliS 

CONTROL BUS 

- IF COMMON 


2. 2. 2. 5 Prefilter and Multi look Integrator 

The general multilook (for four looks) prefilter situation is depicted 
for the input signal to the prefilter shown In Figure 15. The signal 
spectrum, which has been sampled at the sampling frequency f^, is offset 
from the center by an amount fQ^tf The objective is to filter the 
signal into the four individual looks and decimate the sampling 
frequency by a facto«" of four. The clutter lock circuit will detect 
the offset frequency and this will be applied to the overall spectrum 
(point (a)) to center it about zero (point (b)). Each of the four 
individual look offset frequencies f -| , fg. f3» and f^ are then applied 
and the result is low pass filtered to isolate the individual look 
bands. Since the bandwidth for each channel has been reduced by a 
factor of four, the sampling rate can also be reduced (decimated). 
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Figure 15 

Multilook Filtering Functions 


An implementation of the filtering and decimation function is given in 
Figure 16. A 32nd order low pass filter is used to achieve the de- 
sired passband and ripple characteristics. This is obtained from the ® 

approximation given in reference 2, 

In 2 - In ty 
AF« 

) « 

where AF is the passband to stopband transition width, M is order 
of the filter, and ^2 the passband ripple. 
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If we set AF * <05, and 1 2 obtain N^32. In the Implsnen- 

tatlon of Figure 16, the decimation Is accomplished by the serial to 
parallel shift register tap registers operating at a reduced (4 to 1) 
clock rate. The decimation Is exploited by operating the filter 
coefficient multipliers at the Input clock rate and accumulating four 
products prior to the formation of the adder tree. 
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CLOCK 


CLOCK/4 



FILTER COEFF 


CLOCK- 


I • * • 

-“MS 

^ ACC ACC ACC ACC 

T 


ACCUMULATE 4 SAMPLES 
IN SEQUENCE 


CLOCK/4 


Figure 16 

Low Pass Filter 
(32nd Order) Decimate by 4 


An implementation of the multilook Integrator Is shown in Figure 17a. 
The method shown with four Individual look integrators Is preferred 
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over the delay line approach since less raonory Is used. If X is the 
product of range cells times azimuth compression £.amples per look, 
the storage In Figure 17a Is 4X. A delay line approach (Indicated In 
inset) would require 3X + 2X + X » 6X. 



Multilook Integrator 
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2. 2. 2. 6 Time Domain Hardware Summary 

Table 6 lists a hardware summary of the time domain processing 
options. The totals assume 11 MHz clock rates In the correlator 
modules. AlthOLch this speed Is achievable. It Is somewhat high for 
TTL application. We feel It Is practical In this Instance, since it 
is an on-module clock rate. Most of the module Interface rates are 
at one fourth this rate or 2.75 MHz. Nevertheless, the total number 
of modules Is 3 to 4 times as high as the number expected for an FFT 
convolver Implementation. 
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Table 6 
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Quality 


Time Domain Module Summaries 


NUMBER REQUIRED 


FUNCTION 

TYPE A 

TYPES 

REMARKS 

RANGE CORRELATOR 

130 

130 

llMHz CLOCK. ~200 CIRCUITS EACH 

RANGE CONTROL 

1 

1 


PREFILTER CHANNEL 
MODULATOR AND CONTROL 

4 

4 

11 MHz IN. 2.75MHz OUT 

INTERPOLATOR 

m 

4 


AZIMUTH CORRELATOR 

420 

420 

ll.OMHz CLOCK. 4 CHANNELS / MODULE 

CONTROLLER 

1 

1 


INTERPOLATOR AND MAGNITUDE 

1 

1 

4 MODULES /INTEGRATOR 

MULTI LOOK ACCUMULATOR 

16 

16 

500K WORDS /MODULE 

CONTROLLER 

1 

1 


SYSTEM CONTROLLER 

1 

1 


INTERFACE /TEST 

1 

1_ 



577 

581 



2.2.3 FFT Convolver 

Azimuth processing, using an FFT convolver, has been divided Into 
two approaches: Case 1, discussed here, where the range and azimuth 

correlation are operated separately and Case 2, discussed In Section 
2.2.5, where a two dimensional convolution Is performed either Inclu- 
ding, or not Including, the range compression function. Range com- 
pression alone using FFT convolution Is discussed separately In 
Section 2.2.6. 

FFT azimuth correlation Is accomplished by forming an azimuth matched 
filter moving along constant range lines of the radar return. The 
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filter is formed by transforming the data to the frequency domain, 
multiplying by the complex reference of the azimuth focusing function 
«nd performing the inverse transform to convert the signal back to the 
time domain. The frequency transform filter thus formed performs 
a circular convolution of the azimuth focusing function with the 
signal data within the FFT window. 

With multiple look processing, separate frequency filters must be 
provided for each look. This is done as shown in Figure 17 b with 
separate spectral references for each look and separate inverse 
FFT's. Because the frequency band coverage of a single look is a 
fraction of the total, the frequency samples can be decimated (reduced 
sample rate) prior to inverse FFT processing. In the example of 
Figure 17b with four looks, the number of inverse FFT points is reduced 
by a factor of four to coincide with a four to one reduction in 
bandwidth. 


• SELECTED • RETAIN 
FREQUENCIES 2S6 POINTS 


FROM 

RANGE 

PROCESSOR 



LOOK 4 REF 


Figure 17 b 


FFT Convolution-Azimuth Processing 
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The ability to perform range migration compensation In the FFT 
convolver approach Is due to the unique relationship of range 
migration and frequency coefficients In the azimuth processor 
as Illustrated In Figure 18. As Indicated In Figure 8, the range 
migration correction Is a function of beam angle and the frequency 
coefficients are represontative of beam angle. Therefore, range 
migration correction by shifting frequency coefficients In range 
will correct for all of the targets with energy In that coefficient. 
We have shown (Section 2.3) that for adequate performance the range 
migration correction must be refined using Interpolated sample 
points. 


FREQUENCY 


RANGE 

BINS 



Figure 18 

Unique Relationship Between Range 
Migration and Frequency 


We have selected the FFT convolver algorithm for the recoiranended 
AOSP design and discussion of Its implanentatlon Is contained In 
Section 2.5. 


33 



2.2.4 Subarray Processing 

Subarray (or step transform) SAR processing Is performed by con- 
verting segmented received time data to the frequency domain and 
resolving the Individual doppler trajectory of the map elements. 
Partitioning the data Into small groups or subarrays permits range 
walk and other correction factors to be Incorporated. This par- 
titioning also enables a continuous strip map to be processed 
In an efficient manner by adding only new subarray data to the al- 
ready accumulated data and dropping the old subarray data. 

SAR processing using subarrays consists of the following five steps. 

1. Focus data over short subarray time. 

2. Store subarray data over large array length. 

3. Combine subarray outputs to form large array. 

4. Compute new subarray data entering large aperture. 

5. Drop old subarray leaving large aperture. 

Figure 19 Illustrates the step transform algorithm. Since the data 
Is being continuously fed Into the step transform processor, the 
overlapped subarray computation can be performed as a running process. 
Similarly, the output fine resolution can be computed as soon as a 
full aperture of subarrays have been computed. This Is also done as 
a running process where new subarray data Is added while the old data 
Is dropped. 

Factors to be considered In Implementing the subarray approach 
Include Input spectral aliasing, short aperture sidelobes, and 
straddling loss. These and other factors have been analyzed and are 
c?lscussed In Reference 3 . 
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Figure 19 

Block Diagram of Subarray Azimuth Processor 

The functional diagram of Figure 19 Illustrates the subarray approach 
for an Infinite deramping function. The azimuth processing begins 
with the application of the Infinite deramping function following 
range pulse compression. The reference function applied Is Illustrated 
In Figure 20. It Is essentially a continuous frequency ramp 
wrapping around at the sample rate. Total frequency coverage within 
the sample rate will cover the entire SAR antenna beam response 
Including multiple looks and guard band. A single target will 
encompass the full response. 

After demodulation, the data Is collected In a subarray memory which 
forms overlapping subarray segments which are amplitude weighted 
and processed In the subarray focusing FFT. Data Is then placed 
In a corner turning memory whose output Is first processed for range 
walk compensation prior to final subarray combining. The subarray 
combining data Is Integrated along paths as Indicated In Figure 21. 
Coherent Integration occurs across a single look at a single coarse 
azimuth frequency with multiple look, non-coherent Integration 
following along the same frequency. 
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Figure 20 

Continuous Deramping and Subarray Processing 
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Figure 21 

Deramp Coarse Resolution 
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Bulk storage In the subarray processor can be estimated by reference 
to Figure 22. The Integration must occur along horizontal paths for 
each look. As new subarray data Is received Integration Intervals 
as a function of frequency become filled. These are then read out 
and processed for fine re'wolutlon. Since the data can be discarded 
after Integration the minimum storage requirements for each range cell 
will be contained within the four triangular patches with width 
covering a look Integration time and height corresponding to the 
look bandwidth. Thus, If there are M subarray coefficients and N 
subarrays per look the minimum storage Is MN/2 or MN If a double 
buffered Implementation Is used. For the ERSAR case, the PRF Is 2500 
and the look Integration Is 0.7 seconds. The Input time-bandwidth 
product Is then 1750. Using a subarray size of 64 points and an 
overlap factor of 2 to 1 , the number of subarrays per look becomes 54. 
The number of spectral coefficients used Is 4/5th of 64 = 52. The 
minimum memory storage is therefore 52 x 54/2 = 1404 per range cell 
or about 5.7 x 10® words to encompass 4086 range cells. A programmable 
design of the bulk storage unit In a subarray processor Is difficult 
and It may simplify its control if a double-buffer arrangement can be 
used. In this case the total memory would be 11.4 x 10® words. 

The complex computation rate for the baseline can be estimated per look 
interval per range cell as follows: 


Continuous deramp function: .7 x 

2500 = 1750 

First FFT Wtg: 54 x 64 

3456 

First FFT: 54 x 32 x 6 

10,368 

Second FFT Wtg: 52 x 64 

3328 

Second FFT: 52 x 32x 6 

9984 

Total per range cell 

28,886 

Total per 4086 cells 

113 X 10 

Rate per second 

168 X 10 


The rate of 168 x 10® assumes equal computational load for the 
weighting function and the FFT butterfly computations. If the 
weighting is accomplished by combining it with the FFT process, half 
of the weighting operation are eliminated and the net computation rate 
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becomes 148 x 10 complex operations per second. 







MINIMUM STORAGE - MN/2 
DOUBLE BUFFER STORAGE • MN 


Figure 22 

Bulk Storage in Subarray Process 
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2.2.5 Two-dimensional Convolution 

The 2-dimensional convolution approach to SAR azimuth processing 
has been divided into two parts, one in which the range correlation 
is included and the second for azimuth compression only and the 
2-dimensional process is applied to correct for range migraticn. 

The latter approach also includes a hybrid FFT time domain technique. 

2. 2. 5.1 Full Range Azimuth Correlation 2-0 Process 
Two dimensional convolution as illustrated in Figure 23 is an ideal 
approach in terms of computational efficicency for a single reference, 
but the small depth of focus for a single reference spectrum results in 
reduced efficiency. For example, the azimuth reference function 
must be updated every 6 km of swath to maintain single look focus. 

Over a 40 km swath, the reference v/ould have to be adjusted at least 
9 times if a 1024 x 4096 2-0 array were used. The process also rules 
out the incorporation of multilook registration in the process since 
this requires a reference update every 8 range cells. 


DATA SPECTRUM 


RANGE 



TWO-DIMENSIONAL 
FFT PROCESS 








REFERENCE SPECTRUM 


T DEPTH OF 

i FOCUS 


Figure 23 

Two-Dimensional Convolution Using FFT Processing 
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2. 2. 5. 2 2-D Convolution With Ranga Compresstd Data 

A 2-D convolution azimuth correletion procedure can also be anployed 
on the range compressed data to accomodate range walk. This fu.ction 
can be done using the standard FFT algorithms (4) or other more recent 
algorithms such as the fast polynomial transform (5). The fast 
polynomial transform has been applied to SAR processing (6,7) and a 
short summary will be provided here. 

Consider the convolution of the range compressed data ^a(m,n)J with 
reference £r (n,m)^ , each of size MxN. For the SAR processing in 
mind, M is the number of range cells for which a single focusing 
function could cover, typically 128 or 256, and N is the azimuth 
coverage, typically 4096 or 8192. 

From the data array ^a(m,n , we generate (r+1) subarrays, where 
r ■ 1 + log 2 (N/M). The i-th subarray ^a^ (m,n)^ is of size Mx(N/2’), 
i = 1,2,...,! + 1092 (N/M) except array ^a°(m,n)^ which is of size 
Mx(M/2). To generate ^a^ (m,n)^ , we take each row of ^a(m,n)3, 
partition in the middle, take the sum and difference of the corres- 
ponding elements: 


a^(m,n) = a(m,n) - a(m,n+M/2) 

-1 


0 - n<N/2 - 1 

a'(m,n) = a(m,n) + a(m,n+M/2) 

Of the two arrays generated ^a^(m,n)*^ and ^a^(m,n)^. The first one 
is kept and the second is used to generate a^(m,n), row by row, as 
fol lows : 

a^(m,n) * a^ (m,i.)-a^ (m,n + N/4) 

.2 -1 -1 0 ^ n ^ J - 1 

a (m,n) = a (m,n) + a (m,n + N/4) 

This process is illustrated in Figure 23. 


* The notation is somewhat different from references (6) and (7). 
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The process is repeated r times, the last two arrays are£a'^(m,n) 
and^a°(m,n)^ as indicated in Figure 24. The operation takes a 
total of 2MN additions. 

The operation is also precomputed on the reference array ^r(m,n)^ 
and stored. 


N 


■ ■ 



iX: 

I ■ I 

a*2(n,n) «2(ni,n) 


Bith row of a'M»,n) 
Genoratton of a^la.n) 


(n^*i row of a (m.n) 
Ggnoratlon of aMm.n) 


Figure 24a. Generation of ^a^(m,n)^ from ^a(m,n)^ 

of ^a^(m,n)^ from (m,n)^ . 
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Figure 24b. Repetition of Figure 24a Process 


FIGURE 24. Polynomial Transform Operations 


41 


The next step is to perform the polynomial transform of the arrays 
{a^(m,n >] and Jb^(m.n)^ , i ■ 0, l,...,r. The structure of the 
polynomial transform is very similar to radix 2 FFT, except each row 
in the array is treated as a single element. To illustrate with a 

"DiF" structure applied to the (m,n)"^ array, we first take two 
rows M/2 apart, and take the sum and difference of the corresponding 
elements, i.e.. 


a^(m,n) -a^(m+M/2,n) 0 < N/2-1 

The sum is left alone, but the difference is cyclic shifted (or wrapped 
around) by a prescribed amount (depending on the corresponding twiddle 
factor in the FFT structure) and the sign of the vrapped around 
part is changed as shown in Figure 25. These two rows are put bark 
in the array and two more rows are taken out and operated on in this 
manner. After M/2 such operation'^, the first stage is completed. 
Another set of operation then commences, as in the 2nd stage of the 
FFT. That is, rows that are M/4 apart are operated on pairwise. 

This operation is repeated for the log 2 M stages. 




ELEMENT BY ELEMENT 
ADDITION AND SUBTRACTION 


Figure 25. First Stage of Polynomial Transform 
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The "butterfly" operation in the polynomial transform can be pipe- 
lined by inserting the results back into the array in a criss- 
cross manner (8), as follows. For the first stage, the first pair 
Ut.’co.n)] and ^a^(M/2,nl^,the results are, for example fa^(0,n)T . 

The 2nd pair to operated on by the butterfly should be ^a (M/4,n)^ , 
and ^a^(3M/4,n)^ , obtaining the results ^a|.(M/4,n)^ and ^a]_(M/4,n)^. 

Now the result fa^(0,n)^ should be stored in row 0,TaNo,n)J in row 
M/4, ^a^(M/4,n)J in row M/2, and ^a]_(M/4,n)^ in row 3M/4. 

The total number of additions for transforming the (m,n)3 array 
is Nlog 2 M. So the total number of additions for transforming all 
the ^a^m,n)^ arrays is 2 Nlog 2 M. 

Again, we note that the polynomial transforms of the reference arrays 
|r^m,n)3 can be precomputed and stored. 

Me denote the polj^omial transforms of ^^(m,n)J and ^^(m,n)3 by 
^(m,n)^ and £r^(m,n^ respectively. 

The next step is a cyclic like convolution of the rows of^^(m,n)^ 
with . The difference from the ordinary cyclic convolution 

is that the wrapped around part has a sign change. References (6) and (7) 
suggests the use of an idea proposed by Arambepola and Rayner (9) which 
accomplishes the convolution via a transform algorithm that is only 
a slight modification of the FFT algorithm by changing the twiddle 
factors. One could, on the other hand, use the straightforward FFT 
tc compute the mncyclic convolution and simply do a circular shift 
and change the signs of part of the result. This step requires 
2MN(log2N-3) additions and N(^ + log2N-3) multiplications. 

Finally these (r + 2) arrays are combined to give the cyclic con- 
volution Oi the data ^a(m,n3 af’fl the reference £r(m,n)3. This 
step involves the merging of arrays of sizes Mx(N/2^). The merge 
starts with the summing and differencing, element by element, two 
smallest arrays, each of size Mx(N/2'^) (Figure 26). The result is a 
MxM array. This is then merged with the next size array (i=r-l), 
using the same element by element sum and difference. This process 
is continued until a final array of MxN is obtained. This is the 
result of the cyclic convolution. 
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Figure 26 

Merging of two MxM/2 arrays 

The advantages claimed by the fast Polynomial Transform method are; 

1. It requires fewer multiplications than using the FFT. 

2. The process of Polynomial transformation can be performed in 
parallel . 

The comparison with FFT in terms of arithmetic operations are sum- 
marized in Table 1 of reference 6. The number of additions are com- 
parable and the number of multiplications have a ratio of approximately 
5:3 in favor of the polynomial transform. Bergland (10) has shown in 
1968 that there is a 30% saving in multiplications if one uses radix 
4 rather than radix 2. Recently Nakayama (11) proposed a mixed 
decimation FFT algorithm which results in a 15% saving over the con- 
ventional FFT for the radix 2 case. In his algorithm, the two input 
butterfly computation ke» .p’ is retained. The saving in the radix 4 
case is about 5 ^ 8 %. 

The computation of 2-dimensional transforms is conventionally carried 
out by first transforming each row and then transforming each column. 

It has been shown (12) that a simultaneous row-column decimation 
would result in a 25% saving. It is therefore possible to combine the 
simultaneous 2-dimensional decimation with the other schemes mentioned 
previously to obtain a substantial saving over the straightforward 
FFT approach. Indeed, it has been shown (13) that a 40 a^ 0% saving 
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possible. Needless to say. the saving Is achieved at the expense of 
Increased complexity of the algorithm. But the polynomial transform 
approach, requiring operations on row vectors of different length 
also Increases the complexity considerably over that of a conventio- 
nal FFT. Thus, the first purported advantage of savings In arithmetic 
operations fa achievable with other approaches with perhaps a simpler 
structure of the hardware processor. This latter point definitely 
would require further study to establish. 

The second advantage of the polynomial transform Is also not con- 

i 

vincingly demonstrated. The arrays a (m,n) are of different sizes, 
the row sizes of successive arrays are two to one. Although these 
operations can be carried out In parallel. It Is doubtful that the 
saving would be substantial. A considerable portion of the hardware 
would be Idle. Essentially the speed advantage of using parallel 
operation would be about 2 to 1 but the hardware would be Idling 
about half of the time. 

The tentative conclusion we hive reached at this time Is that, In terms 
of simplicity of structure, a processor based on the conventional 
FFT or some minor modifications appears to be most attractive approach 
for a two-dimensional convolution. 

2. 2. 5. 3 Hybrid 2-D Process 

A final 2-1) process ha;, been suggested by Wu and Liu (14). This 
approach. Indicated In Figure 27, employs FFT convolution for azimuth 
correlation and a tapped delay line convolver for the range migration 
correction. Ordering of data appropriately for the tapped delay line 
convolver Is accomplished by doing the azimuth processing with a 
multiplexed FFT structure. This method, which requires a larger 
memory storage, outputs data from the first FFT on a single frequency 
coefficient sequential range basis. Range migration correction can 
then be applied with a tapped delay line of length equal to the range 
migration. Modularity and an incremental growth and implementation are 
also difficult to achieve. 

In principal this method is functionally equivalent to the FFT convolver 
algorithm using range interpolation in the range migration compensation. The 
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memory required is at about 2.5 times the minimum memory storage required 
in the FFT convolver corner turning memory. Since the minimum memory 
storage is about 17 x 10^ words, the total impact on cost is large. 
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2.2.6. Range Correlation 

2. 2. 6.1 Range FFT Convolver 

The forward- inverse FFT frequency domain matched filter, Figure 28 is 
now an established technique for digital convolution for many appli- 
cations (15). Matched filtering is accomplished by spectral domain 
multiplication with the inverse FFT providing the matched filtered 
time domain output. 

The advantages of the forward- inverse FFT matched filter are well 
known and include: 

§ A computation-efficient algorithm. 

• Complete waveform flexibility. 

• Adaptable to modular construction. 
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• Precise, predictable perfortnan'*e level. 

• Implementation can be adjusted to accommodate different processing 
speed requirements. 

• Will benefit from advanced tcrchnology developments. 

• Direct interpolation of output data possible. 



> iiatiuwMtuiocami^ieN 


uicMoniTiitniacaoii ■( 


Figure 23 

FFT Matched Filter System 

The input data can he processed from real sampled data as might be 
acquired with the intention of processing in an optical processor. 

Using the FFT for providing the real-complex conversion is the most 
convenient approach. If this is done the input data window is 
normally twice as long, but two real channels can be processed simul- 
taneously in tne complex FFT window. 

The FFT matched filter (range correlator) will functionally have the 
form shown in Figure 29. Because the range swath interval is a fraction 
of the total radar pulse repetition interval (PRI) the range data is 
normally buffere': after A/D conversion to slow down the data rate to 
the range correlator. 
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Figure 29 

FFT Convolution-Range Processing 

A tradeoff exists In setting the FFT processor parameters of an FFT 
convolver. In general, the larger the size of the FFT, the more 
efficient the computations. For the baseline system range processor, 
the total number range samples to be processed per pulse Is 4000 
(range Image swath) plus 475 (uncompressed pulse length) plus 86 
(range migration) for a total 4561 samples. This can be processed 
as a sliding aperture convolver with an FFT size of 1024 points or 
as a batch processor of 8192 or 4096 points. If the maximum FFT 
size Is 4096, a separate FFT of 1024 must be processed to cover the 
total range Interval . The net computation rate for the three cases 
is given In Table 8 which shows that the most computationally 
efficient size for the baseline system Is the 4096 point FFT. 



Table 8 


Range FFT Convolver Alternatives 


CHARACTERISTIC 

SLIDING APERTURE 

MAX'SWATH PROCESSOR 

f FT SIZE 


— 

COMPUTATIONS 
PER PULSE 

8|si2(2la92l024«2) 

(W.112) 

8192 FFT 

4096(2 l 09 j> 8192*2) 
(114,6881 

4096* 1024 FFT 

2048(210924096 * 21 * 512(221 
(64.5121 

COMPUTATION RATE 
(COMPLEX . . 

OPERATIONS /seel 

225 X # 

8192 FR : 287x10* 
4096* 1024: 161x10* 


2, 2. 6. 2 Step Transform Unsar Frequency Modulated (LFM) Signal 

Matched Filter 

The step transform subarray algorithm was originally conceived as 
a technique for LFM matched filter processing (16,17). The algorithm 
Is shown symbolically In Figure 30. A received LFM ramp Is demodulated 
by a sawtooth ramp with the same slope, producing CW frequency seg- 
ments "stepped" by the frequency span of the sawtooth (^f) and of 
length equal to a tooth (<0T). A spectral analysis of the successive 
CU segment produces a set of frequency response functions that are 
stored In a time-frequency data reorder memory shown In Figure 32. 

Data as read from the time frequency matrix along diagonals results 
In a linear phase shift along the diagonal producing the "fine" range 
resolution at the processor output. Weighting Is applied prior to 
the second FFT to reduce range sidelobes. 

The step transform pulse compression technique provides a natural 
segmentation of data as a function of range. It offers a means of 
applying range migration compensation to small blocks of range 
samples. Its main disadvantage Is that It Is not universally pro- 
grammable for waveform type or time-bandwidth product. 
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Figure 30 

Step Transform LFM Pulse Compression Algorithm 



Figure 31 

Step Transform LFM Range Pulse 
Compression Processing 

(For maximum efficiency, DFTs are implemented 

as FFTs) 
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2. 2. 6. 3 Digital Tapped Delay Line Correlator 

Tapped Delay Line Pulse Compression Filters - The output y(t) of a 
matched filter Implemented via convolution of an Input signal s(t) 
with the Impulse response of the matched filter h(t) « x(-t), where 
x(t) Is the transmitted waveform, Is: 

y(t) -J** S(r)H* (t-r)dr 

00 

The signal and filter functions are generally represented as complex 
Inphase and quadrature samples for sample data operations. Thus, a 
physical realization of the filter function must accomodate the 
complex multiplication operation. 

(a + jb) (c - jd) = ac + bd + j(bc - ad) 

A tapped delay line matched filter then takes the form of Figure 32. 
The total number multipliers Is over 4(TW)‘ where TW Is the time band 
width product of the waveform. With a maximum TW product of 660 for 
the ADSP, the total number of taps required Is 765. 

A time domain range processor Is very simple to use. !t Is only 
necessary to Input the radar waveform of arbitrary length into the 
reference register together with the appropriate number of zeros. 

If a single unit Is used as In Figure 32, however, It must operate at 
a clock rate of 11 MHz. The system would be partitioned Into 130 
modules with up to 200 circuits per module all operating synchro- 
nously at an 11 MHz clock rate. The total number of circuits for 
the range correlator would In this case be about equal to the total 
number required by both the range and azimuth correlators using FFT 
convolver algorithm. 
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Figure 32 

'ime Domain Range Convolver Controls 
2.2.7 Algorithm Study Summary 

A graphical depiction of the SAR baseline systan performance para- 
meters for the azimuth processor are given in Figure 33. It shows 
the relationship of the total beam integration time of 28 seconds 
to the radar PRF. The time-bandwidth product per look is 0.7 
seconds times 500Hz = 350. Using these parameters, a compirison was 
made of the azimuth processor memory size and computation rate for 
the three main processing algorithms as a function of the number 
of lookSr 
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Figure 33 

Baseline Azimuth Parameters 


The results are given In Table 9. It shows results as generally 
expected with the subarray approach having the lowest computation 
rate and the time domain processor exhibiting the smallest memory. 
However, the computation rate of the time domain processor Is much 
higher than the other approaches. This high computation rate 
translates Into hardware as danonstrated In Section 2.2.2 and 
eliminates the time domain approach as a contender. The tradeoff 
between complexity of control and hardware Is a factor between the 
FFT convolver and subarray approaches. 

The 2>0 convolution approaches have not been Included since they are 
similar to the FFT convolver algorithm. In the approaches using 2-D 
convolution of range compressed data the question Is of the relative 
advantage of 2-D convolution for range migration compensation versus 
range interpolation. Simulation results presented In Section 2.3 
indicate that a simple two point interpolator will give an adequate 
integrated sidelobe performance level. A four point interpolator 
is close to the ideal and its implenentation in the selected design 
is not a computational burden. 

Based upon its greater ease and completeness of programmabll ity and 
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Its lower cost for development, the FFT convolver algorithm has 
been selected for the AOSP recommended design. Further details 
of the design and selection Issues are given In Section 3.S, The 
design and architecture have also lent themselves to the selection 
of the FFT convolver for the range processor. 


Table 9 

Look- Parameter Azimuth Processor Variations 
for Baseline SAR System 


LOOKS 

CONVifi 

T 

.VER 

SUBARRAY 

TIME DOMAIN | 

COMP 

RATE 

MEMORY 


MEMORY 


MEMORY 

1 

306 

67.1 

171 

45.6 

QQII 

28.0 

2 

276 

33.6 

160 

22.8 

19,250 

14.0 

4 

245 

16.8 

148 

11.4 

4,813 

7.0 

8 

215 

8.4 

138 

5.7 

1,203 

3.5 


• COMPUTATION RATE IN MILLIONS OF COMPLEX 
OPERATIONS PER SECOND 

f MEMORY STORAGE IN MILLIONS OF WORDS 

• TOTALS EXCLUDE RANGE MIGRATION COMPENSATION 
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2.3 Performance Levels 

The selection of a signal processing configuration for the ADSP 
required the determination and comparison of the relative performance 
of each candidate algorithm. The analysis/simulation of the algorithms 
Is described In this section. The simulation results permitted cost 
versus performance tradeoffs. A design requirement that was particu- 
larly Important In the tradeoffs was the Integrated sidelobe level 
(ISL). A design goal of -20 db ISL In two dimensions (range and 
azimuth) was established as a minimum requirement. 

2.3.1 Performance Procedure 

The organization of the simulation was broken Into two parts. The 
first part Is the generation of a representative time signal of 
Interest. The second was the processing of the signal using two 
different algorithms, the FFT-convolver and the step transform sub- 
array. 

2. 3. 1.1 Generation of the time Signal 

The objective of the first part of the simulation was to generate the 
time response of a point target. The simulated target was assumed 
to have already been matched filtered In the range dimension by FFT 
convolution. The range response was Hamming weighted and was evalua- 
ted exactly at the sampling points. If the target peaked at a non- 
sampled point, the response at the sampling points was generated re 
lative to the peak. The phase of the target was generated using 

(KO= 

where V Is the velocity of the spacecraft (assumed to be a constant), 
t Is the relative sampling tlme,^ Is the wavelength of the transmitted 
pulse, and R Is the slant range to the point source. 

The effect of the range migration was obtained by displacing the 
Hamming weighted function. Range migration is composed of linear 
and quadratic components. The linear range migration component Is 
primarily due to the earth's rotation. The quadratic component varies 
In a non-linear manner as a function of beam angle. Figure 34 shows 
the two components of the range migration as they were modeled in the 
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simulation. 



Figure 34 

Quadratic and Linear Range Migration 

Linear range migration is simply a slope parameter in the simulation. 
The quadratic displacement is given by 

(v^tyAtVCc/As) 

where V is the velocity of the spacecraft, t is the relative time, R 
is the slant range, c is the velocity of light, and B is the band- 
width of the system. Figure 35 shows the relationship between the 
quadratic range migration (35a) and the corresponding linear FM 
function (35b). The maximum range, Rm, in Figure 35a corresponds to 
the maximum and minimum frequencies and -i^^/2) in Figure 35b. 

The minimum range migration corresponds to the zero crossing point of 
Figure 35b. The signal is present along limited portions of the 
spectrum as a function of time and generates a linear FM (LFM) 
response. The generated signal migrates to different range cells as a 
function of time and beam angle. 
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Figure 35 

Generation of the Time Signal 


2.3.2 FFT - Convolver Algorithm 
2.3.2. 1 Range Migration Compensation 

The FFT-Convolver SAR processing algorithm was simulated first. In 
this approach, the azimuth data is transformed to the frequency 
domain along constant range lines. In Figure 35a, the transform is 
performed along the horizontal lines. Since the energy of a point 
source is dispersed into different range lines, compensation must be 
provided prior to matched filtering in the azimuth direction. Figure 
36 shows the Fourier transform of N range cells. The compensation 
for the range migration consists of moving the Fourier coefficients 
into the proper range cells. The number of ran._ migration cells is 

equal to 

where Ais the wavelength of the transmitter, R is the minimum slant 
range, V is the velocity of the spacecraft, R^ is the range cell 
width andxis the Fourier coefficient number. 
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Figure 36 
Range Migration 


The amount of range walk as a function of frequency allows for the 
proper alignment of whe coefficients. 

The simpliest adjustment to the coefficient would be to simply 
slide the data by integral steps (8). For example if the migration 
for a given coefficient were 5.4 range cells then the data would be 
obtained from the coefficient 5 range cells offset. For a migration 
of 6.7 range cells, the data would be obtained from 7 range cells array. 
This itfithod, while simple, would tend to cause some output distortion. 

Interpolation of the data shift eliminates the effects of the dis- 
continuities of integral transfers. Interpolators of two, three and 
four points were simulated. Table 10 gives the equations of the 
interpolators used in the simulations. 
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Table 10 

Interpolation Formulas 

• 2pt: u-pl Xq ^ pX) 

• 3pti xo-i • ii-iA xo * xj 

. 4pt.. X., . x„ . 

X, . ISI^ X2 

2.3. 2.2 Azimuth Matched Filter 

Following the range walk compensation, each of the range lines are 
matched filtered by multiplying the transform of the azimuth focusing 
function. The azimuth focus function has a linear FM slope of 

a » 2'lh^/Ril 

and Ihi gth determined by the maximum Integration line. The simulation 
has the option of applying a non-unlfcrm walghting across the time 
response of the azimuth matched filter. The location of the time 
aperture of the azimuth matched filter depends upon which of the four 
looks Is being processed. 

2. 3. 2. 3 Weighting 

Figure 37 outlines the terms that are used In the discussion of the 
simulation and results. The waveform In Figure 37 represents the 
compressed response In both range and azimuth. The malnlobe of the 
response can be defined In a number of ways. The most common de- 
finitions use widths about the peak to the 3 db points, 6 db points 
or to the first nulls. Sidelobes are then defined as everything that 
does not lie within the malnlobe. The peak sidelobe level versus the 
peak malnlobe (PSL/PML) Is defined as the largest value found outside 
the malnlobe divided by the largest value found within the malnlobe, 
usually expressed In db. The Integrated sidelobes versus the Integrated 
malnlobes (ISL/IMl) Is defined as the Integration of the region not 
found In the malnlobe divided by the Integration of the malnlobe region. 
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again usually expressed in 1 dB. The PSL/PML and the ISL/IML can be given 
for one or two dimensions. The one dimension (ID) case uses a single range 
line containing the mainlobe peak. The two dimension case (2D) is the two 
dimensional plane of azimuth and range, -sing the definitions of Figure 37, 
Table 11 gives a brief review of differe. v. •'ighting functions. Weighting 
selection represents a compromise between resolution or width of the peak 
response and sidelobe levels. This is clearly seen in comparing the increase 
in the mainlobe width (3 dB and null columns of Table 11) to the decrease 
in peak sidelobe and integrated sidelobe levels (the columns labeled PSL/PML 
and ISL/IML). The results of Table 11 were obtained using a signal aperture 
of 64 samples and a total time history of 1024 samples (the rest of the 



PSL^L -- PEAK SIDELOBE VERSOS THE MAINLOBE (dBI 

ISLHML - INTEGRATED SIDELOBES VERSUS THE INTEGRATED 
MAINLOBE (dBI 

• ID- ALONG AZIMUTH DIRECTION 

• 2D -- ALONG AZIMUTH AND RANGE DIRECTION 

• MAINLOBE - POINTS LYING WITHIN THE FIRST NULLS 
FROM THE PEAK 

• SIDELOBI'i - ALL POINTS NOT IN THE MAINLOBE 

• 3dB - WIDTH OF THE MAINLOBE AT THE 3dB POINTS 
t 6dB - WIDTH OF THE MAINLOBE AT THE 6dB POINTS 

« NULL - WIDTH OF THE MAINLOBE TO THE FIRST NULLS 


Figure 37 

Sidelobe Definitions 
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Table 11 

Ideal Characteristics of Weighting Functions 


LENGTH: 64 

FFT: 1,024 

TYPE 

PSL/ML 

RECTANGULAR 

13.26 

HAMMING 

42.48 

BARTLETT 

2&53 

HANNING 

31.50 

BUCKMAN 

5a 14 

25 dB TAYLOR 

24.34 

30 dB TAYLOR 

30.98 

35 dB TAYLOR 

35.93 


ISL/IML 

3dB 

NULL 

9.68 

1.0 

2.2 

3442 

1.5 

4.5 

25.31) 

1.5 

4.3 

32.91 

1.7 

45 

5489 

1.9 

6.6 

1496 

1.2 

3.0 

2455 

1.3 

3.4 

27.69 

1.3 

3.7 


•PSL/ML - PEAK St DELOBE TO MAINLOBE RATIO (dBt 
*ISL/IML - INTEGRATED SIDELOBES TO INTEGRATED MAINLOBE (dBI 
*3 dB - REUTIVE WIDTH OF MAINLOBE TO 3 dB POINTS 
*NULL -- RELATIVE WIDTH OF MAINLOBE TO FIRST NULL POINTS 


2. 3. 2. 4 Time Domain Output 

The final processing of the signal In the azimuth FFT convolver Is to 
apply the Inverse Fourier transform to obtain the time domain response. 
In the simulation, the time domain output Is not decimated as would 
be the case In an actual signal processor. This reduces the measure- 
ment error since the nondecimated output provides a higher output 
sampling rate. 

2. 3. 2. 5 Interpolation and Weighting Simulation Results 

A number of simulation runs were performed to measure the peak sidelobe 
levels and the Integrated sidelobe levels for different order 
Interpolators. The system parameters used In the simulations were a 
constant spacecraft velocity of 7.45 km/sec, a slant range of 850 km, 

22 MHz sampling rate In the range direction, a pulse repetition rate 
of 1.6 kHz and an azimuth coverage of 4096 pulses. The signal was 
placed In the middle of the scan and the scan was processed for four 
looks. The radar system pulse bandwidth In the range dimension assumed 
to be 19 MHz (over-sampling In the range direction by 1.16). The 
transmitter frequency was assumed to be 1250 MHz. 
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Table 12 

FFT Convolver Interpolator Results 

• CASE 1; NO INTERPOLATION 

10 20 

LOOK I SI/ 1 ML ISL/IML PSL/ML 

1 26.3 12.4 19.8 

2 27.8 12.6 19.6 

• CASE 2: FOUR POINT INTERPOLATION 

10 20 

mim mm esim 

1 31.1 28.9 35.8 

2 3aS 32.6 36.3 

9 HEASUREHENTS IN DB 

I HAftllNG WEIGHTING 

Table 12 summarizes the performance obtained by inco 'operating a four 
point interpolator with range migration compensation. Since the signal 
was placed in the middle of the scan, look 3 has the same results as 
look 2, and look 4 the same as look 1. Looks 3 and 4 are thus not 
shown in Table 12. Table 12 shows that the ISL/IML performance is 
improved by the four point interpolator over no interpolator in one 
dimension only (column labeled 'ID') by only 5 db at a 30 db level. 

The improvement is of the order of 20 db however when a two dimensional 
integration is examined (column labeled '20'). The peak sidelobe 
level versus the mainlobe occurs in the one dimensional line and is 
thus the same for either ID or 2D (column labeled 'PSL/PML'). The 
large degradation in the ISL/IML with no interpolator is expected. 
Shifting of the Fourier coefficients without interpolation results 
in discontinuities that raise the level of the sidelobes in the range 
cells surrounding the signal. 

The results in Table 12 were obtained using a Hamming weighting function 
across the azimuth matched filter. Table 13 illustrates the reason 
for using a weighting function to reduce the sidelobe levels at the 
expense of a reduction in the resolution in a simulation. If no 
weighting were employed, the results would be ambigious on the need 
for an interpolator in the processor. The column labeled 'Reef shows 
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only a 1.5 to 2.0 improvement for the Interpolated case. What has 
happened is that the high sidelobes of the rectangular weighting 
function (window) have contributed to the Integrated sidelobe. Using 
a window with greater sidelobe suppression provides results which are 
sensitive to processing differences. A triangular window ( column 
labeled 'Triangle') suppresses the PSL/PML level of the window to 
27 db from 13 db for a rectangular window and also the falloff of the 
sidelobes is at a 12 db/octave rate versus 6 db/octave for rectangular. 
The mainlobe increases by approximately 40%; however, the difference 
due to processing can now be seen. The ISL/IML difference between 
no interpolator and a four point interpolator is approximately 12 db. 
Using Hatmning window which has a PSL/ML of 43 db, a falloff rate of 
6 db/octave and an ISL/IML of 34.42 db, the processing difference is 
seen to increase to approximately 20 db (column labeled 'Hamming'). 
Considering that the theo»*etica1 ISL/ML for a rectangular window is 
9.68 db and a PSL/ML of 13.26 db, the effects being seen in the "rect" 
column for the no-interpolation case are largely due to processing 
(the ISL/IML is approximately 8 db) , while for the four point 
interpolator it is the window (the ISL/IML is 9.7 db). On the other 
hand, both processor's performance is being measured when the Harming 
window is used. 


Table 13 

Effect of Weighting Functions, FFT Convolver 


I CASE 1; NO INTERPOLATION 


LOOK 

ISL/IML (PSL/ML) 1 

2 

AZIMUTH RESOLUTION 
RANGE RESOLUTION 


RECT TRIANGLE HAMMING 

7.8(13.5) 12.4(19.8) 12.4(19.8) 

&1 (12.8) - 12.7 (19.6) 

1 1.4 

CONSTANT 


• CASE 2: FOUR POINT INTERPOLATION 


ISL/IML (PSL/ML) 

AZIMUTH RESOLUTION 
RANGE RESOLUTION 


LOOK RECT 


1 9.8 (13.5) 

2 9.7 (13.2) 

1 


TRIANGLE HAMMING 

24.0 (26.7) 29.1 (35.8) 

25.1 (27.0) 32.8(36.3) 

1.4 1.4 


CONSTANT 
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2.3.3 Step Transform - Subarray Approach 
2. 3. 3.1 Description of Subarray Simulation 

The second algorithm simulated was the step transform-subarray. The 
first step in the process after range compression is the derrmping 
of the complex time signal by the azimuth reference LFM. The deramped 
signal is then sectioned into subarrays. The number of samples in 
each subarray and the amount of overlap between the subarrays were 
systems parameters examined in the simulation. An example of a sub- 
array length of 64 with 2 to 1 overlap is depicted in Figure 38. 

■ » OERAMPEO TIME DATA 

E 3 

— SI— 

C— S2— ] 

[—S3—] 

C— S4— ] 


I EXAMPLE 

t LENGTH OF SUBARRAY - 64 

• OVERLAP U/2) -- 32 

• SUBARRAYS 

• SI a-641 

• S2 (33-961 

• S3 (65-128) 

• S4 (97-160) 


Figure 38 
Subarray Formation 

After forming the subarrays, each subarray is transformed into the 
frequency domain using an FFT. The suoarrays are then phase corrected 
using a correction factor of eltp (-j2^kd/N) where k is the FFT 
coefficient number, d is the displacement between subarrays and N is 
the length of the FFT. The original signal has been transformed 
into a three dimensional array with the Fourier coefficient? on one 
axis, the subarray number on the second axis, and the range cell 
number on the third. Figure 39 shows a two dimensional slice of the 
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three dimensional array with the range dimension fixed and the FFT 
length 64. 

I COMPUTE THE FIRST FFT (EXAMPLE 64 LENGTH! 


1st COEFF - DC TERM 
2nd COEFF --A 
3rd COEFF - 2A 

33rd COEFF - T/2 

A • fs/LENGTHOFFFT 

fs • IFF • SAMPLING FREQUENCY 


Figure 39 

First FFT in Subarray Approach 

Assuming that the signal depicted in Figure 39 has no range migration, 
then the deramped subarrayed signal would appear entirely in the 
plane. If for example, the difference between the deramping reference 
and the signal was exactly 3 A and the signal began at the 256 sample, 
then the fifth subarray would be the first to see the signal and its 
transform would have the fourth Fourier coefficient as the highest of the 
64 coefficients (assuming no other signals are present). The fourth 
coefficient would be high for each subarray in which the signal was 
present. In the case where the signal exactly matches the reference, 
the first (DC coefficient) would be the highest of the 64. 

However, the signal does not stay within a single range cell. The range 
migration traverses a quadratic function in the third dimension. Consider 
three signals, one that exactly matches the reference, one that is 
delayed in time by an amount which corresponds to the signal peaking in 
the 2nd Fourie** coefficient and the third is delayed by an amount 
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equal to the third coefficient. These three signals are depicted in 
Figure 40a as 0, and 2 AT delay from zero time. Assuming that the 
first signal begins at the first sample, then the first coefficient 
of the first array would have to be range corrected the maximum 
amount. The first coefficient of the 2nd subarray would be corrected 
according to the range migration correction formula used in the FFT- 
Convolver approach. The amount of range migration and correction would 
diminish as the higher number subarrays are addressed until the 
subarray corresponding to the zero time of Figure 40a is addressed. 

This point is shown as the first point on the zero correction line 
of Figure 40b. The correction would then increase to the maximum. 

S'f^ze the signal was assumed to be matched to the reference, then all 
of the correction is applied to the first Fourier coefficient. 


Figure 40a 

Migration of Target' Across Subarray 




CUNIUIIK (U 

MAXIMUM 

CUIUItCTIUN 


/?inu 

COMIUCIION 


MAXIMUM 

COHRICTIUN 


SAMt DlIHlientON AFIM It 0 ACROSS SUHARHAV tUMI.NTS 

Figure 40 

Subarray Range Migration Correction 


The second signal in Figure 40a would first appear at a subarray 
number corresponding to the time delay AT. Once the signal starts, 
the correction function would proceed from the maximum amount down 
to zero and back up to the maximum in the same quadratic shape as the 
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first signal. The only difference Is that th« 'orrectlon Is per- 
formed on the 2nd FFT coefficient. The third signal would have the 
same correction performed on the 3rd coefficient. Figure 40b depicts 
the correction across the subarrays as the straight lines labeled "contour 
of maximum correction", "zero correction", and "maximum correction". 

The slope of the parallel lines would be given by 

slope » {-^s^ T)/(Ni^LFMcft) 

where fs Is the azimuth sampling frequency, N Is the subarray length, 

T Is the total time that the signal Is present, A LFM Is the total 
change In the LFM reference and d Is the number of points that the sub- 
arrays overlap. The correction scheme used for the range walk Is no 
different than that used In the FFT-Convolver algorithm. The data 
Is shifted in integral amounts for no Interpolation and Is Interpolated 
between points otherwise. 

The final part of the processing in the subarray approach Is to per- 
form the focusing FFT. The same slope that Is used In the range 
migration correction is used to address the subarrays for the focusing 
FFT. In the case of a non-integral slope, the nearest subarray Is 
chosen. The signal Is then padded with zeros which minimizes the 
measurement errors of the analysis. The FFT Is computed and the proper 
coefficients are selected from the FFT. Each set of coefficient are 
abutted to form the output. The selection of the proper coefficients 
is analogous to a comb filter bank where each output filter is tuned 
to only certain center frequencies and a given bandwidth about them. 

Figure 41 depicts the process of selecting the output values. 
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tNVAllO 

POINTS 


I 2 3 

SECOND FFT | i | 

FFT I fft 3 

% * * 

FF12 FFT4 


FIRST FFT COEFFICIENTS 


VAIID 

POINTS 


mum 

-32— 




umii 




OFTTRMINATION OF VAIIO COEFFICIENTS OF FOCUSING FFT (-| 
t «<f/vlFRFQUFNCY SPACING OF FIRST FFI( 

( d • »anoVI RLAPPING DISPIACEMENTI 

• AF'a (»/«yN‘lFREOUENCY OISPIACEMENT OF FOCUSING FFTI 

• jf • «n'/n 


Figure 41 

Coefficient Selection 


2. 3. 3. 2 Subarray Simulation Results 

The first issue to be examined in the step transform subarray algorithm 
is the amount of overlap required between successive subarrays. The 
results of the simulation of two overlap cases is given in Table 14. 

Case 1 is an overlap of 4 to 1 and case 2 is an overlap of 2 to 1 . 

That is, case 1 represents a subarray spacing of one fourth of the 
subarray length and case 2, a subarray spacing of one half its length. 

The column labeled 'COEFF' is the FFT coefficient output by the focusing 
FFT (second FFT). The parameters for the second or focusing FFT were 
chosen so that both cases had the same number. The 33rd coefficient 
in column 'COEFF' corresponds to the DC component of the good coefficients 
and the 65th the last valid coefficient. These two coefficients 
represent the extremes with the other 32 coefficients lying between. 

As can be seen from the 'ISL/IML' columns, the greater the overlap 
or subarray sampling rate the better the performance. Increasing the 
overlap increases the spacing between the main subarray lobe and its 
grating lobe caused by subarray sampling. However, higher overlap 
ratios also increase the required amount of computational effort. 
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Averaging over the 64 coefficients for case 2, the average per- 
formance Is found to be 25.6 db for 20 ISL/IML which Is well below 
the 20 db required. 


Table 14 
Subarray Overlap 

ID 20 


CASE 

COEFF ISL/IML 

PSL/ML 

ISL/IML 

PSL/ML 

1 

65 3ao 

40.7 

31.9 

40.8 


33 37.0 

AZ3 

37.0 

4i5 

2 

65 21.9 

19.0 

18.9 

16.7 


33 3ai 

37.1 

31.3 

37.1 



CASEl 

CASE 2 



1st FFT 

• 64 

64 



SUBARRAY 

SPACING 

• 16 

32 



# SUBARRAYS 

• 64 

32 



2nd FFT 

• 256 

128 



*BOTH FFTs HAMMING WEIGHTING 
•LOOK 2 

•AVERAGE ISL/IML FOR CASE 2 OVER 64 
COEFFICIENTS 25.6 dB 

The effects of the grating lobes on the subarrays can also be reduced 
by a proper choice of the weighting function across the first FFT. The 
window is chosen to place a null near the grating lobe. Table 15 
shows that a 30 dB Taylor weighting will meet the requirements. 

The table also shows that optimum performance with the first FFT 
weighting function is achieved by a careful balance between mainlobe 
width, which affects the amplitude of the output grating 
lobe, and sideiobe level which sets overall ISL/IML. 
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Table 15 


Effect of Weighting of First FFT 
in Subarray Process 




2D 


TYPE 

COE^ 

ISL/IML 

PSL/ML 

HAMMING 

65 

18.9 

16.7 


64 

18.3 

17.3 


34 

31.2 

37.1 


33 

31.3 

37.1 

25dB TAYLOR 

65 

15.5 

22.9 


64 

15.6 

3.1 


34 

34.3 

41.5 


33 

34.8 

41.9 

30dB TAYLOR 

65 

21.3 

24.0 


64 

21.4 

3.6 


34 

32.7 

38.4 


33 

3a7 

3&1 


’30dB TAYLOR WEIGHTING WILL MEET THE REQUIREMENTS 
2.3.4 Performance Comparisons 

Having chosen the parameters ‘<hich will meet the requirements for 
the different algorithm the relative performance between algorithms 
is shown in Table 16. Table 16 shows the performance with and with- 
out the addition of linear range migration to the quadratic range 
migration. The effects of interpolator complexity is also summarized 
in Table 16. For the FFT/Convolver approach a 2-point interpolator 
would seem to be adequate since 13 db of 18 db of the processing 
improvement compared to no interpolator has been obtained. For the 
subarray approach a 3-point interpolator appears adequate. The in- 
tegrated sidelobe values in Table 16 for the subarray system have 
been averaged over all valid coefficients out of the focusing FFT. 
Table 17 shows that doubling the amount of linear range migration does 
not affect the performance for the FFT convolver algorithms. 

The reason that the subarray case in Table 16 with linear migration 
and no interpolator provides better performance than the quadratic 
only case is that the linear migration tends to offset the quadratic 
in look 2. 
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Table 16 

Linear Range Migration (24 cells) 

• LiNEAR RANGE MIGRATION OF 24 CELLS 

ADDED TO QUADRATIC 17 CELLS) 



LINEAR 


INTERPOLATOR POINTS 



MIGRATION 

0 

2 

3 

4 

FFT 

NO 

12.4 

2a3 

29.9 

32.1 

CONVOLVER 

YES 

12.4 

25.8 

27.6 

30.5 

SUBARRAY 

NO 

12.5 

- 

24.6 

as 


YES 

13.6 

19.1 

24.4 

25.4 


DATA TAKEN FOR LOOK 2 


Table 17 

Linear Range Migration (40 cells) 

• ISL/IML WITH A LINEAR RANGE MIGRATION OF 
40 CELLS ADDED TO QUADRATIC (7 CELLS) 


• FFT/CONVOLVER APPROACH 


2D 

LOOK ISL/IML (PSD 

1 -12.4(19.7) 

2 -12.4(19.7) 

3 -12.4 i 19.8) 

4 -12.4(19.7) 

1 -25.6(30.7) 

2 -25.8(31.3) 

3 -25.8(31.4) 

4 -24.8(29.5) 

1 -29.9 (37.7) 

2 -30.5(38.5) 

3 -30.5(37.3) 

4 -28.2(35.1) 


INTERPOLATON 

POINTS 

0 


2 


4 


71 


2.4 Technology Survey 

Technology Is an important consideration in the design of the ADSP. 
The three principal areas designated; digital integrated circuits, 
digital architecture and software encompass the major design issuer 
related to the system. This section summarizes the study .evaluation 
and projection of these technology issues for the ADSP and future 
SAR systems. 


2.4.1 Digital Integrated Circuits 

2. 4. 1.1 Survey of Integrated Circuit (IC) Manufacturers 


Keys to the determination of IC technology available for the ADSP are 
the current and future plans of semiconductor manufacturers. An 
industry survey was conducted to determine plans and projections in 
the 1983-85 and post-1985 time frame. The companies included were: 


Advanced Micro Devices 

Fairchild Semiconductor 

American Microsystems 

Intel 

Intersil 

Motoifola 

Mostek 

National Semiconductor 


Raytheon Semiconductor 

RCA Solid State 

Signetics 

Texas Instruments 

Harris Semiconductor 

Synertek 

Monolithic Memories 


Information relative to technologies, memory size and configuration, 
LSI-VLSI logic functions, speed, power and cost were requested for 
the 1983-85 time frame. For post 1985 projections on technologies, 
line widths, gate density, speed, and power were reque.sted with 
specific projections on memories, arithmetic functions and micro- 
processors. 

Several of the companies surveyed declined to respond due to company 
policies on the release of long range plans. A summary of the results 
is contained in Trhle 18-20 and includes current technology, firm 
plans -1983-1985, projections 1983-85, and post 1985 projections. 
Random access memcies, EPROMS/PROMS/ROMS, and IC logic functions 
are included. 
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2. 4. 1.2 IC Technology Trends 

The IC technologies in current use span a large gamut of semi- 
conductors. These include the high speed emitter coupled logic (ECL) 
line to the low power, relatively low speed complementary metal oxide 
semiconductors (CMOS) and integrated injection logic (I^L). Bipolar 
transistor-transistor logic (TTL) with Schottky and low power Schottky 
(STTL, LSTTL) versions are the dominant semiconductor technologies now 
used in digital logic. N-channel MOS (NMOS) is the dominant large 
memory technology and silicon on sapphire (SOS) is used primarily for 
specialized military and space applications. 

As digital logic moves to shorter channel lengths we can expect 
changes in the technologies used. Although we can expect the same 
technologies to be in general use through 1985 we will see a definite 
shift toward narrow channel CMOS replacing TTL and ECL because of its 
low power, with high speed, capability. 

In the post 1985 time frame the principal technologies will be NMOS and 
CMOS with CMOS/SOS possibly moving into some general corranercial 
application. 

There is general agreement that CMOS will have emerged as the prime 
digital semiconductor technology because of its inherent low power 
dissipation. As line dimensions of CMOS circuits decrease, the speed 
also increases. CMOS circuits have been fabricated with 2 micron 
dimensions giving average gate delays of 1 nsec. When projecting 
up to 50,000 or more gates per chip, the power dissipation per gate 
must obviously be low to prevent thermal breakdown. At line 
dimensions of 1 to 2 microns and less the speed and power advantages 
of CMOS/SOS over CMOS are diminished, particularly when using oxide- 
isolated CMOS technology. 

Although CMOS has been indicated to be a post 1985 technology, it 
could happen much faster than that. There is a growing realization 
of the speed-power advantage of advanced CMOS in the industry-and 
increased investment. When the cost of the CMOS circuits meets the 
NMOS ’chnology costs, it will definitely be the choice because of 
its lower power which simplifies vhe power supply design, cooling 


73 



requirements and packaging system. 

The impact of the post 1985 technology on SAR processing will be in the 
feasibility of on-board processing provided by low power, dense 
memories, and VLSI with low power per gate provided by CMOS. 

Gallium arsenide (GaAs) can also be expected to emerge as the high 
speed technology In the post 1985 time frame. Research efforts with 
GaAs have produced GHz functional logic elements and A/D converters. 

2. 4. 1.3 Random Access Memories 

Random access memories have been broken down into dynamic and static 
categories. It can be noted that for the ADSP, static RAM's are 
preferable since the synchronous nature of input and output data 
is more easily controlled with static RAM storage. In general, 
dynamic RAMs are less costly on a per-bit basis. Current static RAM's 
are much faster, consume more power and, in accordance with general 
trends, cost more per-bit for the higher speed capability. 

Table 18 

Random Access Memories 




CURRENT 

TECHNOLOGY 

1983-1985 

FIRM 

1983-1985 

PROJECTIONS 

POST 1985 
PROJECTIONS 

DYNAMIC 

SIZE 

T0 64K 

64K 

128X-512K 

64K-1M 


CONFIGURATION 

xl 

x8, x4, xl 

x8, xl6 

- 


ACCESS (nsec) 

100-300 

200 

50-200 

50-100 


POWER (mW) 
(ACTIVE/ STANDBY) 

300/20 

300/20 

“ 

* 


COST/ BIT («) 

0.1 

0.04-0.1 

0.02-0.1 

• 

STATIC 

SIZE 

T016K 

16K-64K 

64K-256K 

64K-1M 


CONFIGURATION 

xl, x4, x8 

x8, x4, xl 

x8, xl6 

- 


ACCESS (nsec) 

50-100 

25-200 

50-200 

50-100 


POWER 

600/100 

200-1.000 

- 

- 


COST/ BIT («) 

0.1-0.4 

0.1-0.3 

0.02-0.1 

- 
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The 1983-85 projections generally continue to follow the "Moore 
curve" of memory development first observed by Gordon Moore of 
Intel Corporation. It shows memory capacity per chip doubling 
every year. Perhaps more significantly; the cost per bit is pro- 
jected to drop significantly, at least by some companies, in this 
time frame. In addition, the memory size firm committments of 
manufacturers is probably conservative. Figure 42 shows the state 
of the art for developmental memory circuits which is currently 
(mid 1981) 256 K bits for bulk CMOS. There is about a two year 
delay between development and commercial introduction. Based on 
the Figure 42 curve, we should see megabit dynamic memories on the 
market ii 1985. 



Figure 42 

Memory Chip Capacity Trends 
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Memory cost trends are shown in Figure 43. Both chip costs and 
installed memory costs are included. Chip costs should be close 
to the 0.01 cents per bit level in 1985. This translates to about 
$50,000 for the chip costs of the corner turning memory of the ADSP 
processor in the 1985 time frame. 



Figure 43 

Memory Costs/Bit Trends 
2. 4. 1.4 EPROMS, PROMS, ROMS 

An essential part of the ADSP processor are the programmable Memories 
used in the control firmware of the processor. Projection of circuit 
capabilities for erasable programmable read only memories (EPROMS), 
programmable read only memories (PROMS) and mask programmable read 
only memories (ROMS) are given in Table 19. 
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Table 19 


EPROMS, PROMS, ROMS 




CURRENT 

TECHNOLOGY 

1983-1985 
FIRM PLANS 

198M985 

PROJECTIONS 

POST 1985 
PROJECTIONS 

EPROMS 

CONFIGURATION 

2K 4K, 8M 
WORDS BY 
2, 4. 8 BITS/ 
WORD 

32K. 64K 


TOIM 


ACCESS (nsec) 

350-600 

200 




COST/ BIT (({.) 

0.2-0.4 

0.1-0.3 



PROMS 

CONFIGURATION 

1 K, 2K WORDS 
BV 1. 4, 8 BITS 
PEri Vl/ORD 

T0 8KX8 


TOIM 


ACCESS (nsec) 

45-100 

50 




COST/ BIT ((L) 

0.3 

0.1-0.3 



ROMS 

SIZE 

T0 64K 

64K-128K 

512K 

TOIM 


EPROMS are fabricated with MOS technology and although there are 
a wide variety of configurations available, they are generally too 
slow for convenient application to the ADSP. PROMs on the other hand, 
are made with bipolar technology and are much higher in speed although 
not as large. However, the maximum size PROM, 16K in up to a 2K X 8 
configuration is applicable to an FFT signal processor. 

A mask progranmiable read-only memory chip (ROM) may be applicable to 
the ADSP. The size and speed of ROMs of up to 64K are applicable, 
and if the set-up charges are sufficiently low to be offset with the 
quantities required in the ADSP it is the preferred direction. 

2. 4. 1.5 IC Logic Functions 

IC logic function projections are summarized in Table 20. Particularly 
noteworthy to the ADSP (and all other signal processing application) is 
the expected emergence of various "FFT chips" in the 1980' s. These 
circuits will make the principal computation requirement in the ADSP 
much less costly to meet. 
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Table 20 

IC Logic Functions 


CURRENT 

TECHNOLOGY 


198M985 
FIRM PLANS 


1983-198S 

PROJECTIONS 


POST 1985 
PROJECTIONS 


• 16 TO 32 BIT 
MICROPROCESSORS 

• 10,000 GATES /CHIP 


OTHER 

• 16 BIT ERROR DETECTION 
AND CORRECTION 

• 8 BIT COMPARATOR 

• PROGRAMMABl£ LOGIC 
ARRAYS 


ARITHMETIC 

MULTIPLIERS 
8 X 8^ 50 nsec 
16 X 16, 100 nsec 

ALUs 

(1981) 16 BIT FLOATING 
POINT AND PROGRAM- 
MABLE ALUs 


• 8T016BIT 
MULTIPLIERS 


• FFTCHIP 

t 1,000-10^100 
GATES/tHIP 

• 5-20MHZCLOCK 
CLOCKS 


• 32 BIT MULTIPLIERS 
AND MICRO- 
PROCESSOR, 20MHz 

• FFT P;10:ESS0R ON 
CHIP 

• 50,000 GATES/ 
CHI#» 


Most of the new logic circuit development in recent years has been 
directed toward microprocessor development. High speed LSI 
arithmetic circuits applicable to signal processors have not been 
developed to the extent possible. TRW has marketed a number of 
LSI arithmetic components which meet signal processing needs. Chief 
among these are a line of parallel multipliers of up to 24 bits. 

The 16 X 16 bit multiplier has a speed ov about 100 nsec. TRW is also 
planning to release a floating point adder circuit by early 1982 
which could be a key component in the ADSP. 

Another circuit of special interest to the ADSP is an error detection- 
correction circuit. Because of the large memory in the ADSP, it may 
be desirable to employ error detection and correction to increase re- 
liability. 

2. 4. 1.6 VHSIC 

Perhaps the most significant development effort now underway for the 
future of signal processing technology is the Department of Defense 
(DOD) VHSIC program - for very high speed integrated circuits. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 

The product of the number of gates on an IC and the clock frequency 
is a figure of merit for integrated digital logic that is particularly 
useful in assessing computing or signal processing power. The long 
term goals of the DoD VHSIC program is to achieve a gates-times-clock 
rate of 3 xlO . This goal is compared to a number of advanced 
technologies in Figure 44a. All of the technologies shown are or were 
candidates for achieving the VHSIC goals. 



QtAt TED - GALLIUM ARSENIDE TRANSFERRED ELECTRON DEVICE 

GtAt FET - GALLIUM ARSENIDE FIELD EFFECT TRANSISTOR 

ISO - ECL - DIELECTRICALLY ISOLATED EMITTER • COUPLED LOGIC 

'CMOS/SOS - COMPLEMENTARY METAL OXIDE SEMICONOUCTOR^ILICON ON SAPPHIRE 

SO - TRIFLE DIFFUSED BIPOLAR 

OCCL - DIGITAL CHARGE COUPLED LOGIC 


Figure 44 a 

Comparison of New Technologies on 
Clock Rate - Gate Basis 

The VHSIC program objectives include a shrinkage of the circuit 
geometric dimensions from 5 microns to 1.5 microns over 3 years 
and to 0.5 microns over 6 years. The scaling theory for circuit 
performance indicates that this size shrinkage will provide lower 
power-delay products. For example, in CMOS circuits, the basic 
geometric performance parameter is the transistor channel length. 
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2 

L_. With power dissipation proportional to L and the delay to 
c c 

1/L-, the power-delay product is thus proportional to L . 
c c 

The projected VHSIC circuit capabilities will permit chip imple- 
mentation of advanced signal processing functions un a chip. All 
contractors selected for the VHSIC Phase I programs have objectives 
which are oriented toward signal processing. 

2.4.2 Digital Architecture 

An overview of digital architecture approaches was developed in relation 
to the functional characteristics of the ADSP. These functional 
characteristics are surranarized in Table 21 for the range and azimuth 
processing functions. 

Table 21 

ADSP Functional Characteristics 


RANGE PROCESSOR AZIMUTH PROCESSOR 


NATURAL PARTITIONING 


PROGRAMMABILITY 


BULK STORAGE 


MAJOR FUNCTIONS 


• RADAR PRI 

• PULSE LENGTH 

t BANDWIDTH 

• TW PRODUCT 

• WAVEFORM 

• WEIGHTING 

• RANGE SWATH 

• NONE REQUIRED 


• RANGE CELL 

• LOOKS 

• STORAGE RANGE MIGRATION SPAN 

• PRI 

• AZIMUTH COMPRESSION 

• RANGE MIGRATION 

• CLUTTER LOCK 

• FOCUS 

ORANGE SWATH x AZIMUTH X 
COMPRESSION X LOOKS x 
OVERUP FACTOR 


• FFT 

• COMPLEX MULTIPLY 

• AMPLITUDE WEIGHTING 

• INTERPOLATION 


• RANGE PROCESSOR FUNCTIONS PLUS: 

• CORNER TURNING 

I RANGE MIGRATION COMPENSATION 

• DECIMATION 

• MAGNITUDE. SCALING 
I INTEGRATION 

• DOPPLER ICLUTIER) ESTIMATE 

• AUTO-FOCUS 


The first functional feature, natural partitioning, refers to par- 
titioning approaches based on the synthetic aperture radar processing 
algorithms or the inherent radar operation. For example, a natural 
partition in the range processor is the radar PRI (pulse repetition 
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interval) during which a single radar pulse is transmitted and its 
reflected signals from the range swath coverage are received. 

Programmabi 1 i ty must be incorporated for the range and azimuth pro- 
cesses as indicated in Table 21. A major requirement is bulk 
storage which occurs primarily in the azimuth processor. The size 
of the bulk store, required for range-azimuth corner turning, is 
determined by the product of the range cells and pulses coherently 
integrated with the number of looks and signal overlap also invol- 
ved. It is also dependent upon the algorithm employed. 

The major functions listed in Table 21 constitute the well defined 
functions for range and azimuth SAR processing. A programmable 
processor capable of performing all of these functions will by its 
nature have functional capabilities not listed. Modifications in 
processing algorithms should therefore be possible. 

Several basic architecture approaches were investigated for the ADS°: 
pipeline, parallel, SlMU (single instruction stream-multiple data 
stream), cross-bar, and various multi-bus architectures. 

2.4.2. 1 Pipeline Processor 

A pipeline processor using the FFT convolver processing algorithm 
is shown in Figure 44b. It has the primary feature of minimizing 
hardware both memory and arithmetic processing elements. All 
elements of the pipeline operate simultaneously and the control is 
fairly straight forward. It can be programmed for changes in the 
waveforms or focusing function via the spectral reference functions 
and the FFT size. The arithmetic and memory can be scaled to match 
the word size as the signals progress through the pipeline. However, 
this approach is not usually followed since net cost reductions are 
generally realized if fewer design variations are used. Techniques 
for reliability enhancement include providing for interchangeable 
pipeline segments and the incorporation of extra FFT stages in the 
pipeline which can be switched into use when a failure occurs at 
any particular stage. Growth to higher capacity systems must be 
done by the addition of complete, parallel pipeline systems. 
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• MEMORY; CTM ♦ RMC ♦ MK 

• FR STAGES; 46 

I PROGRAMMABILITY; COMPRESSION RATES, PIPELINE ADJUSTMENTS 

• RELIABILITY ENHANCEMENT; SWITCHED FFTs OR FFT STAGES 

• MODUURITY; PARALLa SYSTEMS 


Figure 44b 

Example ADSP Pipeline Processor (FFT Convolver) 

2. 4. 2. 2 Parallel Processor 

An approach to implementing the ADSP with parallel processors is 
shown in Figure 45. In this case, a total of 52 programmable FFT 
processing units are used to compute the required FFT computations. 

In the range correlator, each unit is capable of processing a com- 
plete single pulse range scan. Successive pulses are inputted to 
successive processors and the total number is set to operate on the 
total processing load at an internal clock rate of 8 MHz. The azimuth 
correlator is different in that each parallel unit handles a fixed 
segment of range cells. This permits the bulk corner turning memory 
to be divided equally and allocated to the respective parallel 
azimuth processors. An FFT convolver function is inherent in the 
azimuth correlator as in the range correlator, but because of the 
spectral domain range migration compensation, the forward and 
inverse FFT's are separated. The range migration compensation unit 
is then a special processor in between the forward and inverse FFT. 
The advantage of the parallel processor is in the standardization of 
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hardware and high reliability. The current thrust of L5I technology 
permits Implementation of a FFT processor at reasonable cost levels, 
which when constructed In quantity will further reduce the unit cost. 
Spare units can be err^loyed to enhance reliability. 



• MEMORY: CTM 4- I5RMC >238K 

• PROCESSOR UNITS: 52 

• PROGRAMMABILITY: AT PROCESSING UNIT LEVEL 

• RELIABILITY ENHANCEMENT: SPARE PARAlia PROCESSORS 

• MODULARITY: ADDITIONAL PROCESSORS 

Figure 45 

Example ADSP Parallel Processor (FFT Convolver) 

2. 4. 2. 3 Single Instruction - Multiple Data (SIMO) 

The SIMD approach shown in Figure 46 Is another parallel processor 
approach with an important difference. All of the proc.?ssors operate 
in-step on the entire processing algorithm. This approach offers a 
cotimon control approach, but it is not good for the ADSP. Since 
each unit must do the entire process, it must have a full corner 
turning memory or at least a significant portion of the range swath. 
Thus, the total memory requirements of the system are unreasonable. 
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I MEMORY; -25CTM 

• PROCESSING UNITS: -50 

• PROGRAMMABILITY: COMMON CONTROL 

I RaiABILITY ENHANCEMENT; PARALLEL UNITS 

• MOOilLARITY: ADDITIONAL PROCESSORS 

Figure 46 

SIMD (Single Instruction - Multiple Data Stream) 

?. 4.2.4 Cross-bar 

Cress-bar architecture has long been used in high performance digital 
processors. This approach, shown in Figure 47, can be used to form 
a programmable pipeline signal processor. It has particular advantages 
in cases where several different processing units may be used in a 
varying sequence. The need for a variation of processors would arise 
in a signal processing application where computationally intensive 
algorithms must be performed. It may then be cost effective to 
construct special purpose processors for these functions. The number 
and types of processors might vary for each installation. A technical 
problem with the cross-bar approach is the design of the cross-bar 
switch itself. Several methods have been suggested in the literature 
for constructing modified variations of the cross-bar. These approaches 
are aimed at matching the real switching needs with an efficient 
switching mechanism, generally a multi-level switch matrix not unlike 
an FFT flow diagram. However, all of the proposed switching schemes 
do not provide the flexibility of a cross-bar and thus limit the overall 
efficiency of a general purpose processor. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 



$ «CM(NIV> K«cm 

f n^SSINOIMITSi MIXWMMuaANDI>IPQ.l«SUISyS1EMS 

» controi aw seoudscwpihkcssins units 

I RaMllftllVCNHAiaMDir. SPA« PROCCSSINC (HITS 

• MOOUMRIlVt PARAUa SVS1EMS 

Figure 47 

Cross- liar Architecture 
2. 4. 2. 5 Multi -bus Architecture 

A multi -bus architecture can be used and Figure 48 does not do 
justice to the many variations which are possible. Recent efforts 
in this area have been concentrated on heirarchical schemes which 
partition the hardware structure into various levels of operational 
control. The control problems of this approach seem to be overly 
complicated for a signal processing system as essentially well 
defined as the ADSP. 

STRUCTURE ~ VARIABLE - 2. 3 DIMENSIONS 

0 +- 


I MEMORY: > MINIMUM STORAGE 

• PROCESSING UNITS: >S0 

• PROGRAMMABILITY: SEOUENCE. SIGNAL PATH CONTROL 

• RaiABILITY ENHANCEMENT: AOAPTIW FUNCTIONAL PATHS 

• MODULARITY: ADDITIONAL PROCESSORS 

Figure 48 

Multi -bus Architecture 
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2. 4. 2.6 Progranmable Signal Processor Developments 

Several government- Indus try efforts to build high performance, 
programmable signal processors are underway. These developments 
are not near the mature stage where one can buy a system. A brief 
suiranary of these developments follows. 

Advanced On-Board Signal Processor (AOSP) 

AGENCY: DARPA 
CONTRACTOR: RAYTHEON 

OBJECTIVE: DEVELOP SIGNAL PROCESSING SYSTEM FOR FUTURE SPACE 
RADAR, ELECTRO-OPTIC AND COMMUNICATION SYSTEMS CAPABLE OF 
MEETING ON-BOARD SIZE, WEIGHT AND RELIABILITY REQUIREMENTS 
ARCHITECTURE: MULTIPLE BUSSES INTERCONNECFING ARRAY COMPUTING 
ELEMENTS (ACEs) 

CLOCK RATE: 10-40 MHz 
MEMORY/ACE: > lOM BITS 
GATES/ACE: > lOOK 
ACES/ADSP: 56 

STATUS: STUDY/SIZING/SIMULATION/SOME BREADBOARDS UNDERWAY 
S-1 Super Computer 

DESIGN: LAWRENCE LIVERMORE LABORATORY 

SUPPORT: OFFICE OF NAVAL RESEARCH 

OBJECTIVE: TO DESIGN AN ADVANCED VECTOR PROCESSING 

ARCHITECTURE: CROSS-BAR CONNECTING MEMORIES WITH 16 HIGH SPEED 

PROCESSING UNITS (A-BOXES) 

CLOCK RATE: 57 MHz 

COMPUTATION RATE: 43P ;isec/4096 POINT FFT PER A-BOX 

A-BOXES/ADSP: 10 

STATUS: UNDER DEVELOPMENT 

Enhanced Modular Signal Processor (EMSP) 

AGENCY: NAVAL SHIP SYSTEMS COMMAND 

OBJECTIVE: TO DEVELOP AN ADVANCED PROGRAMMABLE SIGNAL PROCESSOR 
FOR NAVY APPLICATIONS TO REPLACE THE AN/UYS-1 
ARCHITECTURE: MODIFIED CROSS-BAR 
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COMPUTATION RATE: 80 M (REAL MULT I PIES //sec 
EMSPs/ADSP: -25 
STATUS: STUDY PHASE 

Other processors such as the MPP - Massively Parallel Processor 
(NASA-6oodyear) and MRP - Multi -mode Radar Processor (Hughes) are in 
various stages of development. The question that must be answered 
is whether any of these developments are practically applicable to 
the ADSP. 

The specialized needs of the ADSP can be best met with a processor 
matched to its requirements. This will also limit long term life cycle 
costs particularly in the area of mission software development. 

2.4.3 Software 

A study was completed of high order languages (HOLs) suitable for 
the executive and control software tasks within the ADSP host computer. 
The languages studied were Fortran, Pascal, SPL/I and Ada. Any one of 
the candidate languages is adequate for the task, but some present 
advantages over the others. It is desirable for documentation 
purposes and ease of understanding to have as much of the control 
software as possible written in the HOL, as opposed to assembly 
language subroutines. This gives languages which have real time 
control features an eoge, since a language without them must resort 
to assembly language for real time control. Another key feature 
for readability is the degree to which a language permits and 
encourages structured code. The newer languages, especially Pascal 
and Ada, are designed to force the programmer into structured 
pregrams. Manipulation of large data arrays is important and all of 
the candidate languages have this capability. Finally, the avail abi'^ity 
of the language is crucial. Is it supported on the selected ADSP 
computer? Is it sufficiently mature to be fully documented and 
readily usable? Is it in widespread use? These are the key questions. 

Table 22 sunmarizes the key characteristics of each language. Ada most 
closely meets the needs of the ADSP. This conclusion was reached after 
a careful comparison of the candidate HOLs. 
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The HOL which has been proven in almost every application is 
Fortran. In use fcr twenty years, Fortran is a suitable substitute 
for almost any language. The newest versions of Fortran include 
provisions for structured programming, including control constructs 
such as IF. . .THEN. . .ELSE. Most versions have fairly sophisticated 
I/O capabilities, although most are unique to one machine and not 
readily transportable to another. Most Fortran versions do not have 
any reel time control structures, leaving these functions to the 
opera MnvT system. Fortran can handle arrays of data, although 
complex file structures can be difficult to handle. Fortran is less 
readable than the nerf>r languages, although its widespread use may 
offset this, in that so many programmers readily understand it. In 
short, Fortran may be an acceptable vehicle for the task, but its 
limitations were precisely what the newer languages were designed to 
improve upon. It may pay to take advantage of these improvements. 

SPL/I, or Signal Processing Language I, is a relatively new language 
developed by Intermetrics, Inc. for Navy programmable signal pro- 
cessing applications. As a language designed expressly for signal 
processing, it is an inwiediate candidate for ADSP. It contains 
many features of use t ' this application, especially in the number 
and power of its real time control constructs. These can effectively 
handle multiple parallel processes and the interfaces between them. 

SPL/I shares the structured constructs and data typing of Pascal, Ada, 
and others. Detailed study of the language shows strong structural 
similarities to othijr HOLs, including many of their weaknesses. SPI/I, 
like Pascal and Ada, makes an effort to minimize dependence on a 
particular machine or compiler to guarantee software transportability. 

As such, I/O capabilities are essentially left to the implementation 
instead of being defined within the language. SPL/I has only four 
simple library functions defined to handle I/O. This is less than the 
other candidate HOLs. 

Another major problem with SPL-I is a lack of available implementations. 
Versions exist for Navy AN/UYK-7 and AN/UYK-20 computers, but outside 
the Navy, it is virtually unused. Much of this lack of use is due to the 
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expectations in the development of Ada, which has most of the same 
real time control features, and at the same tiiro provides for more 
readable code. 

Pascal, a language invented in 1970 by Niklaus Wirth of the University 
of Zurich, is intended primarily as a teaching language. The original 
intent was to provide the necessary language constructs to force a 
programmer to write structured programs. Many of the concepts de- 
fined by Pascal have been used in later languages, especially Ada and 
SPL/I. Pascal is a highly readable, easy to understand language. It 
is a strongly typed language, which means that all variables and 
constants in a program must be explicitly declared as to its type at 
the beginning of a program. Translation between different types is 
difficult unl?»ss their relationship is pre-defined. For example, 
this feature prevents a programmer from accidently equating feet and 
pounds in a calculation. Strong typing requires a programmer to 
think through the process he is defining before he actually codes It. 
Pascal incorporates many program control features which guarantee 
structured program flow: IF. . .THEN. . .ELSE and CASE statements provide 
explicit conditions for branching, making the reasons for each 
branch clear in each case. Pa*‘cal has array handling features which 
provide easy set up and manipulation of complex data structures such 
as linked lists. 

Because Pascal was originally intended as a teaching language, it has 
a number of shortcomings when applied to a sophisticated real time 
application. The I/O defined for Pascal in its original form is 
severely deficient, amd numerous implementations have attempted to 
improve upon its I/O capabilities. This has seriously inhibited 
Pascal's ability to be implementation dependent, and hence transpor- 
table from one compiler to another, and from one machine to another. 

It also makes it difficult to compare Pascal with other languages, 
since some versions of Pascal may have adequate I/O for the intended 
application. 

Another serious deficiency of Pascal is the absence of any real-time 
control structures. Pascal is not intended for real time applications. 
To be readily transportable from ore system to another, many Pascal 
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compilers, most notiibly the UCSO version, translate to a universal 
Intermediate level language (termed a p-code), which In turn Is readily 
translated to the -nachlne language of host computer. This two level 
translation process, while useful, limits the efficiency of a Pascal 
program In terms of the number of machine Instructions, It takes to 
execute a single Pascal Instruction. The two level process also 
makes the execution of a Pascal program very difficult to predict 
and/or minimize. 

Despite these limitations, Pascal has gained w1a.^ispread acceptance 
In the Industry. Some companies, like Texas Instruments, have 
standardized all of their software development around the use of 
Pascal as their preferred HOL. There are many versions available, 
and any final computer system chosen Is likely to have at least 
one commercially available compiler to use. In addition, an inter- 
mediate p-code language could be translated to any embedded 
processor's Instruction set, providing HOL capabilities and docu- 
mentation within the system hardware, provided real-time control 
and I/O problems can be dealt with. 

Many of the limitations of Pascal may be overcome through the use oi 
Ada, a derivative of Pascal first developed in 1979 by a design team 
from Honeywell -Bui 1 , under the direction of Jean Ichblah. Developed 
as a joint Army/Air Force project, Ada Is Intended to be the Depart- 
ment of Defense standard HOL for all major systems of the future. 

As such, Ada was designed to be an all-purpose language. Incorpora- 
ting the key features of Pascal together with real-time processing 
constructs and Improved I/O capabilities. Like Pascal, Ada Is a 
strongly typed language, with all variables defined clearly at the 
beginning of the program. Ada has all the same structured program 
constructs as Pascal, and In addition has features to simplify 
separate compilation of different sections of a large program. The 
I/O constructs In the original Ada definition are adequate, although 
still somewhat less sophisticated as most Fortran Implementations. 

This will probably be Improved when Ada Is actually Implemented. Ada 
is designed to be highly readable: the design Intent was to simplify 
documentation requirements by using a language as close as possible 
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to a Program Design Language (POL) which car. be used for documen- 
tation Instead of flow chart;s. Ada's multitasking capability 
provides clear definitions of multiple parallel processes, and 
their Interfaces. 

The primary problem with using Ada Is that currently Its Implementation 
Is In Its Infancy. The most notable development Is the Intel VLSI 
LAPX-432 microcomputer which Is designed for the Ada HOL. The Initial 
definition of the language has generated considerable Interest within 
the software Industry. A number of projects to develop working 
compilers has begun, both Inside and outside the defense Industry. 

RCA has Identified 24 such projects for various machines. It is 
expected that by 1983 most major computer systems will have usable 
versions of Ada, and that by 1985 Ada will be In widespread use. As 
such. It appears to be most promising candidate for the AOSP. Ada 
combines the best features of the other languages Into a single 
language: the versatility of Fortran, the real-time control of SPL/I, 
and the readability of Pascal. Table 22 sumnarizes the key characteristics 
of each language. 


Table 22 

Key Characteristics of Candidate HOLs 


FEATURC\lANGUAGE 

SPL-I 

FORTRAN 

PASCAL 

ADA 

READABILITY 

FAIR 

FAIR 

EXCELLENT 

GOOD 

M* 

structOres 

YES 


NO 

YES 

ARRA Y DATA 

Hlics 

GOOD 

FAIR 

GOOD 

GOOD 

STRUCTURED 

PROGRAMMING 

FEATURES 

YES 

IMPLEMENTATION 
DEPENDENT 
(FORTRAN 77) 

YES 

YES 

DATA TYPING 

YES 

NO 

YES 

YES 

I/O HANDLING 

POOR 

EXCELLENT 

POOR 

FAIR 

CURRENT USE 

AN/UYK-7. 
AN/UYK-20(NAVY) 
COMPUTERS ONLY 

WIDESPREAD 

WIDESPREAD 

FIRST 

COMPILERS IN 
DEVELOPMENT 

FUTURE USE 
(1983- ) 

UNPREDICTABLE 

WIDESPREAD 

WIDESPREAD 
(ALTHOUGH MAY 
BE REPLACED 
BY ADA) 

WIDESPREAD 
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2.5 ADSP SELECTED DESIGN 

2.5.1 Selection of Algorithm 

The FFT convolver was selected for the ADSP. The selection of a 
processing algorithm was made principally between the time domain, 
FFT convolver and subarray algorithms. The other approaches were 
considered in relation to these three. 

2.5. 1.1 Key ADSP Requirements 

The key requirements for the ADSP are summarized in Table 23. 
Continuous variation in processing parameters, with provision for 
modular construction and future growth are driving factors. In 
addition, the processor must accommodate continuous or burst modes. 



Table 23 




Key ADSP Requirements 

BASaiNE 

VARIATION 

• RANGE CORRELATOR: 

COMPLEX SAMPLE RATE 

IIMHS 



PULSE SAMPLES 

475 



PULSE COMPRESSION 

410 

10-6fi0 


RANGE SWATH 

4000 

500-4000* 


INTEGRATED RANGE SiDaOBES 

-ISdB 


• AZIMUTH CORmATQRs 

PRF 

2500Hz 



LOOKS 

4 

• 


AZLMUTH COMPRESSION RATIO /LOOK 

350 

20-350* 


INTEGRATED AZIMUTH SIOaOBES 

-20(16 



RANGE MIGRATION 

86 RANGE BINS 
14 LOOKS) 



* TRADE-OFF AMONG RANGE SWAIN. NUMBER OF LOOKS, AND RESOLUTION 


• MODUURITY - PROVISION FOR; 

• SWATH WIDTH - UP TO FOUR TIMES BASELINE 

• MORE PARALLEL LOOKS 

• MULTIPLE FREQUENCIES AND POURIZATIONS 


2.5. 1.2 Programmability 

The variation in parameters can be accommodated by the three 
major processing algorithms shown in Table 24. Both the time domain 
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and FFT convolver approaches offer full control over the SAR processing 
parameters s;uch as compression factor. However, the subarray Imple- 
mentation cannot be readily programmed to cover a continuous selection of 
parameters. The subarray formation process most easily accommodates 
a stepped parameter control. A disadvantage accrues to the time domain 
approach with variations in the number of looks. A variable prefilter 
would be difficult to implement in this case. The Type B time domain 
processor is used in the comparison because it is preferred over the 
Type A processor. 

Table 24 

Algorithm Programmability 

• Aa APPROACHES CAN ACCOMMODATE VARIATIONS IN RANGE SWATH 
THROUGH MEMORY CONTROL AND SCALING 


ALGORITHM 

vMAtiSH 

METHOD 

IMPACT ON HARDWARE 

TIME DOMAIN 
(TYPEBI 

COMPRESSION 

VARY NUMBER OF 
ACTIVE CORRELATORS 

MODEST INCREASE IN 
CONTROL FUNCTION 


LOOKS 

• SWITCH CORRELATORS 

• VARY PREFILTERING 

LARGE INCREASE IN 
CONTROL AND 
SWITCHING HARDWARE 

H CONVOLVER 

COMPRESSION 

CONTROL COMMANDS AND 
REFERENCE FUNCTION 

MODEST HARDWARE 
INCREASE AND CONTROL 
FUNCTION SOFTWARE 


LOOKS 

REPROGRAM CONTRA 
SEQUENCE 

MINIMUM HARDWARE 
INCREASE 

SUBARRAY 

COMPRESSION 

LOOKS 

VERY DIFFICULT TO OBTAIN 
CONTINUOUS CONTROL; 
CONTROL IN STEPS 

REPROGRAM CONTROL 
SEQUENCE 

URGE INCREASE IN 
COMPLEXITY OF 
CONTROL SYSTEM 


2.5. 1.3 Incremental Implementation and Growth 

Incremental growth or implementation can refer to variation in swath 
width, number of looks, or the number of processing cha.inels. In 
addition, the capability of the processor to handle data at real time 
rates is a factor. All of these issues are in a sense related. As 
the swath width and number of channels increases, the processor memory 
and computation can be expected to increase linearly. As the number 
of looks increases the computation requirements generally decrease and 
drop dramatically for the time domain processor with a prefilter. 
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The algorithm and architecture selection are related and in anticipation 
of the selection of a parallel architecture that implementation is 


assumed for the 

FFT convolver and subarray algorithm in Table 25. 


Table 25 


Incremental Growth and Implementation 




PARAMETER 


INCREMENTAL GROWTH 

AinORIIHM SWATH WIDTH 

NUMBER OF CHANNELS 

P»».U|LU»KS 

TIME DOMAIN 

Itypeb) 

ADD PARALLEL 
SYSTEM (S) 

ToTOr 

ADO PARALLEL 
SYSTEMS 

FFT CONVOLVER 

ticmimRi) 

ADD RANGE 
AND AZIMUTH 
PROCESSORS 

SIZE AND MEMORY 
STORAGE 


AZIMUTH 
SUBARRAY 
(PARALLEL . 
ARCHITECTURE) 

ADO RANGE 
AND AZIMUTH 
PROCESSORS 

IS MOST EFFECTIVE 
FOR LARGE 

COMPRESSION RATIOS 
(FEW LOOKS) 

ADD PARALLEL 
SYSTEMS 

INCREMENTAL IMPLEMENTATION 



TIME DOMAIN 
(TYPEB) 

NOT PRAaiCAL 
BECAl SF OF POOR 
MODUii MEMORY 
YRAOL 

• PROGRAMMABLE 
PREFILTEP 

• INCREASES REAL 
TIME RAI^ 

N/A 

FFT CONVOLVER 

<p;;rallel 

ARCHITECTURE) 

REDUCE NUMBER OF 
RANGE AND AZIMUTH 
PROCESSORS 

MORE LOOKS REDUCE 
SIZEOFFFTsAND 
MEMORY 

N/A 

AZDMUTH 

SUBARRAY 

(PARALLEL 

ARCHITEQURE) 

REDUCE NUMBER OF 
RANGE AND AZIMUTH 
PROCESSORS 

LOSES EFFECTIVENESS 
AS NUMBER OF LOOKS 
BECOMES LARGE 

N/A 


2.5. 1.4 Burst Multimode Processing 

The issues in burst mode processing fall into two categories; 1) the 
ability of the processor to accorranodate mode changes with different SAR 
parameters with no loss of imagery and 2) the ability to process special 
modes such as a large members of looks in a burst mode (See Figure 5). The 
first case is illustrated by Figure 49. A SAR processor will have an 
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Figure 49 

Burst Processing Delays 

inherent delay (tq) from input to output which is greater than the 
integration time (t^)« Radar bursts can switch instantaneously from 
one frequency to another as illustrated. A processor can have either 
a switched mode control or a traveling mode control. With a switched 
mode control, all of the parameters of the processor change together 
so the data within the processor when the mode swtich is changed is 
all lost. On the other hand, with a traveling control that moves with 
the data no additional invalid data appears beyond that equal to the 
integration time. The traveling mode control is obviously preferable 
and the time domain processor (Type B) is unique in its adaptability 
to a traveling mode control. The lost data can be avoided by increasing 
the bulk memory size by 50% for the FFT convolver and 100% for the 
subarray technique. 

The second burst processing problem can be considered in the extreme 
case where the length of a burst is equal to the integration time 
and the bursts are of arbitrary spacing. In this situation there 
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can be essentially one Image resolution element per look and the FFT 
convolver In Its normal form Is extremely Inefficient. The time 
dofTialn type B process Is well matched to the case because of Its 
mulltpller - accumulator structure. The type B processor can be 
used with each processing module handling a separate look for the 
full burst range swath. The FFT convolver and subarray systems can 
process this mode with some variation In their structure. The basic 
approach Is to deramp the Input signal, corner turn, FFT, correct 
for range all migration and Integrate the appropriate multiple looks 
as Indicated In Figure 50. 

^ LOOKS ► 
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BURST 
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FUNCTION 


FIGURE 50. INTERRUPTED BURST PROCESSING 
WITH MANY LOOKS 


96 


2.5. 1.5 Algorithm Cost Comparisons 

Relative costs were developed for the three principal ADSP algorithms 
and are summarized In Table 26. The design costs of the subarray 
technique overshadow Its lower materials costs, 'laterlals for the 
tirre domain processor reflecting the large number of computations 
required dominate Its cost profile. 

Table 26 

Relative Algorithm Costs 



ojliSyN 

FFT 

COI4VOLVER 

sub- 

array 

PROGRAM CONTROL 
MONITORING AND DESIGN 
ASSURANCE 

0.246 

0.207 

0.220 

DESIGN 

0.416 

0.384 

0.626 

AAATERIALS; SERVICES 

1.214 

0.263 

0.174 

ASSEMBLY AND TEST 

0.419 

0.146 

0.191 

TOTALS 

2.295 

lasi 

1.000 

(1.0) 

1.211 

(1.2) 


2.5. 1.6 Risks 

Development risks associated with the candidate algorithms are 
Indicated In Table 27. The time domains approach has problems with 
reliability and the control of a large hardware system at the nominal 
system clock rate of 11 MHz. The complexity of the subarray processor 
design Is its major risk factor while the FFT convolver algorithm has 
no serious drawbacks. Its only questionable area Is In achieving cor> 
venlent modularity of the hardware elements - a concern that also 
applies to the other algorithms. General algorithm-independent risks 
are also listed In the table and the major one affecting the lifetime 
of the ADSP Is meeting the correct mix between hardware and software 
controls. 
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Table 27 

ADSP Algorithm Development Risks 


SPECIFIC 


GENERAL 
APPLYING TO 
AUALGORITMIIHS 


TIME DOMAIN 


FR CONVOLVER 
SUBARRAY 


• CONTROL OF lARGE 
BUS S^TEML 
MULTIPLE CABINETS 

• U MHz CLOCK 

• RaiABILITY 


§ ACHIEVING MODULE 
PARTITIONING 


• ACHIEVING REQUIRED 
PROGRAMMABILITY 

« DEVaOPMENTOF 


• CONTROL SYSTEM 


TECHNICAL • 


ING CORRECT MIX 
JEEN HARDWARE 
AND SOFTWARE 


SCHEDULE: • TIMELY AVAILABILITY OF 
PURCHASED ICs 


COST: • 


UNPRfflICTABLE INFUTIONARY 
TRENDS OVER LONG SCHEDULE 


• MEETING COST PROJECTIONS OF 
ICs IF PURCHASES MUST BE 
SPLIT IN TIME 


2.5. 1.7 Summary and Selection 

A summary of the general characteristics of the various processing 
algorithms studied Is given In Table 28. Particular algorithms have 
features which may be useful In certain applications. For example, 
the subarray process, which offers a minimal hardware Implementation, 
Is advantageous for a single dedicated processor for low power and 
size on-board applications. However, the FFT convolver algorithm 
was seiected for the ADSP because: 

e Programming for multiple modes Is straightforward, 

• the higher cost of memory and computations (relative to 
subarray) are offset by savings in control system design 
costs, 

e a lower risk is associated with the approach, 

• integrated circuit technology thrust Is reducing the cost 
of memory and computations, 

e it lends itself to modular growth, 

• it provides a continuous selection of SAR parameters and, 

• the ADSP has no requirements for small size, low power or 
special environmental conditions. 
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Table 28 

Algorithm Summary Comparison 


n quaut 1 


OF 


FCOi 


1 

ALGORITHM 

MAJOR FEATURE 

DISADVANTAGES 

! FFT CONVOLUTION IlDI 

1 

EASIER TO PARTITION 
AND CONTROL 

NOT MINIMUM HARDWARE 

BLOCK FFT APPROACH 

NO CORNER TURNING 
MEMORY 

MORE MEMORY REQUIRED 

1 2D FAST FOURIER ! 

' TRANSFORM 

MOST FLEXIBLE 

1 

1 

HIGHER COMPUTATIONAL 
RATE DUE TO 2D OVERUP 

LIMITED BY DEPTH OF 
FOCUS 


2D FAST POLYNOMIAL | LESS MULTIPLICATIONS MORE DIFFICULT TO 

TRANSFORM | PARTITION 


SERIAL TIME DOMAIN • VERY FLEXIBLE REQUIRES MOST 

CORRELATION | EASY REFERENCE COMPUTATIONS 

GENERATION 
BEST FOR BURST MODE 


PARALLEL TIME DOMAIN 
CORRELATION 

NO ADVANTAGE OVER 
SERIAL SYSTEM 

REC. JIRES MOST 
COMPUTATIONS 
REFERENCE FUNCTION 
DIFFICULT TO GENERATE 

SUB ARRAY PROCESSING 

MINIMAL HARDWARE 

MORE DIFFICULT TO 
PROGRAM 


2.5.2 Selection of Architecture 

A comparison of pipeline and parallel architecture is given in Table 2C 
The oreponderance of key favorable factors make the parallel approach 
the optimum selection. A parallel structure is also characteristic 
of a number cf specific processors including the Massively Parallel 
Processor - MPP (19), the Advan^'ad On-Board Signal Processor - AOSP 
(20) and S-1 (21). Table 30 shows how these systems could be con- 
figured for the AOSP. Their use in the ADSP would depend upon the 
achievement of their development goals in a timely manner. In addition, 
they do not have the incremental implementation and growth characteristics 
desired. For those reasons we have recommended a parallel design 
tailored specifically to the ADSP requirements. 
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Table 29 

Pipeline - Parallel Architecture Comparison 


ADSP 

CHARACTERISTIC 

PIPELINE 

PARALLEL 

NATURAL PARTITIONING 

• EACH PIPELINE EUMENT 

doe:> own function 

• RANGE- PR 1 

• AZIMUTH - RANGE SLICE 

PROGRAMMABILITY 

• FRs - LENGTH CHANGE 

• CONTROL -TIMING MUST 
BE ADAPTABLE 

• FIRMWARE PROGRAMMABLE 
(NO HARDWARE RESTRUCTORING) 

MEMORY SYSTEM 

• CAN TAILOR EACH ELEMENT 
TO MINIMUM SIZE 

• MUST OPERATE AT HIGH 
SPEED 

• LOW SPEED OPERATION 

• MODULAR. CONVENIENT SIZES 

RELIABILITY/ 

ENHANCEMENT 

• MEDIUM RELIABILITY 

• EXTRA PIPELINE UNITS BUT 
CONTROL IS DIFFICULT 

• INHERENT HIGH RELIABILITY 

• SPARES ARE INDEPENDENT 
AND EASILY INSERTED 


• TESTING IS DIFFICULT 

• STATUS TEST IS STRAIGHT- 
FORWARD 

MODUURITY 

• MODUURHY REDUCES 
EFFICIENCY 

• MODULARITY INHERENT IN 
DESIGN 

RISK 

• A/IODERATE 

e LOW 


Table 30 

Processor Comparisons 
('FT Convolver Algorithm -- FFTs) 

AOSP 

MULTIBUS S-1 

ADSP MPP* (56 ACES- CROSSBAR 



REQUIRED 

PARALLEL 

est) 

(lOA-BOXES) 

REAL 

MULTIPLICATIONS 
PER SECOND -r 106 

2.240 

910 

2.240 

2.240 

REAL ADOS 
PER SECOND -i- # 

3.360 

4.428 

3.360 

3.360 

APPROXIMATE 

416 

16* ** 

560*** 



MEMORY 
STORAGE -f 106 
(BITS) 

* MULTIPLIES AND ADDS ARE NOT SIMULTANEOUS FOR MPP 

•• EXTENDABLE 

••• ASSUMES 5 MEMORY CHIPS PER ACE 
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2.5.3 Design 

The design selected for the ADSP features a parallel processor 
Implementation of the FFT convolver algorithm for both the range 
and azimuth correlators. It is programmable and adaptable to multiple 
missions. The parallel architecture allows a modular implementation 
suitable for expansion or reduced system capability. The parallel 
design achiever high system reliability and graceful performance 
degradation. 

The system is designed with TTL compatible logic, has an internal 8MHz 
operating clock rate, has floating point arithmetic and employs 64K 
dynamic RAMs for the bulk memory. It is controlled by a host computer 
for mission set up and performance monitoring, but because most of the 
control reqirements are embedded in firmware, a minimum of software 
effort is anticipated for a new mission set-up. 

Figure 51 shows the arrangement of the parallel processing modules. For 
the assumed 8 MHz clock rates, 22 parallel processors are used for the 
range correlator and 15 parallel channels are employed in the azimuth 
correlator. 

The architecture and algorithm permits a maximum degree of mission 
programmability with a minimum of softv/are set-up effort. The host 
computer is not involved in any of the real-time computations as a 
process proceeds buc is used for configuration and test functions. All 
real time controls are embedded in the signal processing elements. 
Commands from the host computer set up the SAR processing parameters. 

2.5.3. 1 Range Correlator Design 

The range correlator is illustrated in Figure 52. Each processing 
module handles a full radar PRI. The PRI/FFT unit, whose primary 
function is FFT processing, is the same module used for FFT processing 
in the azimuth correlator. 

Each PRI/FFT processor consists of an input/output (I/O) memory, a 
working memory for computations, an arithmetic computation unit, a 
reference function store, and a controller. The I/O memory is large 
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Figure 51 

Advanced Digital SAR Processor Design 




103 


Figure 52 
Range Correlator 



enough to hold the full range swath window prior to range compression. 

Data Is transferred from the I/O memory through the data switch to the 
working memory for processing. The computation unit can be prograiraned 
to do the FFT butterfly, weight ing, reference multiplication, magnitude 
calculation (I^+Q^ or ), normalization, scaling, and either real 

or complex FFT computations. 

2. 5. 3. 2 Linear Range Migration Correction 

From the point of view of control simplicity It appears preferable to 
compensate for all of the linear components of range migration separately 
from the quadratic terms. This Is the approach employed In the recommended 
design where the correction Is applied after range pulse compression. 

The Implementation Is simply an Incremental delay plus an Interpolated 
fractional sample delay. Our simulation showed that, with the FFT con- 
volver processing algorithm, adequate performance could be obtained with- 
out the use of a separate linear range migration unit If Interpolation 
were used In the frequency domain correction. However, we have Included 
the function since It may contribute to higher quality Imagery In some 
cases. 


2. 5. 3. 3 Corner Turning Memory 

The corner turning memory module functions are given In Figure 53. A 
PRI buffer Is used to capture the range swath to be processed. The 
net data rate through each parallel channel Is then only about 22/15 MHz. 
This low data rate permits the use of slow 64K dynamic RAMs for the bulk 
store. 

The organization of the corner turning memory modules Is shown In Figure 54. 
With a 4096 point azimuth processing window, the windows must be repeated 
each 2048 points. This is accomplished by using an overlap - save sequence 
in the memories. After the full 4K samples are stored In two memories the 
next 2K of data Is read In to one-half the memory while the 4K window Is 
read out. Operation In this manner means that the total bulk corner turning 
memory storage is approximately thv? product of the range swath and the 
azimuth correlation extent. 
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Figure 53 

Corner Turning Memory Function 
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Figure 54 

Corner Turning Memory Organization 
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2. 5. 3. 4 Azimuth Correlation 

After the first FFT is taken to convert the constant range lines to the 
spectral domain, the quadratic range migration correction (and any 
residual linear component) is applied together with the spectral reference 
focusing function. Bearing in mind that the data rate at the interface 
to this function shown in Figure 55 is less than 2 MHz, the implementation 
of the functions can be streamlined by time-sharing hardware elements. 

Most of the adaptive processes in the azimuth correlator are centered 
here. A separate clutter lock estimator is used. Although the focusing 
function generation is shown on the module, our selected approach for 
this function is to anploy a separate FFT processor module for 
computing all of the references as discussed later in this section. 



MULTI LOOK FFTs 


Figure 55 

Range Migration Correction and Matched Filter 
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2. 5. 3. 5 Mult Hook Integrator 

Figure 56 shows the design approach for the multi look Integrator. Input 
data Is first Interpolated, If necessary, for Image registration or geometric 
mapping correction, or changes In the output Image pixel density. If 
the autofocus Is not accomplished with modification to the focusing 
function, muUllook registration corrections can be applied here. The 
multiple looks are Integrated after the magnitude of the complex samples 
are obtained. The multiple looks are Integrated after the magnitude 
of the complex samples are obtained. The multi look menory store Is 
approximately one-eighth as large as the corner turning menory and all 
of the associated functions can therefore be Included on the module. 


MULTI LOOK 
MEMORY STORE 



Figure 56 

Multi look Integrator 
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2. 5. 3 . 6 Clutter Lock 

A clutter lock technique was described In (22). It relies on a measure 
of symmetry of the mean image energy of the four looks. The normalized 
error is determined from 

El + E 2 - E 3 - E 4 

En'TFTTrr|4 

where E 1-4 are the energy le.els in the respective looks. A modification 
of this approach is suggested here based on the following two observations. 
First, the referenced approach assumes that the actual mean image energy 
of a particular look, say look 3, does not fluctuate too much as the 
terrain is changed. Although this is generally valid, there may be specific 
peculiarities, e.g. transition from ocean to land, that will violate this 
assumption. Secondly, as the frequency offset (Fc-fc) approaches zero, 
the difference between the mean image energy looks 2 and 3 would appear 
to be insensitive to changes in (Fc-fc). 

The suggested modification uses two different segments of the range 
compressed data in such a way that the image of look 1 of the first segment 
coincides with the image of look 4 of the second segment. Furthermore, 
only the spectra of looks 1 and 4 are used to establish the sytranetry, as 
these are much more sensitive to changes in the frequency offset (Fc-fc). 

Two measures of symmetry can be used. The first is the difference of the 
total image energy of each look, i.e. E 1 -E 4 . Another measure is obtained 
by comparing each spectral line in the two looks. Let «|Si(k)^ and -^ 4 (k^ 
be the component of looks 1 and 4 respectively. Then a good measure^ is 

D |Si(k) - S4(k)l. 

This quantity appears to be much more sensitive to changes of the frequency 
offset (Fc-fc), as can be seen from the result of a typical simulation 
presented in Table 31. The disadvantage of using D to perform clutter lock 
is that one has to seek the minimum of D as (Fc-fc) is varied whereas the 
use of (E 1 -E 4 ) needs only to search for the zero crossing. 
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Table 31 

Clutter Lock Simulation Results 


• POINT TARGH 

t 1WO POINT INTERPOLATOR IN RANGE COMPENSATION ALGORITHM 
• YAW IS APPROXIMATaV 1/2 A LOOK 


• RESULTS: 




h’k 




MEASURE 

0.008 

0.006 

0.004 

0.002 

0 

-0.002 

-0.004 

El 

5.057 

5.055 

5.061 

5.070 

5.076 

5.082 

5.089 

U 

5.094 

5.087 

5.089 

5.083 

5.076 

5.070 

5.060 


-0.0364 

-0.0325 

-0.0286 

-0.0131 

-0.0004 

+0.0125 

+0.0287 

0 

3.989 

3.003 

1.964 

1.442 

0.343 

1.468 

1.985 
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2. 5. 3. 7 Generation of Focus Function 

The reference signal needed for the azimuth matched filter is of the form 

{e-ja(n-no)^}^ where the constant a depends on the range as well as other 
system paramei:ers. This dependence on range demands that new reference 
functions be generated fairly frequently for proper focusing, and the 
Fourier transform of the newly generated reference function be calculated 
if the matched filtering is carried out in the frequency domain. The 
required spectral focusing function can be generated by taking the FFT 
of a linear frequency ramp function which has been adjusted for the 
correct slope. While we have initially selected this first method, it 
is worthwhile to consider alternates which reduce the amount of computations. 
One method for reducing this computation task is described here. 

Assume that ^*j3(n-no)^ is the focusing function in use, and the value 
of a is changed to a+A for focusing at a different range. Rather than 
generate ^ Hn-no)^ anew, we rewrite ^ as follows: 

e-j(a-*-A )(n-no)2 = 

= [l-jA(n-no)^- ...j 

If is sufficiently small, the first two terms in the series expansion 
are sufficiently accurate for our purpose. In the transform domain: 

; {-ja(n-no)2/ ,-ja(n-no)2^- 

Thus, one needs to store only the three transforms on the right side 
of the equation and combine them using whatever value of a necessary to 
obtain the transform of the desired new reference. 


no 


To determine one range of A for which this approach offers sufficient 
accuracy, we generated a reference with A varying and used it for matched 
filtering the returns from a point target. The results were compared with 
those obtained when A-0. 

If there is no yaw, then the proposed alternate will produce less than 
O.SdB change in the response when A is kept within jfSO range cells. If 
a yaw corresponding to half a look is present, then A has to be within + 30 
range cells in order to keep the change to within 0.5dB. If the yaw is 
increased to one look, then A has to be within +20 range cells. 

While the alternate method offers significantly fewer computations, the 
availability of FFT processing modules makes the direct reference 
computation cost efficient. With a new reference required every eight 
range cells, two FFT processing modules can be applied to the task of 
computing the references for 15 parallel channels, 

2. 5. 3. 8 Control System 

System requirements which affect the range and azimuth processor controls 
are listed in Table 32. These general requirements impact in various 
ways on the processor depending upon the algorithm and architecture 
employed. 


Table 32 

AOSP System Control Requirements 

Rang« Processor 

Rsnga Cotnprasslon Ratio: 10 to 660 

Swath Width; SOO to 4000 


Aalinuth Procesior 

Azimuth Comprasslon: 20 - 350 per look 

Reference Function Update: Cross track - every 8 range cells 

Along track - every 12,000 radar pulses 


Required Modes 
Dual Frequency 
Dual Polarization 
Burst 
Conti nous 

Variations within Mission In altitude, look angle, resolution, 
number of looks, swath width, and azimuth weighting 
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The FFT convolver Is relatively easily controlled for the variation 
of system parameters. The length of the FFT can be changed to accommodate 
varying time bandwidth products and swath width. This is especially 
easily accomplished in a serial FFT where the length change is by firmware 
control. In addition, waveform and weighting function variations can be 
accommodated by either storing the additional functions in ROM or by 
inserting them into a temporary storage RAM via the control computer. 

A continuous selection of swath widths and waveform length can be 
batch processed with the FFT convolver within the limits of the maximum 
size FFT. The FFT can also be operated in a slidino aperture mode. 

Then successive FFT windows are processed at intervals of the difference 
between the FFT size and waveform length. Table 33 lists a number of 
variations in FFT and waveform size and resulting complex computations 
per output point. The table illustrates the logarithmic increase in 
the required computations as the waveform and FFT size increase. It 
also illustrates the computational advantages of increasing the FFT 
size relative to the waveform length. 

TABLE 33 

FFT Range Correlator Convolver Variations 

RANGE COMPLEX 

WAVEFORM COVERAGE COMPUTATIONS 
SIZE SIZE PER FFT PER POINT 


512 

256 

256 

20 

1024 

512 

512 

22 

2048 

1024 

1024 

24 

4096 

2048 

2048 

26 

8192 

4096 

4096 

28 

8192 

2048 

6144 

18.67 

8192 

i024 

7168 

16 

8192 

512 

7680 

14.9 
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Controls required on the FFT convolver to meet the total prografrarning 
requirements are indicated in Figure 57. Some of these controls may be 
most appropriately placed on a common control module for multiple range 
processors. The tradeoff in this regard is the consideration for 
multiple uses of the module, as in the azimuth correlator. 

Controls for the azimuth processor are distributed to each functional 
element. For the FFT convolver algorithm all required controls for the 
FFT are embedded in the FFT module used in the range processor. The 
primary control centers are then in the corner turn memory, the range 
migration-matched filter module and the multi-look integrator. Linear 
range migration is a separate function in the design as is the clutter 
lock estimator. 



SFT.iv ntiMimt 


Figure 57 FFT Convolver Control 


113 













ORIGINAL PAGE IS 
OF POOR QUALITY 


2.5. 3.9 Test Subsystem 

The complexity of the ADSP dictates the implementation of a diagnostic, 
status test subsystem in conjunction with the operating system. A 
conceptual test system is shown in Figure 58. The system elements 
would be exercised with a test signal sequence and check-sums would be 
obtained at selected locations in the processor. These check sums are 
read back to ‘'he control computer for failure localization. 

In addition to the built-in-test system, a test fixture should be developed 
for testing and debugging functional modules. This test unit could be 
used both for production testing and long term maintenance of the equipment. 

2.5.3.10 Physical Configuration 

The ADSP system was sized for its physical configuration. The module 
complement for the system is shown in Table 34. Large 15" x 17" modules 
were assumed capable of housing 200 or more standard dual-in-line circuits. 
Using an average of 170 circuits per module the total number of circuits in 
the system is about 33,000 including 5% spares. Five standard cabinets 
can be used to house the ADSP with about 2000 watts dissipation per cabinet. 



mn - FURCnONAkMOOUU 
ACC • AeeumiLATOfi 


Figure 58 Test System Concept 
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Wirowrap interconnections were assumed. Because of the large number of 
circuits in the system it is recommended that they be acquired with a 
burn-in specification of 60 hours, particularly if plastic circuits are 
used. A POP-11/23 computer with 256 K bytes of memory, dual disk, and 
hardware floating point was selected for the host computer although any 
of a large number of computers would be adequate. 


TABLE 34. MODULE SUMMARY 



FUNCTIONAL 

OPERATING 

MODULE TYPE 

UNITS 

SPARES 

PRI/FFT Prc essor 

54 

4 

Corner Turn Memory 

60 

4 

Range Migration Correction 

15 

1 

Multi-look Integrator 

15 

1 

Linear Range Migration 

1 


Timing Unit 

1 


Clutter Lock and Focus 

1 


Test Probe Interface 

6 


Controllers 

3 



OTHER SPARES 
(OPTIONAL) 

4 

4 

1 

1 

1 

1 

1 

1 

3 


156 10 17 


2.5.4 Impact of Technology 

The specific impact of future technology on the implementation of the ADSP 
can be deduced from Figures 59 and 60 which show key integrated circuit and 
architectural developments over the ADSP development time. All of the DC 
developments shown should serve to minimize the ADSP cost if they are 
available in a timely manner. Either 256 K bit dynamic RAM's or 64 K bit 
static RAMs will be available by 1983. Wide word configurations of up to 
8 bits will reduce costs. Particular logic developments which can be 
anticipated are highlighted by a coming selection of chips specifically 
designed to reduce the cost of implementing the fast Fourier transform (FFT). 
The culmination of logic will come with the advent of VHSIC technology which 
is primarily aimed at signal processing applications. In the microcomputer 
field the recent announcement of the i APX-432 Ada 32 bit microcomputer by 
Intel Corporation is just the forerunner of an expected set oi competitive 
VLSI computer systems. 
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Specific signal processors should also be completed within the general 
time frame of the ACSP as Indicated In Figure 60. These developments, 
however, as previously Indicated are not Immediately applicable to the 
AOSP unless they have sufficient maturity, which should come at about 
the time the AOSP Is completed. 




1983 


19SI 


198S 


RAMS: 


PROMS: 


LOGIC: 


AOSP 


4 

64K 2S6K 
0,giS DYNAMIC 
</BIT 


4 

8Ki8 

STATIC 


4 4 

FLOATING MISC 
POINT SIGNAL 

ADDER PROCESSING 
ICs 


4 

8Kx8 

4 

FFT 

AU 


VHSIC 

PARTS 


.4 


MICROCOMPUUR: 


INTEL 

4« 

(AOAI 


• AND 
ITS 

COMPCTITION 


Figure 59 

IC Technology Projection 
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16 CONTRACTORS! 
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Figure 60 

Signal Processor Architecture Projection 
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2.6 AOSP SCHEDULE AND COST FACTORS 

The nominal development time for the ADSP has been established as 
three years. However, changing the program schedule will affect the 
expected cost as illustrated in Figure 61 which was derived from outputs 
of the RCA PRICE* program for the ADSP design. However, at shorter 
schedules both costs and risks tend to increase until at about 18 months 
the risks are too high for rational consideration. A 36 month schedule 
provides a higher cost, but the additional time gives a higher success 
expectation. 


REUTIVE 

COST 



Figure 61 

Effect of Schedule on Cost 
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♦PRICE - Programmed Review of Information for Costing and Evaluation 
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If a two phase program is planned for the three year effort, the 
expenditure rate?; in Figure 62 will be obtained. The first 24 months 
of the job includes all of the design and first piece test work plus 
the integration of all the elements into a partial working system 
capable of doing all ADSP functions at a non-real time rate. The 
two-year design costs include all of the cabinetry and control for the 
full performance system. The final 12 month period then involves 
purchasing additional parts and modules and integrating them into a 
full performance systtsn. 



Figure 62 

Typical Expenditure Rates 
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2.7 IMPLEMENTATION PLAN ALTERNATIVES 

The recoiranended schedule for the ADSP is given in Figure 63. It reflects 
the two phase dev-elcpment outlined in Section 2.6. 


YEAR 1 YEAR 2 YEAR 3 



Figure 63 

Recommended Schedule 


The question of structuring the ADSP as a partial (non-real time) 
implementation with a four year development cycle was addressed. Full 
capability would be achieved with later add-ons. Alternatives for this 
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scenario are listed In Table 35. The first option Is essentially to 
stretch the two-year dei/elopment cycle of Figure 63 to four years. While 
this would achieve all program performance objectives, It would result 
In a cost Increase (not Including Inflation) of 20 percent. In addition, 
the four year program could produce a technology deficient design unless 
the hardware design details were delayed until the last two years. 

Table 35 

Alternative Trade-off 


OPTION 

1. NONREALTIME, 
FULL FUNCTIOiyiS 


REUTIVE 
COST 
(4 YEAR) 


0.78 


COMME NTS 

• WOULD DEMONSTRATE ALL PROGRAM OBJECTIVES 

• FOUR YEAR PROGRAM IMPOSES 20% COST PENALTY 

• FOUR YEAR PROGRAM MAY BE TECHNOLOGY DEFICIENT 


2. (a) NONREALTIME. 0.73 
REDUCED 

PROGRAMMABILITY 


• REDUCING PROGRAMMABILITY PJMINATES ONE OF 
KEY PROGRAM OBJECTIVES 

• TEST FACILITIES. BOTH BUILT-IN AND MODULE. 
ELIMINATED 


(b) CASE (a) PLUS 
REDUCED GROWTH 


0.60 • LOSS OF PROGRAM OBJECTIVE OF EXPANDABILITY TO 

FLU REAL TIME SYSTEM WITHOUT ANY DESIGN 
• FULL SYSTEM CAPABILITY WOULD REQUIRE 
REPACKAGING OF HARDWARE 


3. DELAY HARDWARE ? • TOTAL COST I S UNKNOWN -- COULD BE HIGHER THAN 

ALTERNATIVES 

• APPLIES MOST ADVANCED TECHNaOGY 

• HARDWARE DEVELOPMENT BETTER MATCHED TO 
CURRENT MISSION TIMING 


The second option In Table 35 would Involve reducing the level of 
programmability ?nd the built in test and module test facilities. 
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This would eliminate one of the key program objectives, programmability, 
and the elimination of test facilities would make the equipment very 
difficult to maintain. A substantial cost reduction (.40) is achieved 
if the previous two factors are eliminated plus growth capability. That 
is a partial systen, with its limited control, would be designed and 
built. A large disadvantage to this approach (2b) is that future 
expansion would require design effort and thus cause much greater ultimate 
costs. 

The final choice would be to delay the hardware development. While the 
ultimate cost of this approach is not known, it does have the advantage 
of applying the most advanced technology, either custom built or purchased. 
Also since the SAR missions which might use the ADSP do not occur until 
1986 and beyond, a delayed program would be more in line with mission 
timing. 
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3.0 CONCLUSIONS/RECOMMENDATIONS 
General Reconmendatlons 


1) The FFT convolver algorithm is the recommended choice for a programmable 
SAR processor not having specific environmental requirements. 

2) Parallel architecture is recommended since it provides modular 
structures ^fhich are matched to program needs in terms of incremental 
implementation and growth and the thrust of technology. 

Performance Results 

1) Equivalent performance can be achieved with the FFT convolver, 

step transform subarray, or time domain azimuth processing algorithms. 

2) The performance of the step transform subarray technique is a 
function of the subarray overlap factor. 

3) Range migration compensation without interpolation is not a 
recommended procedure. 

4) Adequate interpolation can be achieved with a four, three, and 
possibly a two point interpolator. 

5) Very little range spreading is observed with any processing algorithm, 
although some skewing of the resolution cell can occur in the subarray 
process. 

6) The integrated sidelobe levels are comparable with either the subarray 
or FFT convolver approaches. 

7) The mainlobe width (resolution cell) is determined by the weighting 
function as long as interpolation is employed. 

Algorithms 

1) The FFT convolver algorithm offers a high degree of programmabi 1 ity 
with moderate hardware complexity. 

2) Time domain approaches, while riroviding some good control features, 
require by far the most computations which translate directly to 
hardware size and cost. 
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3) Subarray processing offers a minimum hardware approach, but Its 
control syst^ complexity outweighs the hardware reduction for a 
highly progrunmable, laboratory based systan. 

4) The subarray technique may be most appropriate for a dedicated, 
onboard processor where size, weight, and power requiranents 
predominate. 

Architecture 

1) Parallel architecture offers the desired systan characteristics 
of modularity, modular implementation, modultr growth, high 
reliability, maintainability, programmability, and accommodation 
to new technology. 

2) Pipeline architectures minimize hardware, but are difficult to make 
reliable and do not lend themselves to modular implenentation and 
growth. 

3) Current government and industry advanced signal processor developments 
are not ready for immediate use in a real SAK processor. 

Integrated Circuit Technology 

1) VLSI, VHSIC technology trends will have a profound impact on signal 
processor implementations such as the AOSP in the 1985 time frame. 

2) Between 1981 and 1985 a continuing array of improved integrated 
circuits will become commercially available including: 256K dynamic 
RAMS, 8Kx8 static RAMs, 64K PROMs, FFT chips, and advanced microcomputers. 

3) These developments will make on board processors realizable with modest 
size and power. 
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4.0 NEW TECHNIQUES 

A modular parallel SAR signal processing design was conceived and 
designed which meets the severe performance requirements of the ADSP. 
The system, described in Section 2.5, is capable of being progranvned 
to cover a wide, continuous range of SAR parameters by virtue of its 
use of the FFT convolver algorithm. The algorithm has been implemented 
in a parallel architecture which allows incremental implementation 
and growth, high reliability, progrmmiability, maintainability, and 
low cost. The design is adaptable to subsequent advanced technology 
infusion. 
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APPENDIX - PERFORMANCE SIMULATION PROGRAMS 


PART I - SAP SIGNAL GENERATION 


ORIGINAL PACT- 
OF POOR QUALiTY 




Pt»0SH4>l TO 6eNe«4ft THE B*N6t WALK ANO 
THfc PH436 HI3TOWV OF A TAMGEI IN A 
SAB TYBfc ENVIHONMENT 


VIBTUAL At (4?<»6) >AG(a09b) 

virtual TEl'PAt3H«0) * tempo ( 3048) f ranger ( aa96) 
common /I',3TR/A(SP) 

c 

e 

CALL INST9 
CAUL INSTRA 

CAcL M«p (» I , TE mPA , r E mP8, R anger, LMN) 

CALL auTP(TbMPi, T empo, RANGER) 

CALL FFTR(AI , AO, TEMP A, Tempo, RANGER, LMN) 

ENO 


A1 


ooo oooorioooorioooooor*oooor»orioooor>oooooririorir>oor»rir>ooriooooooooor»oof# 


ORIGINAl. PAGE B 
OF POOR QUALITY 


SiiHQI.pOR 


instruction stT poa the s*« 
c IS TMe speed uf cxght ■ s, e « a ( m/scc ) 

PI ■ tt,a«T*Mt,) 


ten TRSNSNITTEO PMEQUENCV ( MMZ ) FT 

* 13 ) iiiavkCtNUTH OF F T C M ) C/FT L 

M 3 ) VtLOCTiT OF SP*CECH«FT ( HM/ 8 EC ) V 
*(«) RANSt TO TARGET ( SM ) R 
* 15 ) SAMPLING RATE ( MN 2 J F 8 
A(S) PRF ( KHZ ) 

*IT) NO, OF PULSES THANSMITTEO N 

'IS) SAMUiMOTh of a range cell ( NHZ ) 9W 

*;•)} RATIO Or FS/ttM 

*u«) total time of Flight i sec ) n/prf 
A im constant for the range walk crw 
v**3/(a*R) ( M ) 

A(U) RANGE RbSULUTlON RR 

C/Uaaw) ( M ) 

AU3) constant for the range AHPLITUOE function 
PI/RATIO 

*114) constant for the phase correction 
4*pi«crw/wl ( i/sec«*3 ) 

A(15) HAMHInG WEIGHTING ACROSS THE ANTENNA WINDOW 
VbS ■ I NO ■ 0 

A(i«>) nuhger of range Cells lookeo at 
A un range walk tes 4 0 NO • 1 

Alia) SANO-IOTH OF THE SYSTEH ( MHZ ) 

*IIR) range cell that target is in 
*(30) NO. OF SAhPlES per ooppler filter 
A(3t) range correction OPTION 

0 NO CQRNECTION 

1 WU SCHEME 
3 FRENCH 

S STEP TRANSFORH 

*133) RATIO OF THE SAMPLING RATE TO ThE LFH RATE 

*(35) linear RANGE WALK 

NO ■ 0 TES ■ MAI AMOUNT 
A(3b) MAIIMUM AMOUNT OF RANGE WALK 
A13T) TYPE OF INTERPOLUTION FOR FRENCH SCHEME 

3 IS 3 PT 

3 IS 3 PT 

4 IS 0 PT 

A(30) TEMPORT STORaTE UStO FOR THE QTNAMIC RANGE 


suhRoutine instr 


ki 



A2 


ORIQINAL PAGE ES 
OP POOR QUALITY 


common /IN3TM/AU) 

C 

c 

c 

C*U *S8IKN(fe, »TTI ») 

CALL 4SSI(iN(5, *LPI ') 

CALL A89mN(2, ‘SAMS, SMC '«e* ’9CR'»»2) 

C 

C 

c oi8« rue row tnc err waw oa'a 

c 

c 

MNire(6,s) 

S rOHMiT(' NAME or THE OI3K FILC '/J 

CALL A39IGN( t > 'OOMMY»,-l, *NtW*, ,2) 

c 

c 

c 

C ASK THfc UUtSriONS 

C 

C 

wW1TE(6iIS)) 

10 ru»MAT(' twansmitteo rwcooKNCy ( nnj )•/) 

A(l)aia90. 

0 WbA0(3,U} A(l) 

It row.MAr(M0,0) 

c 

•t 

V 

C WAVeLCNTN C/A(t) 

L 

C 

A(2}iS.b*2/A(l) 

C 

c 

wwiTEra»20) 

20 roWMA' (' VtuOcrTY OF T«e SPACF.CHAFf ( KM/SCC ) V) 

A(i)«f ,aS 

0 weAo(«>,n] A(j) 

c 
c 
e 

«wire(3«sn) 

So FOWMAtr WANOe TO THt TAROCT { KM ) »/J 

A la) aaso, 

0 RbA0(6,n) A(«) 

c 

c 

wmiTe(S«ao) 

40 FORMATC sampling RATe IN THE RANGE DIRECTION ( MHj ) »/) 

AC5JA22 

0 RbA0(6,U) A(9) 

C 

c 

Mwi re(6»5«9) 

So FORMAT!' pwr OF beams ( KHZ )•/} 

A lb) At .0 

0 WbAO(b,tt) A(b) 

c 

c 

WWlTEl 6 ,b«) 

bO FORm»t{' total no, of beams lateo down V) 

AlT)A«0Vb, 

0 RbAOCb.U) A(T) 

c 

C 


A3 


CRiGiNAL PAGft 

OF POOR QUALITY 


.vWlTF. (6,?a) 

JkS fUM^Arr BANOrtIUTM \if A HAN&S CtUL ( MMi »'/) 

0 atAU(h.lt) A(B) 

c 

c 

90 foawATC' MAM>nNO v<(eiGNT.N«» aCROHS ANTtNNA WINOOW */ 

1 ' Yt3 • I NO • 0 '/) 

A(lS)«iJ, 

atA0(6iii) A(is) 

c 

c 

MHITE (6»<)0) 

90 Fuai^AT(* NO, OF Range CeULS UOOKEO AT PER BEAN */) 

AU6)*I». 

0 H»AO((»,m A(19) 

C 

C 

waiTEtb.RS) 

99 FORMATC' 00 YOU '^ANT range wauk */ 

I * YE8 ■ 0 NO • I '/) 

AtlT)«3, 

0 RbA0(6.U) A(17) 

wRITE(G»9T) 

97 AO‘’‘’AT(' 9AN0WX0TN OF TwE SYSTEM ( MHZ ) */) 

0 REAOCb.m A(t9) 

c 

i«M! YE («• 49} 

98 FORMATC* cocation OF THE TAH6ET WHICH RANGE CELL V} 

A(19)«S, 

0 n&A3(6,n) A(i9) 

c 

c 

c 

WRITE (b>99) 

94 FUR"AT(* range correction OPTION »/ 

t ' none • 0 wU • I CANO ■ 2 STEP ■ S '/) 

A(2t)«l, 

A(0U*2i 

REAOCb.m A(£l) 

C 

17 (A (2t ) .NE.2.} GOTO 110 
wRtTtCbfIOB) 

105 FobmaTC' TYPE OF INTERPOCATION 5 OR « »/) 

»BAO(b,tn A(57) 

110 continue 

c 

c 

c 

c 

c 

c 

c 

c 

c Ratio of the satpcing rate to the range cell banowioth 

c 

c 

At9)»A{aJ/A(9) 

e 

c 

c 

c total time of the plight ( SECS ) 

c 


A4 


non •“ n nnnnn nnnnn nnnnn nnnnnnn 


ORIGINAL FAG : 

OF POOR QUALITY 


ACI!])aA(M 


constant F09 TKg RANGS MAI.K 
n/56C**3 


A(in«t(3)*A(j}/(2.«A(an 

A(l UaA(in«t.e«3 

RANOc R&SOUUTJON ASSUMING SMALL ANSLt OF INCIOSNCC 
A(l<f)a5.e*j/(<!.*A(ia}) 

constant F0« the amplitude weighted function 

PI»4,«ATAN(T.) 

AU3)>P1/A(<I) 


constant for the PmASE COfiKSCTlON 


Aa4)a4.*p{«Aan/A(a) 

A( 2 ^)* 0 . 

wwite(s« less) 

FUHMATC' linear range walk no ■ S yes ■ MAX AMOUNT */) 
RtAO(6,in A(2S) 


RETURN 

ENO 


A5 


ORIGINAL PAGE IS 
OF POOR QUALITY 


c 

c 

C St*'iM«.fOR 

C 

C 

c 

c 

c 

c 

c •**•*••*•• 

c 


c 

c • « • 

c 

c 

c 

c 

c 

c 

c 

leev 

i 

a 

5 

•t 

5 

6 

c 

c 

1010 

1 

g 

I 

a 

c 

c 

toao 

1 

g 

s 

a 

a 

C 

c 


c 

leso 

1 

2 

i 

4 


aUOROUM^E tNSTRA 


COMMON /lNaTH/*(l) 


W»tTE(3. (A(n,Ial,7) 

►'owMATC TRANSMIT reu ruenue^CY c mmi j ',tf>Ei3.8/ 

' KiAveLENCTM ( M) *,lPEt<>.S/ 

' VtUlCITY OF TMt 3P*CEC«*FT ( KM/SEC ) StPCtS.a/ 
' KAiMGt TO TmE target ( HM ) '»lPEIS,a/ 

• sampling Rate in Range ( MHZ ) *>tPEts.a/ 

' PMF IN azimuth ( KHZ ) 'ftPElM.S/ 

' NO, OF PUUSES LAYEO DOWN *»F7,l/) 


W»ITE(5, 1010) (A(n,l«»,14) 

Fu»mat(' sanorIOTh of a Range cell I mmz ) "itPElS.S/ 

' RATIO OF tne, sampling »aTE In RANGE TO BANOwIOTM OF A », 
' CELL *,IPE 13 ,«J/' TOTAL TIm£ OF LOOK ( SEC ) SlPElS.S/ 

» RANGE mALK constant { M/SEC »«2 J 'ilPElS.S/ 

' RANGE resolution ( M ) '.IPElS.a/ 

' CONSTANT ►OR amplitude FUNCTION 'ilPElS.S,* ( Pl/RATIO ) 
' CONSTANT FOR UUAOHATIC PHASE ’flPElS.Si ' (SEC«*-3) »//) 


«piTe(s,i» 20 ) (A(n,i>t 9 ,iR) 

FO»mat(» hamming aEIGmTING across antenna APERTURE *, 
» TES • I NO • 0 ST'S,!/ 

' NO. OF- RANOt CELLS TO PROCESS »,F5.l/ 

* range FALF yes • 0 NO • J ',F5,l/ 

* banulIDTh of the SYSTEM ( MHZ ) '.tPElS.S^ 

* LOCATION OF TARGET RANGE CELL NO, ',1PE1S.S//) 


OELLPm« 2,.A(3) «*a*An«)/(A(2)*A(4)) 

Of^aIbI/aIT) 

OFiJ»f)eLLFM/oF 

Al2.a)tOFd/«, 

A(22)«A(6)/0ELLFM 

WRITE(S» 1030) OELLFM, A (b) ,0 F , A( 20 ) , A( 22 ) 

FOPMAT(' change In frequency ACROSS NINDOW C RMZ ) ►, 

IPtlS.S,' sampling rate ( RmZ 5 *, 

IPtts.a/ 

' NO. OF KHZ PER FfT CELL ',IPE15,8, 

* NO. cells per OOPPLER filter '.IPEIS.S/ 

» RATIO OF SAMPLING RATE TO LFM RATE '.IPE15.8//) 


•/ 


A6 


o o o o o r»o Ci n 


OF POOR QUALITY 


wwiTECi, i(2t) 

K0«M*T(' COSKtCTION SCHEME NONE • 0 WU ■ I * » 
J » CANO • i SrEP ■ 3 '♦►s.l//) 


!»> (ACeu ,E0,S.) taHITE(S.IOttS) hUI) 

Po»M*T(' TYPE or XNTEMPOL*riON TnO OM F0U>» *»^8.l//) 


*NirE(8» iesu) A(S8) 

Will PUHMATt* LINEAR RANGE WALK «OOEU '/ 

I » NO • 0 ves • M** AMOUNT', FT.2/J 

c 

c 

HtrURN 

ENO 


A7 


ri o o o o •- o ooooonoooo o o o- o o r*o nooo ooonor>oo 


original page rs 

OF POOR QUALITY 




suaRourxNti tu, tempi, tcmpo* r anger# lmn) 


VlNTUAt Am)»4ANGEH(l)|Aa(l) 
VIStU&L TeHPI(i)»TEMPQ(n 
COMMUN /INSTR/A(1) 


TlMt«-AU0)/2, 

T»rbPat.E-3/A(&) 

Inau»-TIME 


FOPWATC B6G 
CUNwaA(U) 

ObUT AMat (<|) 

i>fet.TaaOfeLTARa«, 

CONSRaA(l8) 

CONSiaAlti) 

CUNBaA(ia) 


TlME.TSTtPjTFINAU 
IPE15,»#' STbP •, 


IPElS.a,* END 


MPEIS.8) 


NPULSeaain 

iMANGbaAUb) 

I«AtKBa«(l7] 

laPTSTaAClR) 

HUOPfa*(iI) 

IUtN)«aA(2S} 

IT (lUINM.NE.O) C0NLB-a,aA(25)/{A(T)aTSTEP) 


iwae 

continue 


compute the AMOUNT OP range NALK 

TSO«TlMEaTXMe 
PANMtKaid, 

IP (1«ialAH,EU,1) goto 9 
MAuaaCONuaTSQ 
PANMi.Aa>A(.K/CONSR 
5 CONTINUE 


A8 


rtnrtorto o nn« o> oooo«n n> or* oonnrt o o o o r» o 


POOR QUAUTY 


c 


IF (lL(M<tN£« 9 } ttANMLi^aHANMUK * C 0 nU*T 1 ME 


fI^o the phase 

PHASE«rSQ*CONS 
tnllnl ♦ t 
Al(lri)*PHASe 
AU(lw) arfANwCK 

TiMtiTiME ♦ Tsrep 
IF r'.W.LT.NPulSE) Coro 1 


FlNO TMb MIN AMO max RANGE MAES 

XMAXAAQ(I) 

XH(NaAQd) 

UMAXal 

LMlNal 

00 4 la2, NPUI.se 

if(aq(i}.lt.xmax) goto 2 

XHAXaAQ(t) 

LMAxal 

IF(AGU).GT.XMIN) goto 4 
XMlNaAO(I) 

UMlNal 

CUNTlNUE 


NOPMALUfc THE map TO THE MIN/MaX OF THE RANGE WALK 

NPANGEatnANGE 

LMNal 

XFCxmax.EO.O,) goto 9 

XXMAXalFIX(XMAX) 

IF (XXMAX.Ng.XMAX) XXMAXaXXMAX * I. 

NR ANGEbnHANGE * IFIX(XXMAX) 

IF (FM(N.ut«!9.) goto S 
XSa-XMlN 
XXMINaiglx(XS) 

I^IXXMIN.NE.XS) XXMXNaXXMZN A 1, 
NRANGbaNSANGE * XXMXN 
tMNa-xXMXN ♦ 1 . 

IF (I«uOMT.EO,l) GOTO 9 
NH ANQgBNRANGE ♦ 2 
UMNaUMN • X 
CONTINUE 


wPITE(E»i>(dSe) LMIN. XHIn.LMAX.XMAX.lMN. XRANQE.NRANGE 
PORNATI' range information »/» MlN X(»« IX. IPEXS.S/ 

X ' max * . Xb. tx, XPEIS. 8 /' offset ', 16 .' 
e » NO. wanted ', 16 ,' NO PROCesStO ', 16 //) 


COMPUTE the signal for EACm SEaM POSITION 


A9 


ORIGINAL PAGE IS 
OF POOR QUALITY 


C <S « FUNCriO'J OK »AN6ti 


c 

L»t 

ti9 COMriNue 

c 

iwtl 
I maUMN 
PHASiiaAKl.) 

CSaa'SlPHASE? 

SSs31N(PH*Se) 
hAl'iWk.KaAQCU) 
li CONTINUE 

C 
C 

TeMpniw)t0, 

TtMpailW)a0, 

IP (Iv^.OT.NMAnSE) aOTO 20 

c 

*lwtuw . IOFF3T 
XXtWaXlw » MANMLrt 
CAUU SINX (VI • XXlWf CONSX} 
xxiwaxxlw ♦ OEtTAP 
CAUL SINX(v2,XXlta,CONSX) 

XXlvaaXXlM - 0EUT2 

CALL 31NX(V3,XXIm,C0N8X) 

C 

c 

c 

OOavi * 0,tt2ba(V2 ♦ V3) 
lli^) »OOaC9 
Tt‘*PO(XW)aOO*SS 
C 

TMalw * t 
llwaXiM a 1 
SOTO X2 
C 
C 

20 CONTINUE 

C 

c 

c 

C 

c 

MRITgCS) (Tempi (J) ,J al,NPANGE) , (TEMPO ( J) ,JaX , nRANGE) 
C 


L»L ♦ I 

XFCL.LE.NPulSE) goto 10 

c 

c 

100 CONTINUE 

C 

c 

C REARRANGE THE RANGE 
C 

L*NPUL3E/2 
UU 110 J81,NPUL9E 
L»L ♦ 1 

IP (L.f.T.NPULSt) L«1 
TtMPPaaQ(j) 

RANGER (L) aTEMPP 
110 CONTINUE 

C 

e 

c 

A (30) aNKANGE 


AlO 



ri n o o 


OmGINAL page 

poor, 


wOlTE(o>/«'<9) LXLK>'VH*N6e(NPUL3e 
7091 P 0 **«AT(//» NO, P» 0 C SX 6 »* NMANG »,X 6 ,» NPUC 5 E 

C 

PC TURN 
ENU 


suenouTiNE siNx(v,x,cuNsrj 

v»t. 

oo*coNsr*i 

0«3S3*AHS(00) 

lP(0A83S,GT,t.E*E) V*SXN(00)/00 
RETURN 
ENO 


All 
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c 

c 

C SARliJ.^OH 

c 

c 

c 


G 

c 

c 

c 

c 

e 


c 

c 

10 


It 

c 

c 

c 

c 

li 

z 

c 


suannurtNE outp(ax>ao»nanger) 


vtaru*t ai()0),4O(S0),r«n(»er(|) 

at*L i*(a 0 ) 

COH.'iQN /INSTR/Ad) 


MMITE(6»t0) 

FORMAT c' Plot out the time ReaneNSE TP.i ■ 0 no • 

IHAXtl, 

REAO(b.lt) KHAX 
FORHAT(M0,0) 


IF CVHtx.NEtO.) GOTO 10000 


00 la t«t#a 0 

IA(l) ■ 100000 , 


RERINO 0 


c 

c 

c 

c 

c 


NP0tSb«4(T) 

NRA'M&E«*(30} 


FliMU The haxIHUH KEER THE aZIHuTh anq RANGE BIN 


c 

c 

c 

c 


XHAXB0, 

00 i00 jBifNPULaE 

REAO(0) (Aim ,la||NRANGE} > (AQ(I) , lat .NRANUE) 


UO 50 l«l,NRANGE 


Ai(n«Ai(i)*«a ♦ AQci}*«a 
IR(Aim,t.T,XMAX) GOTO 50 
XMAXIAI (I) 

KJaJ 

Alai 


C 

50 CONTINUE 

C 

C 

100 CONTINUE 
C 

e 

c 

C SCALE ALL The POINTS TO MAX PUT IN LOG SCALE 
C 

c 


1 V) 
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c 

c 


t 

i 

e 

e 

c 

c 

c 

c 

c 

c 

c 


c 

1S0 


tvea 


c 

^e^ 

c 

c 

c 

c 

10000 

c 

c 


NH tNa^ltANUC 

IF (N 4 «N,GT,^ 0 ) NMAM 80 


4lHire (5|900) KI«KJ»4Ht<, (|,|al,NH*N) 

P0H'**T(//* HANtte ••la#*’ AZIMUTH '*I6» 

' OF MAIIMUH value ( 0 Oa a£FeH£NCfi ) 'itPElE.d// 
9l,00(0X,l2,2l),' WALK '/) 


PCMINU 2 


JJ*NPULSe/2 
00 200 Ja|,NPULSe 

ff£A0(2) (AI (1) »Xal,NMANGft)» (Aa(Z)»Ial»Nf<AN6e) 
00 190 IF1,NKAN 
IA(I)al000O0. 

U0a(AI(l}aa2 t AO ( I ) aa2) /XHAX 

IF (OO.GT.e,) IAma.i0,*ALO610(OO) 

coNriMue 

jjajj ♦ 1 

IFUJ.r.T.NPULSfc) JJtl 

j» (iA(n,iai,a 0 ),HANGea(jj) 

FU«mat(IS, 20(1X,FS.2] , tX«F0.9) 

Continue 


CONTINUE 


RETURN 

ENO 
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c 

e 

c sAMuairofi 
c 
e 
c 

c«***«* 

C WRITES our > ^lUE IN THE ruLUOWlwC ^gflMAT 
c 

C ReCORO «1 

C aRRAT a initial CONPISURaTION ( real ELCMeNTS ) 

c 

C RECORDS • e TO • range CELLS » 1 

c 

C I COMPONENT AHHAlf OF FFT *n0 Q COMPONENT ARRAY OF FFT 

C ( FOR THE STEP THANSFOHM IT IS NUN FFT DATA ) 

C EACM array is of length 0 OF SEAM POSITIONS REAL AHRATS 

C IF A flF reams TS not a power of 2 THEN THE LENGTH 

C OF THE ARRAYS IS nEiT HIGHEST POWER OF i 

C 

c 

C RECORD p RANGE CELLS ♦ 2 
C 

C Range walk for each seam position 

e 

c 

C**A****AA*A«A«AAA 

c 

SUBNOUTINE FFTRIAIi A0»XXX,XX0,RANGER»LMN) 

c 

Ct*«*«AAA*AA*A**AA 

c 

6 

VIRTUAL XXKSOrAO) »YXU(SO,40) 

VIRTUAL RAN&ER(t),AI(l)»AaU) 
real Xl(90)fXU(S0) 

COMMON/INSTR/A(50) 

C 

C 

C bring in each range row to make up 
C the azimuth array 

c 

c 

NPULSEaS (F| 

IwALKaAO/) 

IWINUaACtS) 

nRANGE«A(I6) 

ISCHHtA(2l) 

iNTPiA(2n 

IRANGEaA(30) 

RaTT0*A(28) 

XHAXo AA (2b) 

*NPAFtOAT(NPULSE/2j/K»TIO 
jnPianP * 0,S 

0 WPITEft. T9RRRJ JnP,NPULSE»XNP, ratio 

TB999 FURHAT(» JNP ',I6»* NPULSE ',Ib»' XNP »,IP£13,S,* RATIO *»lPeiS«a) 

c • 

c 

c 

C0NSTaA(2)**2*A(4)/(8,*A(3}*«2*A(12)) 

CONST»CON8T*t .£• J 
FSrtPaA(b)«l,E*3/FLOATtA(yj) 

0 WHITE (biSbbe) CONST, F 9 TEP 

5S68 FORMATC COST *,lPEl*i.d,» FSTEP *,IPEIS,8//) 
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c 

CaNL«-2.*«(£S)*A(2a}/rt.U«T(NPULS£) 

c 

c 

c 

c 

c wntre uuT rne configuration 
c 

URITEU) A 

0 WMirt (o.acRi/) NPUL XRANUE» A(|) 

9080 FUWNat(I6,U|X6>1x XRblStS) 

C 

C PUT NE range mALK amount in to file for the step 

c 

IF (ISCHM.E'.f. ' .RXTE(i) (MAN 6 ER(X)»X*ltNPUl.SE) 

C 

C 

00 100 i*l>|RAN 6 E 
C 

REWIND S 

e 

DO S 0 Jal.NPULSE 
C 

READ (a ) (XI(N) ,Ha|,lRANGE) , ( XQ (K ) , K • t , XRANGE} 

AX(J}«XI(n 

AI 3 (J)iXQ(X) 


C 

90 

C 

C 

c 

c 


c 

c 


c 

32 

c 

c 

c 

39 

C 

c 

c 


CONTINUE 

WINDOW THE DATA 

IFdwINO.NE.n GOTO 93 

XKKW0.0006TR9 

!(NN«-NPULSE/2 

00 92 jalfNPUUSE 

XXaXNN*XKK 

IF (AttS(XX) .61.0.0000 VA(.a(SlN(XX)/XX)**2 
Al(J]8AI(J]*VAt 
A0(J)iA0(J)*VAL 
XNNbxNN * 1 . 

CONTINUE 


CONTINUE 

COMPUTE THE fft nr ' -nce celc 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


IF (ISCmm.nE.SJ FFTVINPULSE, Al,AO) 


WRITE OUT THE AHRATS TO TmE FlUE 

WRITE n ) (Ai (J) , J«t .NMULSE) , (AQCJI.jal.NPULSE) 
WRITE(b.a 090 ) NPULStfl 
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t I-) 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

U9 

c 


e 

c 


CUNTINIJt 


PIMO THb UANGE walk f'ON EACH ANGLE TAG ZT 


IP ;iSCHM.kQ,0) GOTO 1000(9 
IP (ISCWN.tO, J1 GOTO 10000 


0 EW 1 NO e 
NPsNPUtSfc/O * 1 
NWA!iai 40 

NPHlNMULSE/NPASS 
NOVEHiNPUtaE - NPPtNPASS 
ISTA»T*1 

IP (NPP.LT.l) GOTO 150 
iMPiptNPucse ♦ 1 
NPP?«NPIP ♦ 1 
ISTS«~I.MN ♦ i 
UNPalSTP ♦ NPANGE - 1 
iHNGPalPANGt • NPANGE * I 


00 120 J»1#NPP 

PtwiNO I 
KtAU(l) 

00 106 Kat^IPANGE 
I8»:STA0T 

PEAOCll (AI(L),L«l*NPULSe},CA0(Ll,Lal»NPUL3E) 


DO 109 LLal>NPA33 
C 

XXI(N,LL)«AI(X3) 
XXa(K,LL) aA0(I3) 
I 3 'tS ♦ 1 
C 

109 continue 

c 

106 CONTINUE 

c 

c 

C COPPECT THE APRAV IN RANGE 
C 

c 

00 110 K«1,NPA33 

c 

JJalSTART 


103 


104 

c 


wA|.KRa 0 , 

IP (IWALK ,E0,n GOTO 103 
IFCJJ.GT.NP) JJaNPP2 - JJ 

VALUfcaCPLOAT (JJ - 1) apSTEPl »*2aC0N3T 
‘«ALKHavALUE 

IFULINW.EG.O} GOTO 104 
IFIISTAhT.GT.NP) jja ISTAPT - NPUL3E 
VALa«CONL*PLaAT (JJ - IJ 
WALHPaMALKR ♦ VALE 
CONTINUE 

IPdSCHP.bQ.E) GOTO 103 
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t'J 6 b 


107 

0 

6060 

lOS 


c 

c 

le 4 

C 

9019 

C 

C 

tlO 

c 

100 

c 

c 

t jo 

c 


c 

c 

I J9 

c 

c 


c 

c 

c 

c 


c 

161 

c 

c 

c 

c 

c 

las 


i-L'»*d3iw6U<«> ♦ 0,6 ORIGINAL rACu fS 

POOR QUALITY 

le*ItNO * IwU 
If^dS.GT.O) GOTO 1066 
18*1 

lEiNMANGfi 
UOTO lOT 

lF(U,ce.IH*NGf>) UOTO 10? 
lb*19«NG£ 

ISalKNGR 

w9tTb(0) (XXI(t.U#K},l.l.aIS.lE), (XKQ(U»K),l.L«X8,Xe) 
<<Nirb(b,b06O) l8r*MT,w«LK9,UU,l8,lE 

?u9mat (lb, lx, ii>&i8.o, j(ix, isn 

GOTO 109 
CONTINUE 

KCLiK 

CALL C*NO(XXl, XX(9,N«LH9,NHiNGE,Kl.Li I^TP, ISTN, INANGE) 
t»IUTt (S) (XXX (LLiK) ,Uttl,NH*NGE} , (XXQ(LL,K) ,LL*1 ,NNAN6E) 
coNTXNue 

wniTb(6,9019} 1START,JJ, VALUE, VAL2#MALX9, X8,XE 
?ORNAT(0(tX,Xtt) , Snx, 1PE16.S) ,211X,X9)) 

XSTARTaXSTART * 1 

CONTXNUE 

CONTXNUE 


XT (nOvEN.lC.O) goto 139 
NPP«1 

npasSanOVER 
NOVE9*0 
GOTO 119 


CONTINUE 


REwXNO 1 
9EA0(1) 

DO 100 Kai^NRANGE 
PEnXNO i 

00 131 JairNPULSE 

»E«0(a) (XX (t) ,L*1 ,NHA''6E) , (X0(U N9ANGE) 
AIU'axXlK) 

AQ(J)axO(K) 

continue 

m9XTE(1} (AX (L},IpIiNPULSE), (AQ(L},L«1,nPULSE) 


CONTXNUE 


1U0OJ CONTINUE 
C 


rt 0 «t 

c 


WRITE (6,eaai j npulse 
fosnatc la3t one 


call CtOSE(l) 

U 

c 


PtTUWN 

ENO 
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c 

C SAKGaA.I'Qf} 

C 

C 

C 

SUSiRQuriNe CAN 0 (A 1 « INrPi tsrH,IRAN 6 it} 

c 

virtual Ai(*>a»a(9}»Aa(90,40) 

c 

c 

irauaffu 

OtCTA»»« - FUOAT(IRW) 

U»l 

Jaisrs ♦ IRW 
If {INTH.eg.tt) JtJ -I 
If (If.rp.£g,3) JaJ •! 

IMOiiLTA,LT,0,) JaJ - 1 
If (OtlTA.LT.a,] OEUTAat, * DELTA 
IfCJ.LT.lJ J«l 
IRAN(>sIRANi)b - N 

0 MMire (3,5’i'ill) K,J>ISTR>MW 

0 MHlTECSibtiaS} (JKf AI(JK»K;iJKat«lRAN6e) 

5M50 FO»M*r (3(I'>, t *) . IPElS.aj 

SUSS Format ( 13, IX, IPbiS.a, IX, 13, IX, IftilS.At >X, Ii« IXi IPEIS.d} 

tf ClNTP.gQ.a} CiOTO 1E90 
If (XNrp.gy.3) GOTO 300 
1RAMG«IPAMG • 2 
If CJ.GT. (R ang) JalRANG 
0tLTA3»0f LTA aOELTA 
0bLTA3a0gLTA2*0ELTA 
C 

AMl»-oELTA*(OELTAa - 3,*0tLTA ♦ a,)/6. 

AlDa(0f.LTA3 > a.aObLTAa - OElTA * 2tJ/2, 
AI»-13tLTAa(OELTAE • OELTA -a,J/2, 

AaaOELTAa (ULLTAa - l,}/b, 

C 

c 


I 

c 

c 

c 


c 

c 

c 

c 


c 

iian 

c 


tul 

c 

c 


CONTINUE 

VALUETaAMI •*! (J,K) ♦ A0 »AI(JaJ,KJ ♦ AlaAl(J*S,«) 
VALUEntAM|*AO(J,K} * A0*«U(J*1,K] * A I • A3 (JtS , >< ) 

AI(L»K}avALUhl 

AQ(L,i<}avALUCQ 

JaJ ♦ I 
L»L * I 

If (L.Ue.N} GOTO I 


all done 


GOTO 10000 

continue 

OELTAaai , - OELTA 
If (J.GT.IRANG) JalRANG 
CONTINUE 

VALUEIaOELTASaAl (J,K} ♦ CELT A a A X( J« I , R } 
VALUEOaUELTAaaAU(J,X) ♦ DELTA* AO (Ja 1 , K) 


a A3»Al (J*3,A) 
a Aa«A0(3a3,K) 
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at(U.x)»v'Ul 6 l 
1 'VACUtO 
C 

JaJ ♦ I 
L«U ♦ I 

GOTO 10t 
C 

GOTO 10000 
300 CONTlNUe 
C 
C 


ORIGINAL I-/':-' 

OF poors QUAun' 


iHAMr.alHANG - 1 
IF(J.bT,lHAN 6 ) JalRANS 

e 

AMl»OeLTA» (OELTA - l,)/ 2 , 

40*1 , • 0 tLT 4**2 
At ■ 0 EtT 4 *( 0 bUTA * t.)/ 2 , 

C 

101 CONTIMue 
C 

v 4 i.uei* 4 >n. 4 Hj,K) ♦ ao*ai(j*i,io ♦ ai*ak«u 2 ,k) 
VAUUEQaAMl*A(3CJ|K} * A0 « AQ ( J * 1 , H ) ♦ A1*AU(J*2|K) 

c 

AKUfOaVAuUei 
AOtUfK) aVAtUbO 
C 

JaJ ♦ I 
Uau ♦ I 

IPCL.UE.N) GOTO 301 
C 
C 

C ACL DONE 

C 

C 

10000 continue 

0 WMITK (S, (JK| AX (JK,K} , JKalfN) 

HETUON 

ENO 


Ca 

C» FTTV.FOH 

Ca 

Ca 

Ca fFT ALG0T«IT^^^ FUR VIRTUAL ARRAYS 

Ca 

C» 

sushoutinc; ^ftv(num«r,xt,tt) 
virtual XT(1) ,YT(1) 
c 

c NUMUR IS The no, 0 ^ POINTS OF liOUO OaTA 
C IF NOT 4 POnER of two PAO TO NEXT MIOMfeST 
C WITH ZtROS 

C 

c 

NaPOWal 
I TfeHPaNUM9H 

I ITEMPalTfeHP/a 

IF cneHP.ut.n goto 10 

N(JPOW«N2POW ♦ 1 

goto 1 

10 continue 
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c 

(N.Gg.NUM|M«} GOTO 19 

M«M*2 

^l^POl«■N^t*Ori ♦ 1 
19 CONTINUE 

(NUMGM.ce.NJ GOTO 20 

c 

NIbNUMHR ♦ 1 
00 13 |.^'■N1>N 

c 

VT(UM)«0, 

c 

19 CONTINUE 

c 

c 

20 CONTINUE 
C 

C 

M«N2P0W 

00 b00 l.0al»N 

tMXa2**(M»L0) 

L 1 X« 2 «LN| 

SOU «b. 209 1 a3/M.0AT(t. IX) 

on (»||)0 Uf^alfUNX 

* 4 G*(UM>U* 5 CU 

C*COS(AHG) 

SaSIN(AR 6 ) 

00 000 UI>UIX|N,U1X 
Jl«U-UIXaUM 
J2aJl«UNX 
TiaxT(Jl)-XT(J2) 
T2avTtjn-VT(J2) 
xr(jnaxrtji)axr(j 2 } 
VT(jnaYT(jn*rT(j2) 
XT(J23aCaTla8aT2 
rt (J2)acaT2-S*Tl 
000 CCNTlNUb 
Nv2aN/2 
NM}aN>t 
J«1 


OU 035 laj 

IK(I.GE.J) GO TO 631 

TiaxrCJJ 

r2«YT(JJ 

xru)sxr(X) 

rr(j)«YT(i) 

xTa)«n 

YTinafa 

6^1 43 NV 2 

620 CONTINUE 

If (K.GE.J) GO TO 633 
JaJ«i( 

»a«/2 
GO TO 620 
695 Jaj*K 

NUMbRa^ 

RtTUMN 

EnO 
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MNUOV,(<0H 
V1I4TUAC VEXStON 
ttlNOOWlMG FUNCTION 


ARGUNENTS 


AHRAy IS THE ARRAT TO ttE xll.iUORbO FIRST UOCATION 

ISTART is TnE first Ei.ErEnT of the ARHAV to »E WlNOUWEO 
NUHSR IS THt NUHdtR OF ELtHtNTS TO UE WINOOWEO STARTING AT ARUVE 
lOPT IS TYPE OF RlNOOR 


Id NO wihdorinu (Rect function ) 

I HAMMING WINDOW 

e sartcett window (triangle) 

3 HANNING WINDOW 

4 DLACmman window 


SUDROUTINt WlNODVlAHRAY, IST aRT|NUHRR« lOPT) 


VIRTUAL AHHAV(t) 


IF (lOPT.EOf Dl GOTO l<£idigO 
XNlNUMHH - 1 
ISTaT«ISTART 
IF (I0Pr.EQ.^) GOTO 2ea 

PIia.tATANCl.) 

PI2I2. *P1 
PHASEuPia/XN 
IF (10PT.EQ.4) GOTO lOfl 

Cl'dl.S 

CgaU.S 

IF tlUPT.EO.i) GOTO 10 

Cl*e.S4 

C2B0.46 

Id CONTINUE 


DO 20 I«i,NUMOR 
«N«I - I 

wuaCl • C2«C0S(PHASE*XN) 
ARRAY (ISTAnawO* ARRAY (ISTAT) 
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c 

c 

c 

c 

c 

c 



c 

aia 

c 

c 

c 


c 

aai 9 

c 

c 

c 

liaaeia 

c 

c 


I 8 T*Ttisr 4 T ♦ I 
cuNTiNoe 
GOTO 
CUNTINUti 
PH«Sgaat>H* 3 e« 2 , 

00 laa x«i|NUfi 8 R 

»\«I - t 

wolu.aa - a,sij*cus(Pn«Sfc*xN) * 0.oe*co3(PHAs'ia*iN) 
*I«'UV CISTftT) awl}*AHHAY (tSTAT) 

XSTAfalSTAT ♦ X 

contr XNUK 

GOTO xeoeo 

NUMhR is Tne NUHt)b»< OF EttMeNTS XN THb MXNOOM 3 TAt 4 TXN 6 
FWOM THt FXHST AOU*<bSS OF AHRAV 


.CONTlNUe 

lUP*NUM0R/2 

00 ail« XalilUP 

WDst-LOATC^Ad - 1))/XN 
AMiiA f USTAT)aMO*ARNAV(XSTAT) 
ISTAT»X3TAT ♦ X 

CONTINUE 

XUPaXUP ♦ I 

00 23 a XaIUP,NUF' 9 P 

woaa. - floatc3*(x - nj/xN 

ARAAVIX3TAT; aMO*ANNAT (XSTAT) 
ISTATaxsTAT ♦ X 

CONTINUE 


GOTO 10000 

continue 

RETURN 

ENO 
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PART II - FFT CONVOLVER AZIMUTH CORRELATOR 




P*<0QKi*M TO *>R0CtS9 THC 5A'« DATA FROM SANG PROGRAM 


VIRTUAL AST (AURS) » A3C(4ia<>b) 
VIRTUAL P 1(01^46) 

common /INSTR/AOO) 

COMMON /TfcMP/A7.a(8) 


Call '“stn 
call GgNCriiPQ) 


00 I t*l«4 

CALL OOPPLRCPI.FO,!) 

call SlGNALU9If AS0,FI,T(3,I) 

CALL 0UTP0(TQ,F1,AS{« ASU) 

CONTlNUfc 


ENU 
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Mart 


Of PC OR 


QUAUT't 


tNSTsucTioN ser fan the s*i» 
c IS rne speeo op wight • j. e * s ( m/scc ) 

PI • tt.aATAHCl,} 


A(i) THANSHiTTeo PneuueNcy ( hhz ) pt 

*(«) i«*VtW6HGTM OF FT ( H J C/FT W 
AIJ) VtLQCTtT OF SPACeCTAFT ( KM/86C ) V 
AC4) SANCfe TU TAHGfeT ( km ) R 
A(S) SAMPUNU NATE ( HHZ ) FS 
A(6) PMF t KHZ J 

A(T) NO. Of PULSES TkanSmITTEO N 

A(S) SANOMtOTH Of A HANGE CELL C 'ThZ ) BW 

A(9) RATIO Of FS/Bw 

At|0) total time of flight I SEC ) N/P«F 
A(U) COMSTa.mT for the range walk CKH 

( H ) 

AUE) range Resolution rr 

C/(E«Bh) ( H ) 

AUJ) constant for THE RANGE AMPLITUDE FUNCTION 
PI/RaTIO 

AUQ] constant for The phase correction 
4*PI*CR«/HL ( l/8EC»*a ) 

A(ts} hamming weighting across the antenna wINOOM 
TES ■ I NO • 0 

Alts} number of RANGE CELLS LOOKED AT 
A(tT) range MAlk TES * 0 NO * I 
Alts) BANOWIOTh of the STSTEM ( MHZ ) 

AllR) location OF The target 
A 130) NO. OF CELLS FOR EACH OOPPLER FILTER 

A(3t) range correction option 

0 NONE 
t WU 

a canoxans 

AIE3) ratio of the sampling hATE TO THE LFM RATE 
A12S) linear FM hALK AOuEO no ■ 0 TES ■ MAX AMOUNT 

A(3T) TTPE wXNOOR ACROSS FILTER 

At39) STARTING LOCATION IN FREQ OF THE FILTER 

A130) STOPPING LOCATION OF THE MATCHEO FILTER 


SUBROUTINE INSTR 


COMMON /1NSTR/A130) 


CALL ASSIONtS# 'TTI ') 
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c 

c 


pa:..; : 

OF POOR QUALi; 


J 

V 

a 


MRireCb, I) 

1 FCRM4TC NAME 0^ ClNePRINfEM OUTPUT PILE */ 
i ' tPt UKiSSNAME MTOi name */) 

C Caul *SSUiN(5» »OUMMf »,-i) 

CAUL A98l(iN(5,»l.P«*J 

c 

2 FUbMATt* NAMb OP PILE POM OUTPUT ( OAI OR MTBiNAME ) '/) 

CALL A S 8 ION (A I 'OUTP.SCRSU* 'SCR*»# 2 ) 

CAtU A'(81SN(3» 'OU7P2.SCHM9, *SCR*m2) 

C 

C 

c NAME OP 3 I 3 K PILE TO OE PRUCbSSEO 

C 

C 

t>HirE(6»S) 

5 PORmaTC name OP 0I8H PILE TO Bt PROCtSSEO */) 

CALL ASSIGN! J I 'OUMMY*,-!) 

C 

c HEau In THb configuration 
c 

REAum A 

C 

c 

e 

c 

e 

WRITEIS.IRUH) (ACD.ltt.n 

iaJI0 FORMAT!' TRANSMITTED PREGuEnCY ( MMJ ) » iPEiS.S/ 

I ' «AVEl.£NGTM ( M) 'flPEIS.a/ 

6 ' VELOCITY OP THt SPACECRAFT ( KH/iEC ) *»tPEtS, 8 / 

5 * range to The target { AM ) »,iP?j5,e/ 

4 ' sampling rate IN Range ( mh 2 ) ',lPElS.a/ 

b ' PRP OF Tme aIIMuTh t RM 2 ) 'jlPtlb.a/ 

6 ' NO, OF PULSES LAYEO OObN *,FT,l/) 

c 

c 

MRlTEtS, 18)10) (Am,i»e,ia) 

1010 FQRMAT(» MANOWIDTM OP A RANGE CELL ( MHZ ) '.IPElS.a/ 

1 ' Ratio of the sampling rate in range to RANOwIOTM op a 

2 ' CELL %|Ptlb.a/' total TI"E op lour ( SEC ) '»tPE|S,a/ 
i ' Range naLR constant ( M/SEC»«a ) '»IPE15,8/ 

4 ' RANGE HE80LUU0N ( M ) IPElS.a/ 

b ' constant POR amplitude function 'dPEIS.a,' ( PI/RATIO ) •/ 

b ' CONSTANT FOR OUAORATIC PHASE *,lPE|b,S, ' C8EC**-2) »//) 

c 

c 

MHiTEisaoaa) (A(i),i«ib,ia} 

1020 FORMAT!' MAMMINO pEIUMTINO ACROSS ANTENNA APERTURE 
I ' YES • I NO 1 0 ',Pb,i/ 

e ‘ NO. OP range cells to process '»ps,i/ 

3 ' range hALR yes « 0 NO ■ I 'fFS.l/ 

4 ' BaNOWIOTm op the system ! MHZ ) '.IPEIS.a/J 

c 

c 

OeLLPMa2.*A!3}«*2*A!t0)/!A!2 ) aAC4]) 

OL»A!b)/AII) 

A(20}iOELLPM/OL 

Ai20)aA(20)/4, 

C 

C A90VE assumes hE have FOUR OOPPLER FILTERS 

C 

WRITE!?, J0l0) OELLPM, 416) ,OL,A!20J 
10)0 FORMAT!' change IN PWEOuEnCY ACROSS wINOOw ( RHZ ) *, 
t IPtlb.a,' sampling rate ! khz ) 
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ooouo 


i 

4 

c 

c 

1040 

I 

c 

c 

10S0 

c 

U53 


1060 

1 

2 

0 

1061 

C 

1060 

1 

2 


OHIGINAL PAGE S 3 


ipe.l3,»/ 

' NO. OF *i«z Pe« FFT S4MPUE '.iPsn.a, 

' NO, OF ctcts PtR OOPPtES FltTEH »,lP61S.a//) 


HRI rE(9, Il«40) 4(21) 

F0MM4T(' HANGe COMHECrlON F«CT0K NONE *0 WU • 1 *» 
• CSN40I4NJ • 2 


wRITe(5, 1030) 4(2S) 

FUNmaT(//' linear f» i»4lK AOOEO no • 0 VE8 t NAX ANOgNT *|F6.2//) 
NMITEta, 1133) 4(22) 

fobm4T(//* location of the TARTGET »,F3,1) 


UH1TE(6.I060) 

FORMAT!' TYPE OF wInOOR ACROSS MATCHEO FILTER •( 
• NONE ■ 0 hamming ■ I aAMTLSTT ■ 2 HANNING ■ J 
'RLACKMAN ■ 4 '//) 

*(2T)«1, 

REAO(6.106I) A(2T) 

FURMAT(F10,0) 

«HITE(S. 1062) A(2T) 

FORMAT!' TYPE OF WiNOOW ACROSS maTChEU FILTER '/ 
' NONE ■ 0 hamming ■ 1 0ANTLETT ■ 2 HANNING • 3 
'BLACAMAN • 4 ',F5,2//) 

RETURN 

tNO 
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ORIGINAL PAGE IS 
OF POOR QUALITY 




SUeROOTlMK CeN(ASI|«SQ) 


VICTUAL «9I(1) **5Q(t) 
COMMON /|NSrR/4(l) 


GENe»Are our tmi> ot^i-xENce signal 
rOH MtTCMEO FILTEM 4 N 6 In the ANKLE OOMAXN 
LlNtAM l»AVE ^QRM 


NPULSE*«m 

TSTEM»J,t-3/*(fc) 

TlMt«-ACte)/2« 

TFInal«-TIhE 

C0N9«*(14) 


continue 

IFU.CT.NMUtSE) GOTO 100 


T80»TImC*TIME 

PHA3b«TSU«C0NS 

XI*C0S(PMA3E) 

ASI(t)«XX 

ASOCnaXX 

TlMgsTIMt ♦ T8TEP 
l»l ♦ I 
GOTO I 


00 


CONTINUE 


WPITE out TME TOTAL WAVE TO A TEmPOPT FXlS 

WHITE (2) (Asnj)*J>WNPuL3E], (AS0(J],J«t,NPULSE) 
C 
C 

PETUPN 

ENO 
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ORIGINAL PAGE 13 
OF POOR QUALIT/ 




aijhaoutine oupp(.n(^l,Fa, tr*6) 


VISTUAL FI(n»PCl(l) 
COMMON /INST»</l(t} 


NPULSEaA(7) 

NOMdaNPUtSt/a 

ti«iNow«A(er) 


PCMINO i 

Pt*0(d) (FI(I) .lai.NPuCSe), (FQ(I),Iat,NPUW9E) 


NUUt OUT THE TIME nISTOPY OF lYO XNTEpEST 


Ktal 

If (ITAG.EQ.n GOTO 20 


<l)a(IT»G • l)*NUM8 
»Ua«U ♦ 1 
00 10 Kai,KU 
C 

FX(H)a0, 

F 0 (»)aa, 

C 

10 CONTINUE 

C 

C 

IF (XTA6,Ea.«) GOTO S0 
C 

20 CONTINUE 

C 

C 

KualTAGaNUMO a 1 
00 28 Kat(u#NPULSC 
C 

FI(K)aO, 

FQ(K)ae, 

c 

28 continue 

C 

C 

S(9 CONTINUE 
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OF POOR 




c 


c 

c 

c 

c 


PUT 


C 

c 

c 

C FIM) 

C 

C 

oast 

c 


ipciMiNjw.eo.in GOTO oast 


HtMMtNS '/jbtGHTING ACROSS r^lLTeR 
|ST*«T»«L 

CALL v<l*iOOV(r' I,ISTAHT«NuMa,XrilNUW) 
t$rAHr»Ki. 

CALL W1N00V(FU> (STAHTiMUMItflMINOM) 


TH6 SPECTRUM OF THE REFERENCE 


continue 

NUM«R*NPULSE 


CALL FFTV(NUM8«»FI,F0) 

C 

loeoo continue 

c 

KKMMU • KL * t 
WHIT£(6»/S0} KKK»KU«KL 
T»9 FORMAT(A(lX,Ib),* DEL UP LOm *) 

CC 

C STORE awaT TmE starting ANU STUPPPING LOCATIONS 

c 


Ataoj*! 

A(iO}ANUM8R 

c 

c 

c 

c 

c 


RETURN 

ENU 
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ORIGINAL PACr, 

OF POOR QUAtn a 


c 

c 

C MAT4,PUR 

c 

c 

c 

c********** 

c 


c 

c 

c 

c 


c 

c 

c 


c 


c 

c 

c 

c 

c 


SUOROUTXNE SXGN*L(*S{,*Sa.M,P0iITA6) 


VIRTUAL ASI(1)|ASQ(1) 
VIRTUAL Flm 
COi-imon /InSTR/AU) 
CUMHUN / rEMP/TtM(*(l) 


NPUL8e«A(7) 

XRAMUE«A(lb) 

ISTARTaA(2<0 

ItNUaAOd) 

REWIND I 
NKAUU) 

REWIND 3 


RMAXaD, 
I ARsD 


C 

c 

c 

sad 0 

c 

c 


INSaO 


00 laa XalfIRANGE 

wRiTE(b,seaa) xtag.x 
FORraT(» SIG OOPP »,X6,» range ',163 


REAOm (ASl (J) « JalfNPULStl , ( ASQ ( J ) # J t i , NPULSE) 

c 

c 

C MULTXPLV bV THE REFERENCE SIGNAL SPECTRUM 
C complex CONJUGATE the RESULT FOR A SUbSEQUENT 
C XFFTV 
C 


C 


00 10 jaiSTARTd XENO 

XXia»si(j)*fil(J] * ASU(J)*F0(J) 
XXOa-*5t(JJ«FQCJ3 ♦ ASQUJaCXCJJ 
ASIlJjaxxl 
ASQU) a-XXa 


C 

10 

c****** 

c 

c 

c 

C NULL 


CONTINUE 


OUT THE NON-OVERLAPPING REGION OF THE SPECTRUM 


A30 


«30R QMl .f y 


c 


c 


c 

11 

c 

c 

2 t) 

c 

c 

c 

c 


c 

21 

c 

c 

Je 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 


u (i 3 Ti<}T,eu,i) boro eie 
IS115T*«T - 1 
00 II J>l»iS 

<SI (J)«0. 

4S9(J)«0. 

COMriNUE 


CONTINUE 

1T(IE^0.SQ,NPULSE} GOTO JO 

IfcalENO ♦ 1 

00 21 JaZE^NPULSe 

«SI(J)aO, 

ASQUJaO. 

CONTINUE 


CONTINUE 

COMPUTE TMg IFFT 


NUMaaNPUl.SC 

CAUU FFTtf (NUM8| ASif ASQ) 


sTorte THE magnituoe souareo away on oat slot * 3 

XMAXaO. 

IXMao 


DO AO Jai^NPUtSe 

c 

ASI (j) a*si (J) **2 ♦ »S 0 (J)aa 2 

IF (A 3 I (j) . lt.xmax) goto ao 

XMAXaASKJ] 

XXMaJ 

C 

uo CONTINUE 

c 

c 

6 T WRITE(J) (ASI (J) > Ja| ,nPUI.SE} 

NRITE( 3 «> 1000 ) I» 1 XM,XMAX 

WHITECb* 1000 ) ItlXM,XHAX 

1080 FOAMATC RAnOE AZ ', 16 ,' VALUE ', 1 PE 15 ,S) 

IF (XMAY .ur .HMAX) GOTO 100 
RMAXaxr^Ax 


lARal 

INGalXM 

C 

C 

c 

c 

100 continue 


c 

T»;>if'(|) aNPULSt 
TCMP(a)BlAR 
TtMH ( J) aiNG 
TEMP(4) aSMAX 


C 

C 


return 

Eno 
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c* 

C* Ftr^tFOH 

c* 

c* 

C« FfT *trfOTRITMM FUR VIRTUAL AKRAY3 
C* 

c* 

SUdROUTlNg FFIV(NUMttR,XT,yTJ 
VIRTUAL XTIDirTU) 

C 

C NUMbH IS The NOt O'- FOINTS OF eooo OATA 
C I^ NUT A PQMtR UF TWU RAO TO NtXT HlUHeST 

c MiTH Zeros 

c 

c 


1 

NaPONBI 

HEMPaNUpIBa 

ITEHP8lTeHF/2 



IF (ITeMF.LE.U 

GOTO 10 


NaP0WBN2P0M * 

1 


GOTO 1 
coNTiNue 


c 


NaaaaNaPQM 
IF (N.bE.NUMbR) 

GOTO IS 


NaN#a 

NRPOMaNaPOM ♦ 

t 

tj 

CONTINUE 



U (NUHdR.CetN) 

GOTO 20 

c 


NlaNUMBR t 1 
00 13 LH8NI.N 


c 


XTCLNJAOf 

YTCLHJABt 


c 

IS 

CONTINUE 


c 

c 

20 

CONTINUE 



c 

c 

HlNEPOW 
00 NOR LO*t»H 
CH-LO) 

Li*»a*LH* 

SCL’h.S'JSies/FLOATlLrr) 
00 6 UR LMaULT^X 
AnUa (UH-n asCL 
CaCUSCARG) 

SasiN(AMG) 

OU NRR LI«LIX«N|LIX 

Jl»Ll-LIX*LH 

JOaJlaLMX 

Tia»TtJn-*T(Ja) 

T<i«VTCJU-TT(ja) 

XT tJn*xTun*xT(ja) 
TT(jnBVTUn*YT(J 2 ) 

XT uaj aCaTJ ♦SaTa 
VT (jaj acaTa-5*Tl 
buo CONTlNUe 
Nva*N/a 

NHl aN«I 

J >1 
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r\a> 


ORIGIMAI i' 
OP POOR 


00 03^) UliNMl 

ir u.C'iitJ) SO ro ^il 
n«M (j) 

T'f'rT CJ) 

»r(J}»*T(i) 
v( u j •vTcn 
m (1) »n 

VTlU'Ti 
631 KINV2. 

6ti<9 CONTlMUt 

(K.SE.J) 60 TO 633 
JaJ-K 
KaK/e 
GO TO hii 
633 JaJ^K 

MtJMOfl aKt 
9tTU«N 

CnO 
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ow 6 »»‘- 1 ;*® 

OF POOR 


MlNUOVtPON 
VHTUAC VSKSION 
MlNUQwINli PUMCriON 


ARUUMgNTS 


AHrtAY IS THt ARXAT TO OK MlNOUWtO FIHST LOCATION 

ISTART is TnE FIRST SLENENT OF THt ARRAY TO St wlNOOWEO 
NUMSR IS THE number OF ELEMENTS TO SE MINOOREO STARTING AT ASOVE 
lOPT IS TYPE OF mInOOH 


0 NO WtNOORING CRECT FUNCTION ) 

1 HAMMINU xlINOOW 

a 9ARTLETT M|NOOW (TRIANGLE) 

3 HANNING RINUOM 
a slacaman riiNOOW 


SUHROUTINE WtNOOVCAHRAY, I ST AR T , NUMdR , lORT) 


VIRTUAL ARRaYU) 

c 

c 

IF (IOPT.EO.0) GOTO 10000 
XNANUMaR - t 
ISTAT«ISTART 
IF (lOPT.EQ.E) SOTO 200 
C 

Platt. «ATAN( i .) 
pxa«a. API 
pHASEaPlE/XN 

IF (lOPT.EQ.tt) GOTO 100 , 

C 

Ci«0,s 

ca> 0 ,s 

IF (10PT.E0.3) SOTO 10 
C 

CiaO.Stt 

CiaOtOE 

c 

10 continue 

c 

c 

00 ao lai.NUMttR 

e 

XNai ■ I 

HOaci - Ca*LOS(PHASE*XN) 

A-'BAY USTaT J awOAAHHAV tISTAT) J 
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c 

da 

c 

c 

c 

c 

c 


c 

uo 

c 

c 

c 

c 

c 

c 

2(4 «9 

c 

c 

c 


c 

2tia 

c 

c 

c 


c 

220 

C 

C 

C 

10000 

C 


OWGIMAL PAGH 
OF POOR Q JAUrY 


isrATiisrAT * I 
CONTiNUt 

Goro 10000 
CONTtNUt 
PH«S):2aPh«9e*2« 

00 120 ItlfNUMaR 
KNSX “ 1 

« 0.50*COS(PMASt*XN) * 0,ee*COSCPHASK2aXN) 
AltRAY (isra n awO*Ai4Rar (ISTAT) 

18TAT»ISTAT ♦ 1 

CUNTXNUfi 

GOTO 10000 

NUH0R IS THE NUMBEP OF EUEMENTS IN THE wlNOO^'l STAHTING 

P»OM the first address op AKRAr 


Continue 

IUPaNUMHR/2 

00 210 lal^IUP 

W0aFL0ATC2«Cl - 1J)/*N 
array (I STaT] anO* array (ISTAT) 
ISTATalST*T ♦ 1 

CONTINUE 

lUPaxUP ♦ 1 

00 220 lalUPfNuneR 

m Ft,0AT(2»(I - IJJ/XN 
array ( 1ST AT) aWO*ARRAY (ISTAT) 
ISTATalsTAT ♦ I 

CONTINUE 


SOTO 10000 

continue 
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PART III - SUBARRAY AZIMUTH PROCESSOR 


OT.POf? 


PROGRAM TO PROCESS THg SAW DATA SAWS P»0(iRAW 

Tils program PRUCbSS THt UATA USING THfc STEP TRANSFURM 


VIRTUAL ASt (4»9fc) > A3(3(«UR») 
VIRTUAL FG(«0V6J 
common ^InSTR/A(30} 
common /TAY/T(ia&) 
common /TfeMP/TEMP (la) 


CALL XNSTR 
CONTINUE 

IF (TEMP(«n ,gT,TEMP( 8) ) GOTO 10 
CALL GtN(FI,FQ) 

CALL UtRAMPCFIjFQ, ASlj ASO) 

CALL CORRCT (F1|Fq,a3I,ASQ) 

CALL FlNt(Fl,FQ) 
call 0uTP 0(FI,ASI, AS(i;FQ) 
TEMP(6)«TEMP(6) ♦ U 
TEMP19J»TEMPCRJ * TEMPda) 

GOTO I 

a continue 

Ave»TEMP(7)/TEMP(6) 

WRITE (6, laaa) temp {6} , TEmP(T) ,ave 
WHITE (5. 1000) TEmP(G) , TEMPI?) ,AVE 

FoRmAT(//» no. ',FS.0,' SUM »,IPE15,8/ 
I ' AVS SIPEIS.8//) 


ENO 
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ORIGINAL I A:.L 
OF POOR QUALITY 


DM, PON 


iNaTROCrrON StT PON THE S«R 

C IS THE speed OP U6HT • I, E » S ( H/SEC ) 

PI ■ a.«»TAN(|,) 

COMMON UCUCK INSTN 

*(|) THiNSMlTTEO PREUOENCV ( MHZ J PT 

*13J -*vei,tN6TH OP PT ( « J C/PT t 

*(53 VtCOCTIT OP SP*C6CR*PT ( KM/SBC ) V 
*(a) PANOli TO TARGBT ( KM ) R 
*(3) SamplxnU RATE ( MHZ ) PS 
*(6J PHP ( RHZ ) 

*(T) NO, OP PULSB8 transmitted N 

A(d) RANOWIOTH QP a NANGb CELL ( MHZ ) SW 

A(4) RATIO OP PS/Sh 

AtiJij TOTAL TiMfe OP Plight ( sec ) n/hnp 

Aun constant por the range walk crw 

( (1 ) 

A(U) range RESOLOnON RR 

C/(R«SW) ( M ) 

A(|i) CONSTANT PON the range AMPLIToOE punction 
PI/RaTIO 

A(ia) constant por the phase CORRECTION 
a*Pi*cRA/wL ( l/sec**2 ) 

A(1S) hamming weighting across th? antenna window 
YES ■ 1 NO ■ D 

IP YES OPTION OOESTIOn TYPE OP WINDOW 


A(1E) number op Range cells looked a7 
AUT) range >*AtR YES 10 NO ■ t 
A(i«) SANOWIOTH OP the system ( MHZ } 

A(19J 

A(20) NO, OF CELLS POR EACH DOPPLER FILTER 
A(Sn RANGE CORRECTION OPTION 

0 NONE 
t WU 

e CANUlANS 

1 STEP 


ttee) RATIO OP TMB sampling PREOUENCY to the lPm sweep preq 
A l8i) CORRECTION OP RANGE WALK a PT InTERP • i NONE ■ 3 
3 PT INTERPOLATOR • 3 3 PT INTER • 4 
Algu) STARTING SUSARRAY WHICH LOOK 

I IS FIRST LOO» SUB ■ 1 
3 IS ENO LOOK SUS ■ TOTaL/4 * I 


A(E5) SIZE PP THE FIRST PpT 

AlEb) NO. OP POINTS IN the OVERLAP OP SUBARRAYS 
A(E7) NO. OP SUBANRAYS por THE ENO PPT 
A(E8) SIZE OF the End OR PlNE PPT 
A(E9} NO. OP subarrays CONTAINING OATA 
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ORIGINAL PAGL lli 
OF POOR QUALITY 


*(3li9) MtXlMUM NUMaCft OF NANGC CKLLS f>RCC£SSE0 OUE TO mALK 

***************** 

COMMON f»U0C« TtHF 

******************* 

TtMM(i) 3TCMES TH5 Nu, OF FINAL AIMUTH SAMFUES 
TfeMP(?J TmE location in range OF THE MAX 
TbMP(3) TmE location in AZIMUTH OF THE MAX 
TMEP(a) Tm6 maximum VALUE 
TtMPISJ 

TtMP(fe) THE NUMUes OF OIFfENENT STARTING TIMES HAVE OONE 
TbMP(T) TmE running AVERAGE OF THE INTEUFATE MAINLOBE TO 
TmE 310ELOSE FOh OIFFEmeNT aTARTiNG TIMES 
TtMP(e) TmE last fiMg FOR different siawting times 
rtf'FOI TMC RUNNING OELAT FOR DIFFERENT STARTING TIMES 
TcMPtlS) THE increment for THE DIFFERENT STARTING TIMES 


subroutine instr 


common /INSTR/AI30) 
COMMOs/TAV/TArLOR(lJ 
COMMON /TEMP/T£MP(IJ 


CALL ASSIGNC6, *TTl») 


RRITEIs,!) 

FORmiTC NAm£ of LINEPRXNTER OUTPUT FILE *f 
I ' LPJ OROSNAME MTtJt NAME */J 

Call ASSIGN!?, *lpi ') 

CALL ASsir.Nia, »outp 3 ,sar', 0 , '-jNRNowN', ,a) 
Call ASSIGNC3, »OUTP3.SAmS0» 'unknowns »23 
CALL ASSIGN u, 'OUTPa,SCw',0,'SCH»,,a) 


name of disk file to 8E processed 


w»lTECS,5J 

Fohm*t{* name of 0I5K FILE TO BE PROCESSED '/) 
CALL ASSIG Ul , 'DUMMY', -1) 

CALL ARSI6N(t, 'ST12«SAR',0) 

READ IN THE CONFIGURATION 

REAO(I) A 

c 

c 

A(as) aSB, 

AianaSB 

A(PB)al 28 

C 
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3 

r 

s 

0 

9 

0 


la 

0 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

loaa 

i 

i 

i 

tt 

9 

6 

C 

C 

tala 

i 

a 

5 

a 

5 

b 

C 

C 

1020 

1 

i 

i 

a 

C 

c 


^ftlGINAL 
OF POOR QUAUTY 


«HT rf (<S(bJ 

8126 or TMt MW8T 6FT »/) 

AfeAi}(b«/) «(29) 

FllNMlt 10,0) 

Ill'll riu&i a) 

6o0m»t(' no, of points of The OveWLSAF '/) 

*(2b) 

M4trk(b«9) 

FO»N*T (• NO, OF SUdAQNAVS USED IN 2N0 FFT »/) 
Pt4i)(b,7) A(27) 

MNirii (b,)0) 

FO«m4T(* no, of points in 2nd FFT */) 
PbA0(6«T) «(2S) 


WXITE(9| 1000) (A(I)>lal»7) 

FOH'<*T(' TPANSMITT60 FSeOUENCY ( mnj ) ',tPEtS,8/ 

' WAVEUeNGTH ( M) ',}P 61 S,a/ 

' vetOtlTT OF THfc SPaCEC«IAPT ( KM/SEC ) 'flPElS,8/ 
' HANGfc TO THE target ( AM ) »»IPE19,8^ 

' SAMPLING MATE IN PANOe ( Mh2 ) ',1P£15,8/ 

» PRF OF TMft A2IMOTH I KM2 ) ',lPtl9,8/ 

' NO, OF PULSES UATEO OQwN ',F7,I/) 


wpir6(s»i0i0) (A(i) ,i>a,i4) 

FORMATC BANOwXOtH OF A RANGE CELL ( «RZ ) ^,IPEIS,8/ 

' RATIO OF TMt sampling rate IN RANGE TO BANOhlOTH OF A *, 
' CELL ',lPfcl9.8/* total TIME OF COOK ( SEC ) ',tPSt3,8/ 

' range RALA constant ( M/StC»«E ) *,lPElS,8y 
' RANGE RESOLUTION ( M ) ',lPEl>i,8/ 

» constant ►or AMPLITuOt FUNCTION *,IPE15.8,» ( PI/KATIO ) 
* constant for quadratic phase '|IPE19,8» ' (SECaa-2) »//) 


RRITE(S>1020) (A(I) ,Ial9,lS) 

FORMAT!' HAMMING WEIGHTING ACROSS ANTENNA APERTURE •, 
' YES a J NO a 0 ',Fb.l/ 

' NO, OF RANGE cells TO PROCESS ',FS,1/ 

' range gala yes a 0 NO a t '»F5,1/ 

' 0ANOWIOTH OF THE SYSTEM I MM2 ) 'ilPElS.S/) 


OELLFMa2,aA(3)**2aAa0)/(A(2)*A(4)) 

0LaA(«i)/A(r) 

A (20) tOELLFM^OL 
A(20)aA(20}/4, 


C 

C ABOVE assumes we have FO'.R OOPPLER FILTERS 

c 

WHI TE(9f 1030) OELLFM, A (6) >OL, A(20) 

1030 FORMAT!' CMANGE in FREQuEnCt ACROSS WINDOW ( AH2 ) '» 

1 IPE19,0,' SAMPLING RATE ! A«2 ) », 

2 IPtlb.a/ 

3 ' NO, OF NmZ per FFT sample '»iPeiS,8, 

4 ' NO, OF CELLS PER OOPPLER FILTER '.lPE15,8//) 

C 

e 


*! 
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Acat) 

lv>44 FUt)^<AT(' HANGt bO*«»<eCTlON FaCTQH none • a WU a I *1 
I ' C*^AWIANS a e STE»» a i ',^9.1//) 

c 

c 

c 

e 

e 

c 

c 

c 

c 

c 

MwiTE(>>» taea) 

IffeO FUI9 maT(//» TfK£ OF /JlNOOW ACB05S feAC-i APE9TUK6 */ 

1 * I a HAMrtlNU a a 8*BTI.tTT 3 • HANi'flNtt */ 

a ' a a ucacaman 8 as os tavcor*# 

3» 6 S<3 08 TATUOK T 35 Ott TAYLOR'//) 

c 

*(l5)ai, 

0 RtA0(6,T) A(19) 

C 

wf»lTE(5, 10T0) AU5) 

iOTO FUR*^AT(//' type of WlNOUM ACROSS EACH SUBARRAY '/ 

I f 1 hamming a RARTLETT 3 MANNING *, 
a ' 4 a 6LACRMAN b 35 08 TAYLO«'» 

9 ' b 3a OR TAYLOR T 35 06 TA YLOH ' , F9 , I/) 

c 

c 

c 

MHXTECi. (A(j) , jaab.aa} 

tUSM F0 Pm4tj,vp mo, Qf COEFFICIENTS IN FIRST FFT ',F6,a/ 

I ' NO, OF POINTS OF OVtRLAM OP 8UUABRAY3 ',F&.a/ 
a » NUMBER OF SUBAHRAY3 USEU IN FINE FFT ',Fh,2/ 

3 * NO, OF POINTS IN THE SECOND FFT 

C 

c 

loaoa CONTINUE 

c 

A(a3)aA(SJ/0£LLFM 

wRiTE(5«3aa0) A(sa) 

aoeo FORMAT!' RATIO OF sampling FREO to LFM hate 'rlPEIS.a/) 

c 

c 

c 

c 

WRITEC6»10R0) 

ia90 FORMAT!' CORRECTION SCHEME FOR RANGE RALK '/ 

1 ' a PT INTERP a j nUn‘E a a 

a 3 PT INTERP a 3 3 PT INTER a 4 '/) 

*(33)ai, 

«EAO!b,T) AiaSj 
HPlTE(b, 1095) 

J095 FORMAT!* SUBARRaY STARTING POINT «HICH LOOK '/) 
A(341al, 

REA0!b,7) A!34) 

C 

C 

w»ITE{5. 109E) A!a3),A!34) 

I09b FuRmaT!' CORRECTION SCmEME FQR RANGE 1 a A PT a a ru f 
I F5.3/' SUOAHHAY starting RQINT »»FS,3//) 

c 

c 

r, CALCULATIONS FOR THE DIFFERENT STARTING TIMES 
C 

c Slope time over lfm expanoeo 
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Wlffi-Mi, J, 

Poor; Q,. 


C 

C 

C I^CHtMfeNf e*cn 0 EU 4 V TIME 
c 

C 8 LU ‘'6 • >»'»»'/ (Size UF 1 ST Ft-T • SIZE OF THE Snfi FFT/t) 

c rnt /E IS uo ro uvexsampuing paooing »nH tienos 

C TMtCE ueUMTH 

e 

Tbi^FUeiaStOFtAAI^I/CACESlAACEai/E.) 

c 

c 

c STAMTING TIME 

C 

C TriR /« IS Fon TmE PAOOZmG hIZh ZEROS 

C 

C 

C AMU EMUZnG riHE 

c 

c 

re><P ( 8 ) • TEMP ( 10 ) •A( 26 )/ 4 , 

WWI re(6»HJ97) 

1 b 9 Z FoBMATC starting time C 0 - up )*/) 

«Xl> 0 y 

ReA0(8*T) xil 

TEMP{9) a-TEmP( 8) ♦ X*I*TeHPU0) 

c 

c 

C PUT In for EXPEo 

c 

TeMp( 8 )« 0 , 

TeMP(ia)arEMPa 0 )« 8 , 


c 

c 

C RUNNING sums and counter 

c 

c 

TEMP(T)»0, 

TEMPt6)»0, 

C 

C 

WRITE(S,IE0P) (TEmP(IJ,I»e,10) 

U00 FO»*^AT(' COUNTEHS S 2(U.fS.0)/ 

t ' ENUING ANO starting times * , E (1 X » 1 PEIS , 8) / 

i • Increment each time %ipeis«6//) 

c 
c 
c 

XSTaAOS) 

IFFTiaA(rfS) 

IMIST.GT.0) CALL WT(TAYL0RaP^Tl,I5T) 

RRlTE(b«7a4t) A(30) 
a«ITet5,T89l) A(30) 

I6<tl FORRATt//’ NO. OF range CELLS PROCESSED *LL TOT '.FS.8J 

C 

C 

C 

C 

RRITE(6,SaS7) A(83) 

3«S7 FoRMAT(//» type of INTERPOLATION »,F5.1» 

I ' I A PT a NONE 3 2 PT « 3 PT *.V) 

WRlTe(3.30ST) A(23) 

C 

C 


RETURN 
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C 

c mcK2,i*0H 
e 
c 
c 

c****#***.*. 

c 


SUDNOurXNt GEN(FI,rO) 

c 

C*«*<****«t«*«****** 

c 

coi“«ON /i^srt</A(i} 

COMMON /'TEMP/ r£MP(l J 

c 

e 

C GCnEHATE out the PEPettENCg SIGNAC 
C ^OH MtrcMEO FIUTEMInG in the anuce domain 

C CINEAH FM WAVE TQNM 

c 

e 

c 

c 

NPUCSKaAIT) 

TSTbP»l,t-3/AC6) 

0ELAf»rtMP(9J 

TlMt»-4(lCV)/a, 

TMNACA-TIM6 
TIMEaTImE ♦ OeCAY 
7f>lNAC«TFlNAC ♦ oeCAV 
CONSAAdO) 

c 

Wi<l 76 ( 9 / 1000) TEmP( 6 ) ,TSTgP,oeCAV,TlMe,TFINAC 
WMI Tg ( 6 , 1000) TEMP ( 6 ) , rsrEP,oeCAY, TIME, TFXNAC 
1000 FohmatC//' step /lO, TIME STEP 'ttPEIS.S/ 

I ' CD9H3P0NUING OELAT (SEC ) 'rlPElS.S/ 

a ' initial time ',»PE19,8,' final time 

c 

c 

c 

c 

I CONTINUE 

c 


IF (I.GT.NPULSE) GOTO 100 
C 

c 


TSOiTlME*TIMe 

PMASEaTSNaCONS 

AXaCOS(PmaSE) 

FI(I)«»» 

«i«SIN(PmasE) 

F3(I)aXX 

C 

TIMEaTIME ♦ T 8 TEP 
lAl ♦ I 
GOTO 1 

c 

c 

100 CONTINUE 

c 

s 

® PETUPN 

ENO 
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ORIGINAL PACE B 
OF POOR QUALITY 


OlCKS.POfi 


SUttMOUriNE 06R*MP(^I,f0,SX*3Q) 


viRnjAL Pi(t]«Fon)>si(n» 3 Q(i) 

COMMON /INSTR/40) 


NPULse«*(/) 
IWlNOa* ( 13 ) 
NMANCEa* (33) 
iHANCeaAClb) 


RgWXNO 2 
RtMiNO I 
KEAO(l) 

H£<10(1] (SI (J) , J*1 ,NPUUSt) 

MHITECb.2341 Sicn ,S1 (3 1 2) « SX ( 1324) « SX (33aa) » 
1 SI (3U<*b) ,31 (409b) 

34 FORNAUiU*, IPetS.S)) 

00 300 1*1,NR4NGE 

GET THt SIGN4L FROM TM£ ITM RANGE CELL 


REAO(I) (SI U) , Jai ,NRULSE) , (SU ( J) , J* 1 , nPUCSE) 
OERAMP THE SIGNAU 


00 3 J*I,NPUt.3g 

XXI* ♦ S:(J3»PI(J) ♦ SOtJ)«FO(J) 
XXO* - 3UJ)*Fa(J) * SG(J)*FX(J) 
SI (J)*XXI 
S0(J) «XX(4 

continue 


PERFORM The FF? store RESULTS AWAT ON DISK FILE ♦ 2 

USE . A S3 X Q'»!gRtAP 


PUT Tm| phase correction on oue to overlapping appertubes 

call 0VER(SI.«SQ) 
c 

330 continue 

c 
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0ICK3*,K0R 


v9 0#*fl***t*4 

SU9R0UTINE 0VEi«(9Ii3Q) 


VlRTUAt SI(l)>30(t) 

A t (Ud) , «U(l2a) 

CUMMUN /1N3TR/A(i) 

COMMON /T4t/TAVC0H(l) 

Pe^iFOHM THE FFT ST0n£ RESULTS AnAY ON DISK FILE ♦ 2 


PUT The phase CoRHECTIOn ON OuE TO oveplappino appertures 

CORRCTIUN FACTOR IS 

EXP(J*2«’?I«HATI0 OF overlay*arrav no. - I) 


NPULSe*A(7) 

lLeN«At23) 

ILEN2«A(2S) 

lMlN0«A(i») 

ISai 

lEalLEN 

NOMSaO 

overlap factor for phase cornection 

D«A(3b)/A(2S} 

Pl»«,*ATANa,) 

piaaa.APi 

P120*PI2*0 


CONTINUE 

Kai 

OU 10 J«IS»l£ 

AI (K)aSl(J} 
AQ(K}iSQ(J) 
K«K ♦ 1 

0 CONTINUE 


ISaiS * ILEN2 
lEalE ♦ 1LEN2 
C 
C 
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c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

ir 

c 

c 

e 

c 

c 

c 

c 

e 

c 


c 

la 

c 

ae 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


fur wlNUOrilNG AC*<0SS THE SU6*))I4AYS 

IfCwiNO.Ea.a) GOTO IT 

ou la K*i,it.EN 

A1 (K}aAl(K)*TAYtUH(K) 
AUCOtAaCKJaTAYLURCO 

continue 

CONTINUE 


CAUL FFT(1(.EN,AI,AU) 
lK(NUMe.£Q.,U) GOTO 30 

PHASe»PXa0»Ft.0AT(NUM8) 


00 10 K*a»Xl.EN 

XK«PhasE*FLOATXK-IJ 

SNaSXN(XK) 

CSaCOS(XK) 

X*laC;i»AX (K) ♦ 0NaA(J(K) 
XXGaC0*AQ(iO - SN*Al(K) 
AX (A) aXXX 
AOCKJaUXO 

Continue 


CONTXNOE 


sto»(E amav on Oat slot * a 

HBXTttE) (AICH) .Kai.UENJ , ( AQ (K ) , K a J , XUEN) 
NUtiSaNUM0 ♦ X 

XFUE.LE.NPULSej GOTO t 

A(d 9 ]aNUH 8 

RETURN 

ENO 


A45 


oonoM o ri o o o o r>o o o o o o o o o o rtn o onoo n rtn r> 


ORIGINAL IPA3L’ c3 
OF POOR QUALITY 


OTb.rOM 


COKft&CTIOl^ OP MANGE walk XN THE STEP TNANSPORH 


SUePOUTiNE CORHCTCHANCEPf A1,AQ,0UHHV] 


VlPTiJAU AI (120,332 «A0nEdii32} ,PANGEK(i) »0UMHY(1) 
PtAL FF(2S<») 

COMMON /INSTR/A(l) 


NO PANUe COTMECTION ? 

iF(A(iT),EQ.i.) GOTO tooea 
f*ATlO« .(22) 

RM0v«BATI0*A(25)/A(a7) 

NHULS£aA(7J 

MAIS0a*A(29) 

MXSMAXSUS - ( 

MINS JH»l 

NMCaEK»A(25) 

nm 2 »NmC 0EP/2 


hEmIno 1 
pEAoa) XXI 

RcAO(l) (OOMHV CD » I»i *NPULSE) 


PEOPOEM THE APRAV 

Pino the min 

XMlNwOUMHYCn 
L»NPULSe/2 
00 3 l>I,NPuLSe 
L«L ♦ I 

IP (L.GT.NPOLSE) L«1 
PANliERU) aOUMMYCL) 

IF (XHIN.GT.UUMMY (L) 2 XMIN»0UMMYCU 
CONTINUE 


Pino the location op the ieho range cell 

LMN«1 

IP (XMIN.GE.O.) GOTO « 

XS»-XMIN 

* JIMIN» 1 P Ix (»S) 

Il» (XXMtN.Nb.XS) XXHINaXXMIN ♦ I, 
LMNl-XXMlN ♦ I, 

4 LMN«LMN - 1 

C 
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OHIGtNAL PA ^■p , 

fVOR quaut^ 


NMANGb«*(3'3) 

lH*NGbaA(tO) 

lSCt;MbaA(« 3 ) 


Rbi^I^O 4 

00 ia»i iat>M*xsua 


NEAO IN THE 1 TH SUQARRAY fOU EACH RANGE Cbtl ANO THEN 
00 INt CUMRtCUNG WRITE THE RbStJtTS OUT TO ♦ 4 


rewind 2 

n«i - i 

irdi.iT.n GOTO le 
00 s J«UII 
REA0(2) 

CONTINUE 


DO 2ia Js1,NKANGC 

RE AO (2) (AI (K,J] ,Kal,NHC0EF)» (AQ(K, J) ,Kal,NHCOEF) 

If (J.EQ.NHAnGE) goto 29 
00 12 All, MX 
REA0C2) 

CONTINUE 


00 The correction whole ABhaY at a time 


c 

c 

c 


STARTal 

XRNGaNR'^NGE « IRANGE 
IRNG2aINNG * I 
IRNGLalHNG 
IRNGalRNG • 2 
lRNGt3*IRNGL - I 

ISTRa-CMN ♦ 2 

00 39 Jat,NMCOEF 

ISTARTaSTART ♦ 9,S 

If tSTABT.LT.a.) ISTAHTaSTAHT - 9.3 
If (iST*RT.LT,n GOTO 35 
JtLH»NMe;*lSTART 
HANwUAaWANGER (lELM} 

IFClISCtME.tC.n, OR, (13CEME. £0,3)3 GOTO 24 
If CISCbfE.EU.a) GOTO 24 
IHI.IAHS IRANwtK) * 9.3 
IF (SANwtA.LT.a.) IRwa-IR|i( 

«i<«:.sth ♦ iRw 
IFlAA.LT.n GOTO 28 
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c 


se 

24 

C 


C 


c 

c 

2 b 

C 


C 


C 

9(480 

C 

C 

c 

c 

c 

9(461 

c 


c 

9(34(9 

C 


!(• (KK.ar.lHNUt!) KKa|HMG2 
00 22 K«l,lNAN0e 

Al(J,K)iAl(J,KK) 

A0(J,K)aA(3(J»Klt) 

KKiKK A 1 
CUNTlNUt 

cure 28 

CONTlNUt 

iKIWifiANi'^tK 

oeuTAiHAiM'^cK - FuOArtXHW} 

K«I3TH t 14W 

IFdSCtHt.EQ.n XSK - I 

lF(lSCtMt.6U.4) K*K - t 
IF lObUr A.UT.a.) KIK - 1 
IFCO&UTA.UT.O.) ueuTAat, A DELTA 

IFCK.UT.l) «•! 

Ual 

IP(lSCEMt.E0.3) 6(JTU 6083 
IF(If’Ct'iA.,gQ.4) OOT(3 6(940 
IFIK.GT.I'^'^G) KaiHNG 
ueLTA2a3tLTAAOtLTA 
OtLTAjaOELTAaaoeLTA 

AMta.ObLTAA(oeLTA2 <• 3.aOCLTA A 2.)/b, 
AOa(ogLTA3 - 2.AUELTA2 > DELTA A 2.)/2. 
Aia-oeLTAA(OeLTA2 - OELTA - 2.J/2, 
A2a0bLTAA(0ELTA2 - i,)/6. 


CONTINUE 

VALUEIaAMltAI (J,KJ a A^aAICJ^KaU a AlaAICJ^AAi) 
I &SaAI(J,HA3} 

VALtjEaaAMl*AC3(J,X) A Ai) * AQ C J , K A } } A AtaAa(J,KAa) 
1 A2 aAOU»Ka3) 

A1(J,L}«VALUEI 
AG(J,U ‘VALUEQ 
L«L ♦ 1 

KaK A I 

IFtL.Lb.If’ANGE} GOTO 2b 
GOTO 28 

CONTINUE 
2 FT INTERP 

0ELTA2ai. . DELTA 

IFCK. GT.INNGL) KaiRNGL 
continue 

value I»0ElTA2aaI CJ.A) ♦ OELTAa*! IJ,KA1) 
VAtUt:JaOELTA2*AQ(J |,K} A OEL T A A AQ ( J , K A 1 } 

AI IJ.KjaVALUEX 
AOU,A)aVALUE(J 
LaL A 1 
KaK A I 

IF CL. LE.IxanGEJ goto 6081 
GOTO 28 

C(JNT1nuE 


A 

A 
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g 

AMlaOEuTAaCDELTA - l,)/ 2 . 
AOaj, . UELTAaaB 
AlaOEtTA*(OEtTA * l.]/ 2 « 

PAS2 iS 

OF POOR QUALITY 

8091 

C 

CONTINUE 



Vi(.UbI«AMl*«I (J,K) * »l3*AI(J,K*n ♦ 
VAUljeC9«Aril.A0(J,K) * Ai)*AQ(J«K*l) * A1*AQ(J,K*2) 
AKJfUaVACtJbl 
AQ(J<C)«VALUEU 
L*(. * I 
* KiK « t 

xi'a.tt.ittANGe; goto eo4i 
c 
c 
c 

28 CQNTINUe 

c 
c 
c 

5TABT»3T*»T - UMOV 

c 

30 COMTINUe 

c 

c 

c 

c 

c MRxre our rnb newucrs to ♦ a 
c 

39 CONTXNUe 

c 

e 

00 4V) KaXflRANGC 

40 wRXTftCa) (AI(J,K),jat,NMCOEF}, (AQ(J,K),Jal,NnCO£F) 

c 

too CONTINUE 

c 

c 

C WRXTe the file 8AC4 TO ♦ 2 In THE FORM IT KANT8 

C 

c 

ISKIPaIRANGE » I 
MMagaNMCOEF 
REhINO 2 

c 

00 200 KAlflRANGE 
REWIND 4 

IFCK.EQ.n GOTO 110 
IIAK • i 

00 los 1*1, n 

105 REAO(a) 

110 continue 

e 

c 

00 190 lat|MAXSUB 
C 

REAU(a) (FF CJ) * j*l (NM) 
wfiITE(2) (FFCJ) f J"1»NM) 

IF U.EQ.haxSUH) SOTO ISO 
00 lao jalflSKlP 
IRO REA0(4) 

c 

ISO continue 

c 

c 

200 CONTINUE 
C 

10000 CONTINUE 

c 

RETURN 

ENO 
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OF POOR QUALITY 


UTSt^ 0 » 


SUbROUTlNtl MNE (PX,FO) 


VI 9 fU*W FX(n>Fi](l) 

common /CO&FS/AI (1924) 
common /lNSTR/t(l) 
common /TtMP/TEMPCl) 
COMMON /MtM/H*M( 296 ) 


COMPUTE The fine mesouutxon FFT 


IRANUEaA( 14 ) 

9 *TI 0 .A( 22 ) 

NMCOtFtAOS) 

NMEaNMCOEt^ *2 
IFINfe«A(27) 

MAXSUaaA (29) 

MlNSueat 
XPAUaA (2d) 

«M0Vaa4Ti0.A(7)/(A(2b)aA(25)) 

STARTSal , 

IP (A (24) , to. a.) STAHTSaA(T)/(4,aAcae»)) ♦ 1 , 
WMI rtvt. 789) 8TAKTS»A(2a) 

F0HM*T{' STaRTXNG 8U8A9WAY % FS , 2 , 1 X , PS ,2) 


set The hinoOM WEIGHTS 
ISTaX 

CALL wT(ham,IPINE,XST) 

w»ITE(6,S00a) IHANGE.NMCOtF, IFINE, XPAO, 

1 MIN8U8>MAXSUS 

5000 F0»maT(//» pine IHANGt ',lb/ 

I ' 1ST PFT %I6.' 4 OP SUBS IN 2ND *,IS» 

a » 2N0 PFT ',18/' min SU8 ',X6,' MAX SU8 ',16//) 

c 

wwlTEC6,500n RATIO, 9M0V 

5001 FOSmatC' patio of SF/LPmF ',1PE15,8/ 

1 ' incpement of sue ',iPti5,a//) 

c 

RMAXa0, 

lAXae 

XRNGa0 


ooorioooooooo o n lu or**- oo(-io oo ooooooo ooooo® r?ooo 


f*OOH QUAL/Ty 


c 

c 

c 


RCmIM) 2 
rCminO 3 

ou eats Kit » ifiANcE 

REMIMO a 


ratNiftH DATA FOR THE K TH ttANQE TO .DAT SUOT ♦ 4 

00 4 tit.MAXSUt) 

aeAo(d) tAi(L).tit«NM2) 

«>R;Tt(4) (AX(U.Ci|»NM2) 

CONTINUE 


00 TME coefficients t TO HALF ♦ I 


lARiO 

CALL P0SHUF(FI,MMUV,START8f lAli.COFNAX.ICOF) 


00 me NbGATlVE COEFFICIENTS 


CALU NEGHLFtFI,HMOV,STA«T8,IAR,C0FMAX, ICOF) 


».RITE(3) CFKJl.Jit.lAR} 
WRITE IT OUT TO TME wP 


wRiTEis. leae) icof.coFmax.k 
ICOF.COFmax.K 

FORI^ATI' MAX IS at M6,* value OF » IPEtS.B.' RAN ',16/) 


IF (COFMAX.LT.HMAX) GOTO EltfO 
HmaxiCOFMAX 
lAIiICQF 
IRNGiR 
CONTINUE 


TEMP(1)iIa 1» 
TEMP(E) iIHNG 
TEMPtSJUAZ 
TEMP(U} iRMAX 


TEMP(l)il4a0, 

TEMP(a)ia, 

TEMP(3) i3i, 
temp( 4 ] la.aoRAASiEES 


return 
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0T*J*,F0i» 


5Ua»0UTlNg POSHLY (FI,9M0V,5T4*:rS, £*P,COPiiA)l,ICUP) 


V1BTU4*. FI(t) 

COMMON /COdK'i/ Al (138) ,*U(t3A} iWOMKI (3S&) fKiOPKO(aSb) 
COMMON /IN3TH/ACU 


NMC0tF«4(3b) 

lFINt»*\(3T) 

MAxSU&"A(3V) 

MlNSUM'i^l 

IPA0cA(36) 

NUMdRalHA0*IFlX(A(3«il)/NMC0fiF 
WilITt(<»,U5b) NIJMOS, IPAOiNMCOgP, A(26) 
F0HMAT(S(Ib,U],tP6l3.8) 


00 THg coepFicifiNTa I TO m»i.p • i 


ICConmCOET/3 ♦ I 
PUNS0M39TART3 
i3TAHT»START3 
ItNO'lSTAHT * IPlNg - I 


cofmax»o, 

ICQPaO 

00 30 l*tfXCC 

IP (ISTAPT.GT .MAXSUb) GOTO 60 
RbwiNO a 
1G«1ST4RT - 1 
IPdG.LT.n SOTO 10 
00 S JalfIG 
Rt4u(a) 

0 CONTINUE 


IE»ieN0 

IP (ieiGT.M4XSUB) l£aM4XSU8 

c 

JJ»0 

00 13 JaiSTAflT.ie 

c 

P8 40U) (41(L)>Ual>NMCU£P), { 4Q (U) , I. • 1 « NMCOEF) 

JJtJJ ♦ 1 
MOMK 1 (JJ] 141 (I) 

WOPPQ (JJ) ■40(1} 

c 
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cQNriNu€ 


' ^ ^ • V , I 

’ \t ■ 

» W / <» ' « 


OF PCO^'- 


i!> 

C 

c 

c 

c 

c 

c 


c 

IT 

c 

c 

iv 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

e 

c 


c 

c 

3e 

c 

c 

c 


Put IN TRHOS WHtSt NO T*«(.€T ft*18T5 


{F(JJ.eO,ZTlN£) GOTO 20 

jvajj ♦ 1 
00 ir JaJXilFINg 
MORKl(J)lO, 
MORXQ(J)ae, 

coNTiNue 


CONT INUK 


PUT MfiMtNO xtRIGHTINa ACROSS SUOARHAVS 
NULL TH^ RbST OF AHRAV 
COMPUTt TMt ►INE structure FFT 2NQ 
MNO TMb magnitude Op THE FFT AU0N6 UAR6EST lECT 


CAUL PA0FFT(F1,C0FMAX, IC0F,IPA0| IFtNE.IAR^NUMSR) 


RUNSuS’RUNSUS * RHOV 
1ST ARraMUNSUO * 0»S 
ItNOalSTART ♦ XFlNg - I 


CONTINUE 

CONTINUE 


RETURN 

ENO 
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or^u.rot 


SUbROUTt^e NEGbLr> (Pl,RI^OV,sr*RTS, X«ifl,COFMAX»ICOR) 


VIRTUH Ficn 

common /cotes/ *I(Ua]«* 0 ( 12 e),l^ 0 HKH^S 6 )»W 0 RK 9 (aS 6 ) 
common /IMSTR/*(1) 

NMC0tFa*(2S) 

IF |Ntat( 27 ) 

MkX90t5*«(2V) 

MlNSuSat 

lP« 0 «a( 2 a) 

NOMdRaIPA 0 «IPlX(*( 2 b)}/NHC 0 &F 


00 The negative COtFFICXENiS 


l(*aFU0AT(NMC0EF/2 - DaSHOy 
RONSUOaSTARTS • XX 
AxaABSCRUNSua) a O.S 
XF CR0NSUS.tr. 0.) xxa-xx 

XSTARTaxy 

XtNUaXSTART a XFInE • 1 


lQaNMCUEK/2 a 2 

00 lOVI XaXO.NMCOtP 

REwXnO a 
ISaiSTAKT 
IFCIS.tT.l) XSal 


XFlXfcNO.tT.n GOTO te 

JJa0 

IFUS.bO.lSTART} SOTO 10 
IOlFea»ISTART a \ 

00 5 Ja;»X0lFF 

W0RAI(J)ae, 
MORKO(J)a0. 

continue 

JJaXOXFF 
C 

10 CONTINUE 
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Q 

C 

C 

DO 

Re AO (4) (*I (U iL*t r^MCOEF) > CAU(U) «l«i iNMCQCF) 

JJeJJ • 1 

WORrtQCJJJiAOCn 
59 CONTIMJE 

C 

C*tu PAUFFTCFliCOFMAXilCOF, IRAQ, tFlNe.XAR.NUMBR) 

e 

C 

6ia CO»tTlNU£ 

c 

c 

e 

RUNSlJtiaRUNSUa * HMOV 
«XaAa 3 (PUNSU 9 ) * 0(9 
IF (RuNU'i9,LT,0.) XXa-XX 
X3TAR1 aXX 

lENOalSTART • iFlNfc - 1 

C 

100 CONTINUe 
C 

c 

c 

RETURN 

feNO 
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OOUOOOOU OOO 0 1»00«J O VO O O Cl €9tnOOOOO o «uo«uooooo oooooooooooo 


original page is 

OF POOR QUALITY 




auaROUTlNI?. P*Ol'FT(PI»CUFn*lf !C0F,IPA0, JJ, IAR,NUH8R) 


VIRTUAL Fl(t) 

C 0 H* 40 N /C 0 bP 5 / Aia 2 a)^Aan 28 )»w 0 Rill( 896 ),M 0 RK 0 ( 296 > 
COMMON /MAM/HAM( 1 ) 

PUT MtMMlNG wblGNTING ACROSS SURARRAVS 
NULL THg HfcST Of AHRAV 
COMPUTt TNt MNe structure PPT 2 N 0 
PINO Tne MiGNtTUOb OP THt PPV ALONG LARGEST »ECT 


PUT THt WEIGHTING across THE MtNOOM 

00 S JatiJJ 

aORNI (J)3nAM(J)«M0RKl(J) 
wORKQU) aHAM(j) *mURKU(J) 

continue 

NULL OUT THE rest OP THE ARRAY 


IP (IPAO.LE.JJ] GOTO 23 

JfJJ * I 
00 2^ JaJXilPAO 
MORKI UJae, 
WONAQ(J) a8, 

CONTINUE 

Continue 


compute the pine structure 

call FPT (IPAOf wORKI,mORK9) 


Compute the magnituoe op vme souareo 

REORCER THE ARRAY 

SINCE half the fine COEFwicItNTS ARE TO THE LEFT ANO 

The other HftLP TO T"E RIGhT FOR EaCh GROSS COEFFICIENT 
Over sampling BY factor OP a to I thror out half the coeff 
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MICROCOPY RESOLUTION TEST CHART 
NATIONAL BUREAU OF STANDARDS 
STANDARD REFERENCE MATERIAL 1010a 
(ANSI and ISO TEST CHART No. 2) 








noon 


original 

POOR 


PAQS is 

quality 


c 

c 


c 


c 

a<t 


JX^lHAO - JU *\ 

DO 6H JuliNUMDN 

xxtwORKX (jx) **2 ♦ MuaHa(jx)**e 

JX1JX * t 

IFlJX.UT.IPAO) Jxel 

IAR*X*M ♦ I 
PI(UH)*XX 

iPCXX.UT.CUPnAX) GOTO 24 
COFMAXaXX 
ICOP*UN 

CONTINUE 


RETURN 

tNO 
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c 

e 

e 

C OtCN<l«l>OR 

e 

c 

c 

c******t*****« 

c 


sutiHuuriNt oori«0(*ti«2i*3«*4} 

c 



c 

VICTUAL 

CU*4M0N /IN&TN/**(|) 

COMMON /TlMC/teMPtn 

c 

c 

c pAif^r our TMt o*r« sctct ro tme mak pa!iseo in arguments 
c 

NHANGbaAA(iei) 

NPO(.3birEMP(i) 

iA«reMP(i) 

XMAXlTEMPCA) 

e 

NNITt (b, Jg(50) IN, JA, XMA*,NP(Jt.SE 
MMITcCS. toIJit) Iff, lA, lMAX,NPUUSe 

lOiSO XORMAT(//' NANGt cetu ',I6/' ANGCt VALUE *,IPE15, 

I • NO. PULSt S16//J 

C 
C 
c 

C FINO THE LIMITS OP THE MAInlOUE ALONG TmE At’.IMUTN SLICE 

C 

CALL SLlCfc(iR,lA,NPULS£, ILF , T M T, XMA A , At) 

c 

c 

c REUPGANIZC the FtLt IN THE AElMUlTM VERSUS RANGE 
C RaTmEN then »N THt NANGt vEhSUS azimuth OXM 
C 

C also put cm pC^n in THt MlOOLC OP THE PILE 

C 

C 

CALL HEORO(NPULSE ,Ia,nRANGE, A t , A2, AS) 

e 

c 

c 

C find The sum IN THt mAINLOsE ANO The SIOlOOE FOR 
e HANGt ANO azimuth OIRECTIONS 

c 

c 

c 

CALL VALUt (NPULSE.IR, lA, ILF, INT.NNANGC, At , XMAX) 


o o o Ti r> o c n o w> o r> o o r> r> r> onN>n» nnnnoon net net r>nr»nnr>nr> 


ORIGINAL PAGE 
OF POOR QUALITY 


dXCK««,FOH 


AtJBMOUTXNe 9UCbnR»X*iNPWt.9ef tCr,t(iT»XM*8,«t) 


VtBTUAU At(t) 


MNO T^e UM1T9 OP r* e M«i.MLoat «no the maincobe 
9 UM THE ANUUb aiice 


HtWINO 3 
K«ia - i 

IPIK.CT.I) SOTO 2 
00 t J«I|K 
Me*o(3) 

CONTINUE 


HE*0(S) (A1 (n,I«t|NHUUSE! 


uo axgmt side opt the peak 


KU«XA 
KAiXA * 1 

SUHHae, 

CONTINUE 

lA tNC.GT.NPUU9S) KC«1 
t)> (KP.GT.NPUL9E) HNat 

lP(AX(R»).Gb AX(KU)) GOTO 4 
9UM4«SUHI» • AX(KU) 

RL««l. ♦ X 
RHaRK ♦ I 
GOTO S 
CONTINUE 

XDTaRU 

EUNHa&wHA « AX(KL) 


PINO The wept 9I0E POINT 


e 


RRaXA 
RLalA - X 
9UNt*0. 
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9 

CUNTiNUd 

c 

t»'(K9.UT.n n4aMPU(.Se 
lf(R(.t(.T.U »i.«^iPUL88 

c 

If (U (KW).Gb.fUdP)) GOTO « 

c 

SUMV.«9imu « *l(KR) 

KOaKM m \ 

RtRRU - ) 

GOTO 9 

c 


c 


t> 

CONTlNue 

c 

ILfaKn 


SU**WaSU^L * At (HP) 


SUMaZisuMt « auMP > AlUA) 

c 

c 

C f !N0 

e 

Tne sxoetoae 

c 

TOTtaO, 

00 ifit I*l«NPUUSe 


itaAUX) 

207 

TOTC-TOTU ♦ «» 

£ 


C 


C 


C 

TOTSWRO. 

PtRKStaO, 

IPCAAa# 

C 


C 


c 

If (ILf .GT.XRT) GOTO 209 
If ULf .LT.2) GOTO 209 


Ka;LF - 1 


00 20a Ia(,K 

c 

TOT8L«TOT31. ♦ Aim 
If (Al n),LT,P6AHSi.i GOTO 


P&;.A8C«RUn 


XPfcARal 

c 


2oa 

CQNTXNUI 

c 


2«9 

roNTXNue 


e 

♦ 1 
C 

kkinpulsK 

If {Ulf.ST.lRT) KRrIUF *l 
00 I •«,«!( 

TOTSURTOrst. ♦ *tm 
If (»nn .uT.PfefNU) GOTO 2ta 
RfcfKGUaAl (!) 

XPt*R«I 


c 

2 ie coMTiNue 


c 

c 

c 

c 
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c 

c, 

c 


R&S*rOTL • SONiZ 


3 F rcc-fs 


c 

c 

c 


I 

t 

c 

eio 

t 

n 

5 

a 

!3 

C 

C 

c 


c 

I 

t 

iid 

i 

3 

i 

4 

c 

c 

c 

c 


OM 440 , 

IP(SUI^«Z.GT.0.) Oa*« >M. *41061 0(NZS/SUM*2) 

0di*4si), 

IF («M4'«.6T,0.} Ot5P4*|0.*«La6tO(HtS/XM4X) 


M»irE(6i0t0) ILF. X4,lNr«SUML»SUnit»SUMAZ,rorL» 

RES»Ott 4 (O 0 P 4 

MNITE(Si 0 tt>) XCF» I 4 |It<r,SOMi.,SUM 4 f SUMAZ,rorL» 

Hb3,0H4,0aPA 

F0PM4Tf» M4|,vL0Ue Of 4N«Lt tEFT POINT *,Vi/ 

* HIU POINT ',19/* NHT POINT SOM LEFT StPblS.S/ 

* SOM SihT '.IPElS.e/* M»jNl.O»e SOM *,lPei5v8/ 

' TOT*t SOM '»tPfclb,8/* SXOtLOat 80M 'ilPElS.S/ 

* SlOELOSb TO MAlNUOdt ( 08 ) SXPE1S«8/ 

' SXOF.LOab TO PEA4 ( 08 ) %lPtl5,0//} 


OO*0e 

XFtSOMAZ.rtT.O.) Oa*Xe.*ACOttl0(rOT3L/SUMAZj 

odSL*e. 

IF (VMAI ,eT,0.} O83(.*l9t*AUO&10(PEAitSt/XMAX} 

WHITE (8# 288} SOMAZ(TOrSl,PEAKSl.,IPEAK,XMAX|XA> 

OBtoasu 

WNirb(S>280) 80MAZ,T0}SL«PEAKSL,IPFAK,XMAX,XA« 
08(08SL 

F0«M4T(//» ANtfUt SUM ',lPtl9,8,' MAJN SOM *.IPE1S,8/ 

' P644 Sloe tOdt *»lPtlS.d,' COCATXON »,X6/ 
» maximum value MPElS.a,* LOCATION »,X6/ 

♦ iNTfeOMATEO (OBJ »,1PEI3,8/ 

» PSL/PEAK tOc!) *«tPEXS,8F/} 


RETOPN 

CNO 
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aicKaa.POM 


SUttSOUTIME NeOHn(NHUUSt» XA|N»itNC&,*t, A8» A3) 


VtMTUAt At U).A8(8e,8aa),AS(2O,800) 


REUi^GANtlE the AlLt IN AZIMJTN VEMitUa RANGi 
RATHfeR TmAN HANCe AllHUTH 

Pur The PEAR IN THE MIOOLE OP THE PILE IN AZIMUTH DIRECTION 


REHINO 4 

IMIOaNPUUSE/8 
ISTARTalA « IMIO 

IP(ISTART.UT.I) ISTARfatSTART ♦ NPUUSC 

NPASaaaeo 

NPAS 8 ANPA 9 S /2 

NPAa 28 «NPA 92 

NPaNPULSE/NPASS 
NQVbRaNpuUSE •> NpA 99 aNP 
ir^ (np.lt.i) GOTO daa 


00 6HB ja|,NP 


c 

c 

c 

c 


c 

BUS 

C 

c 

c 

cc 

c 


c 


REWIND 3 

xaaisTART 

00 Sto lat,NRANGE 

REAO(S) (A| (R)|Ral,NPUkSE) 
ISRaxS 

DO aOB RatiNPASZ 

IPdSR.GT.NPUUSE) ISRat 
A£(I.R)aAt(lSH) 
ISRaXBR * 1 

CONTINUE 


IP (NPAS22.UT.1) GOTO BI0 


00 PSA Ra|,NPA922 

IP(I9H,6T.NPUL9E) XSRal 
AACI,R}aAi (XSR) 
I 9 Hal 9 R • t 
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9Ui> CONM'OUe 

c 

•10 CONTlNUk 

C 

e 

00 BIO 

• is »»aTb(it) («in.K)»lil«N4*N0E) 

e 

c 

(ioro sir 
c 

00 916 Ml«NP*S<2 

• 16 NPtT6(«) («S(t»rtJ,l«t,NII«N(ie) 

e 

c 

81? I9T6flT«I9R 

C 

C 

880 CONTINUE 

C 

C 

8S8 IMNOVGR.lE.0) SOTO IS« 

e 

NR«l 

IF(N0VeR.Sr,NPA98) SOTO US 

c 

NP* 88 aN 0 VGR 
NP« 988 a 8 
NQVER *0 
SOTO I 

US CONTINUE 

C 

NPtS82aN0VeR - Nr«98 
NOV 6 R «0 
SOTO I 

|a0 CONTINUE 

C 

e 

RETURN 

ENO 
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SUdHOUTlNE V»t.ue(N(>UtSe» IHANGEi I*, IST* nRANGI, AU XMAI) 


VIRTU*!. *1(1) 


RfcMlNO « 


IM 10 «NPUC 3 E/a 
IDbLT**)* « ttP 

!>’ UatlTA.UT.Q) 10&LT*al0eUT* * NPUtSC 
UlML»ir*lO - lOeUTA ♦ I 
lO(5l.T*»lST - I* 

!*■ (UtUTA.ur.iS) IOEl.T*alOEl.T* ♦ NPUtSE 
UIMNilHlO • lOELT* ♦ I 

IPHRauiMH * 200 

IPRLaUWU - aO 0 
ir {iP"U,l.T,l) IPRUai 
!► UPRR.GT .nPULSE) IPR«*NPUL$E 


toRtTECft.TaRR) IMIO.MPUUSE, lUE!.T*,U(Mi.*UIMP 

v*uut '» 6 (lS,n)//) 


TOTACae, 


TOTSaO, 

TOTUae. 

vstotae, 

lVSqa 0 

lVS*a0 

DO 12B JaiiNPULSE 

PEAU(4) (*lin»lalfNR*NGE) 


TAG A PAlNCOae WITH XT a ) 


lTa 0 

IPCU.bE.UIHU .AND. 
t U.lE.U'^'O} ITai 
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C 


c 


Its 

c 

c 

c 


c 


113 

C 

C 


C 

Ita 

c 

C 


US 

c 


c 

US 

c 

c 

c 

c 

c 

c 


u 


00 tu ItliNHANGE 

T0T«li*l(X} « TOTAL 
CONTlNUt 


suMT«e. 

S 0 f*LL» 8 , 

(IT.E'J.O) goto ua 
IV«|f(«NGE 
SUMtaO. 

IVL«t V 
IWtialV * I 
CONTINUE 


ir (IVR.aT.NKAMGC) GOTO 114 
IF(Ai(Ivni.UT.AlUVU) GOTO 114 
3UML«!»UNL o AKIVL) 

IVfTtltfR t 1 
1VL*1VL * 1 
GOTO 113 

CONTINUE 

SUNCaSUML * AlCiVL) 

IF (I VQ.GT.NQaNGE) SUNLaSUNL * AUNNANfiE) 
IVNN*X VL 
SUNRtld, 

IVL*IV • 1 
IV4»IV 

Continue 

IFa^^L.LT.l} SOTO 114 
IF (&UIVL) ,6T,A1 GOTO 116 

SUMRaSUMN a Al(IVR) 

JVNalVN > 1 
IVUaiVL - 1 
GOTO US 

CONTINUE 

SUNRaSuNR a Al(IVN) 

IMIVL.LT.IJ SUMSaSUMR a 41(1) 
iVLUalVH 

SUMTasUMU * SUNN - Al(XV) 

IFUVLL.LE.n GOTO IS 

IVtalVLL - I 
DO Its laxUVL 

SUNtLaSUNLL a A1(I) 

IF (Al (1 ) .lT.vSXOE) goto 10 

VSIOEiAl (1) 

IVSRal 

IVSAaJ 

CONTINUE 
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c 

la CONTINUE 

c 

ir (IVHtt.n&.MHANGI') GOTO 20 

c 

IVRalVRH * I 
00 IT IiIV*«,NWAN6e 

c 

9UMtL«SUMtL * Aid) 

(AU11.LT.V8XO&) GOTO 17 
V!»]ue«*l (X) 

iVSMal 

IVSAaJ 

C 

IT CONTXNUb 

c 

c 

20 LOnTiNUft 

GOTO !!■? 

C 

iia coNTXNue 

c 

c 

00 50 1«I«NG*NCE 
C 

SunUL«SUMUL ♦ Aid) 

XFUl (X) .LT.VSlOt) ttOrO 10 
C 

wSlOfcaAt d) 

IVSAaJ 

C 

50 r ONTifoue 

c 

c 

110 continue 

c 

c 

TOTSaTOTS ♦ SUMT 
TOTU»TOTU ♦ SUHtL 
C 
C 

C PRINT f/UT The MAXNLOOe AREA 

c 

c 

ir((JAT.XPKU).OR.(J.GT,X»‘0R)} GOTO 92 

c 

00 99 X'liNRANCe 

e 

VACa«999, 

XIAAI d) 

XXaXV/XMAX 

If (XX.GT.O.) VALA-10.aa|.OGie(XX) 

A1 d)av»|. 

C 

89 CONTINUE 

C 

MRXTEll.dOSl) J»(Ald)datiNRAN6e) 

8081 PURN«Td6« 10(1X,F7«2)} 

C 

92 continue 

C 

C 

C 

C 

120 CONTINUE 
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C 


c 

c 

c 

c 


6000 

t 

i 

i 

Q 

9 

C 

c 

c 

c 


c 


e 

0062 

1 

i 

C 

c 

c 

c 

e 


:F(Tf}Tc,i>r.e.} otitt0.*«Loet0(ToTc/TOTS) 
ir (roru.sT.a.} oe»*at0,*«i.O6te(iroTt./iM«v] 
TOTSt'TOTO ♦ TOTL 


>^NirE(5>«0d0) T0TiU»r0TSC»r0T9»TQT|.,08f0l«P 
w*4lTfc(6f B0O0) TOT*C» r0TbL«TOrSiTOrt*UO»OBP 
fOhHtU//* rUTAk in KANGt AND ANGLO *«1P01S.0/ 
* total In mango anu angLO OUa *»|P0l9,0/ 

' total Of SOM 'ilMOtO.O/ 

' total OA 0X00 StPOtO.O/ 

' SlOt/OUM (oo; SXPE19.0/ 

' SIOO/MOAK (06} SlPOtO.O//) 


vsoaavsioe/xNAi 

XPtVOUa.GT.e*) VSoeal0,aALOGi0(V9O8) 

bRXTE(S»6o02) TVSR, (VSA, VOXOSiVODO 
MRXTGCb.OAOa] XVSP» XVOA, VSXOO«VSOO 

PU»MAT(//» PEAK oxoe LOOO LEVEL '/ 

'LOCATION IN MANGO SX<t«' LOCATION IN AZI ',10/ 

• MAGNiTuuff ',iHex 5 , 6 #* ( OB j sipexs.o//) 


retukn 

END 
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c* 

C* F»-T.FO« 

G* 

C* 

C* FFT •tSOTRlTNM 

e* 

c« 

8ubMouTx^b FFr(NUM0M,)ir»vr) 

OXMbNSION xrCDfVTd) 

NUMan IS the no. of points of uooo o«r* 

IF NUT A PONbH of rwo PAO TO NbXT HIGHEST 
WITH Zeros 


NtfPOwat 
I f E hp«NUMbR 

I XTEHPaiTtMP/S 

IFdTEHp.UE.n GOTO lU 
NiiPOWaNSPOW a I 
GOTO I 

te CONTINUE 

C 

NaSaaNSPOW 

IF (n^ge.nuhbr) goto is 

NaNaS 

N^POwaN^pQW a I 

IS continue 

IF (nupep.gc.n} goto as 

G 

NiaNUMBP a I 
OO 19 (.**«n 1,N 
C 

UTlUHJaB, 

TTCtMJaB, 

C 

19 CONTINUE 

C 

e 

ao continue 

c 

c 

HaNgPCW 

00 6(919 teat,H 

U'4laaa«(H*CI9) 

LlVaaaUNX 

SCt*6.aB3ie9/Ft0AT(UX} 
00 bl 6 <e LHaI,LN« 
aw6aCLH-l)asCL 
CaC0S(*H6} 

SaSlN(tPG) 

OO 6(9(9 UaUIX.N.LIX 
Jl»l.r-CI*aLH 
JaaJl^tHX 
TiakTUl)-«T(ja) 

. TdaYT{Jl)-TT(Jaj 
*T(JlJaiT(jn*XT(ja) 
TT(Jt)aYT(Jn*TT(jaj 
IT (ja) aCaTltS«ra 
VTUajaCaTa-SaTI 
600 continue 

NvaaN/2 

NHlaN*t 

J*I 
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00 

so TO bSl 

Tl aiT (J) 

T<*»T(4) 

If (J)«xr(|) 
»T(j)»yT(iJ 
xT(t)*ri 
TT(I)iTil 
6St KIMV; 

ft2i9 CUNTliije 

lX(K.bi\J) 60 TO 639 
JiJ-K 
H«K /2 
60 TO 620 
6i9 JaJ*A 

NUMbRaN 
At Turn 
eNO 
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FACIE I'- 
QF poem Ci;AL!l“Y 


WT,^0» 


MXN00WIN6 FUNCTION 


AR&UMeNTO 


*HR»y 19 THe ARRAV TO STORE THE HEIGHTS XN 
NUMtt X5 THE NO, OF SAMFtiNG EtRNfeNTS IN im?l 


XOPT IS TYPE OF WINDOW 


^ NO uiiNOowiNu ( REcr Function ) 

1 harming UXNOOW 
e 9ARr(.6TT WINDOW (TRIANGLE) 

3 manning WINUOW 
<t SLACKMAN WINDOW 
3 23 09 TaYLOK wiNOOW 

6 3» 03 TAYlON window 

7 33 03 TAYLON WINDOW 


SUOROUTINC wr (ARRAY, NUWaP, lOPT) 


c 

c 


c 


c 


c 


C 

10 

c 

c 

C 


REAL ARNAV(I) 


IF (XOPT. EG. 0) GOTO 600 
ININUN3R " t 
IF(10PT.EQ.2) GOTO 200 

PI«a,AATAN(l,) 
PI2A2.API 
PmA8EaPI2/AN 
IF (lOPT.EO.tt) GOTO 100 
IF nOPT,EQ,6) GOTO 300 
IF (tOFT.EQ.E) GOTO 400 
lF(IOPT.EO,n GOTO 300 


C1A0.9 

C2»«.3 

IF (I0PT.EO.3) GOTO 10 

CIAO.54 

C2*0,46 

CONTINUE 


OU 20 Xal,NUM0R 


ANSI • I 


A70 


e 

e 


c 

iMt 

C 

e 

c 


c 

u« 

c 

c 

c 

c 

c 

c 

e 

c 

c 


c 

■m 


C 

e 


c 

Ud 

c 

c 

}»« 

e 

e 

C «9 

e 

c 

c 

c 


c 


c 

sie 

c 


C’.» I 


'' im fjn. 




w(;uCt - Cii^COACHASlKKN) 

CUNriNUK 

GOTO 149(99 
COMYINUi 

PM* 9 l 4 «f>N«Se* 4 , 

00 129 ia(«NUnSN 
XN»| - I 

«Ua9,«2 •> 9.99*C09(Pr(*.le*XIM) * 9t94*C0S(^NAa£2*XN) 
AKHtyUJiftO 

CUNTINUK 

GOTO 19009 

NUMrtA 18 TNe NUHoen or ettntNrs in tnc minoow staktinq 
rMOM TMfc riNST AOOMtss or an 9 at 


CONTINU& 

lOI»aNUN 99/2 

00 21(9 XalilUK 

«<iUirL0*r(8«(I • U)/XN 
ANNAV(I) *W0 

CUNMNUft 

lUPalur * t 

00 2 i» l«I’JP,N(JNS 9 

wo» 2 , - rtoAT(a«u - uj/xN 

AHHAV(|)tMO 

CONTINUE 

uoro 10999 
CONTINUG 

09 YAV(.09 lnelGHTINa 


8 r*«T»-Pi 
00 319 I«l»NUNa9 

AYyaCOSCSrAKT) - 0 . 908 3 '0*COS ( 2 . *ST ARY ) > 
9.9996ai*C0S(3.*8YART9 * 0.99U*l2r*C0S(a,tSrART) 

XRPAYlDat, ♦ a.aMO 

STiHYaSTAWr * PhASC 

CONYlNUt 
GOTO 10999 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


ya CoNTiNub 


id Ue TAVLOn WtlGHTlNG 
OU 

MU«ut.tf4^b>)6*CU9(ST*RT) - e«igi9TB3a*C0S(ii.*ST*NTi ♦ 

ANHAY({]«1, * 2«M0 

9T*NT«ST*nt « PHA9& 

IB continue 

GOTO IMBae 
UB CONTlNOft 


J9 09 T»VtON weiOHTlNG 


9T4HT*-Pl 
00 91B lalfNUMSR 

»lO«B.5ttU33B*CU9(ST*HT} • B ,Bt Sl<)«^?*C09 (8 . *9? AHT) ♦ 

I Q.B0a<*Tan«CU3(i.*STAHT) - B . BOB T 3'«9S 1 *C09 (tt . «9T ANT] 

AMNAVinalt * S*N0 
iiTANTaSTANT ♦ PHASE 

CONTlNUt 


GOTO IBBBB 
MB CONTINUft 


ONIPUNH weiGHTINQ 


c 

c 

61B 

IBBBB 

C 

c 


00 MB lalfNUMSN 
AHNAVdJai, 

CONTINUE 

continue 


RETURN 

ENO 
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